水下机器人毕业论文

1、目录1 绪论 . 11.1 研究意义 11.2 AUV介绍21.3国内外AUV研究动态21.4 动态仿真研究现状 31.5 本文主要研究内容 42 AUV总体设计和三维建模 62.1 形体的选择 62.2 设计内容 62.3 三维建模 73 AUV的动力学分析103.1 坐标系 103.2 定义运动参数 103.3 受力分析 114 推进器动态仿真 164.1 ADAM仿真164.2 FLUENT仿真204.2.1 理论基础 204.2.2 仿真前期准备 214.2.3 FLUENT 数值模拟计算 235 总结 29参考文献 30致谢 31附图 321绪论1.1研究意义今天的人类正面临着人口、

2、资源和环境三大难题。随着各国经济的飞速发展和世界 人口的不断增加，人类消耗的自然资源越来越多， 陆地上的资源正日益减少。为了生存 和发展，人们必须寻找新的物质来源，海洋应当是首选。海洋是一个巨大的资源宝库，开发蓝色国土，拓展生存和发展空间，“人类重返海洋”将成为全球经济发展的大趋势。世界沿海国家和地区正在进入全面开发利用海洋的 新时期，美、俄、中、英、法、日、加拿大、韩、印度以及东盟诸国等140多个国家相继制定海洋科技发展和海洋开发计划， 采取具体措施加快抢占海洋科技的制高点， 海洋 开发已成为全球产业进步的重要标志，海洋经济已成为全球经济发展的重要增长点。海 洋是强国之本。谁掌握了海洋，谁就

3、掌握了经济发展的未来。没有强大的海洋科技事业， 没有强大的现代海洋经济，就不可能成为真正的经济强国。我国是海洋大国，但不是强国，人均占有陆地面积和资源量都远远低于世界平均水平，研究开发利用海洋是顺应世界海洋开发大潮。为了推动海洋经济持续快速的发展， 科技部、国家计委、国家海洋局、 农业部联合推出了 “科技兴海”计划。中国海洋21世纪议程把“科教兴海”作为海 洋经济可持续发展的重大意义的战略选择。 海洋经济的开放性和带动性，强渗透力、宽 辐射，可以拓展新的经济增长空间，优化经济结构，建立开放型特色经济，促进我国更 快走向世界，对实现中华民族的伟大复兴具有重要政治经济意义。21世纪是人类向海洋进军

4、的世纪。深海作为人类尚未开发的宝地和高技术领域之 一,已经成为各国的重要战略目标，也是近几年国际上激烈竞争的焦点之一。 水下机器 人作为一种高技术手段在海洋开发和利用领域的重要性不亚于宇宙火箭在探索宇宙空 间中的作用。然而水下机器人如果直接进行海试或湖试， 则需要承担极大地风险，由此计算机仿 真应运而生。仿真系统提供一种有效的试验来检验系统和设计，这种试验可以为设计者、 使用者和购买者揭示和预测一些有价值的信息。仿真系统有助于使各系统和子系统之间 的协调。另外，仿真系统可以模拟各种场景，用于各种操作训练，节约实地训练的费用， 最大限度地降低训练风险。通过 ADAMS/MATLAB立仿真模型，可

5、以实现对其控制路径 的模拟，降低试验风险。在海底工作的动力来源就是推进器，所以通过对推进器的流体仿真，研究其所受的的阻力以及叶片周围流体所受的压力变化，对于AUV的运动控制的精确性具有极大的指导意义。1.2 AUV介绍当前水下机器人的种类很多，其中载人潜器、有缆潜器（ROVs）和自治水下机器人（AUVs）是三类最重要的潜器，自治水下机器人AUVs是英语“自治水下潜器” （Autonomous Un derwater Vehicles）的缩写。AUVs不配备主缆和系缆，因此它又称为无人无缆水下机器人（Unmann edU ntethered Vechiles缩写UUVs。这类水下机器人携带能源，

6、依靠自身的自治能力来管理自己、控 制自己，以完成赋予它的使命，自治水下机器人也就因此得名。由于微电子技术、计算 机技术、人工智能技术、导航技术的飞速进展，再加上海洋工程和军事活动的需要，国 外产业界和军方再次对无缆水下机器人发生了兴趣。许多研究表明，无缆水下机器人是 一种非常适合于海底搜索、调查、识别和打捞作业的既经济又安全的工具。与载人潜水 器相比较，它具有安全（无人）、结构简单、重量轻、尺寸小、造价低等优点。而与 ROVs 相比，它具有活动范围大、潜水深度深、不怕电缆缠绕、可进入复杂结构中、不需要庞 大水面支持、占用甲板面积小和成本低等优点111。1.3国内外AUV研究动态海洋机器人在过去

7、几十年间为世界各国的海军、石油开发和救援打捞开辟了崭新的 活动领域。用这种高度计算机化、有的已开始了不需要人进行控制的机器人进行探索海 底，可提供关于海底的大量数据。当今世界各国的一些主要的海洋研究中心倾注很大的 精力正在研制或使用数十种可深潜的海底机器人。 在技术方面，美国的水平领先于世界， 欧洲各国其次，而日本要落后于美国和欧洲。这些机器人的造价与载人潜水器相比，造价低得多，但更加安全，而且可长时间在压力很大的海底工作。从上世纪90年代中期以来，自主式水下航行器（AUV）在海洋科学调查以及军事领域得到越来越广泛的应用。 截止到2005年，世界上共研制了约70艘AUV AUV已经多次成功地应

8、用于海底石油与 天然气、天然气水合物、大洋多金属结核和热液硫化物矿床等海洋矿产资源的探测。AUV的最新应用主要体现在水雷搜索，冰下探测以及水下作业中。目前，美国在开发工作方面正朝着无缆预编程式（即自主式）海底机器人的方向发 展。在过去的20年里，全球各主要国家所研制的 AUV其中美国占了一半。就其性能 而言，应数美国海军的水下搜索系统(AUSS和新罕布什尔州大学的小型AUV即试验性 的自主式EAVEEAST主要用于水下管道和平台的检查。美国AUV的研究和开发经费90%来自国防部的财政支持，主要通过国防高技术研究 项目局(DARPA和海军有关部门给予拨款。早在1988年，国防高技术研究项目局提出

9、了 开发既像无人潜水器UUV那样自由浮游，又像登陆艇那样能离开水域进入内陆的两栖性 自主式机器人。美国海军控制和海洋监测中心于1983年推出了先进的自主式无人搜索系统(AUSS)，本系统的外形尺寸为510cmK 75cm,有效负荷不定。系统用银锌电池作为 动力源，系统配置了侧扫声纳、前视声纳、35mm静物照相机、CCD水下摄像机。导航系 统由多普勒声纳、先进的精密陀螺仪和周期性更新的水声通道组成；水面母船借助LBL发射应答网执行跟踪任务。现在 AUSS为第二代先进的作业型自主式海底机器人，潜航 深度6000m用于评价深海海底搜索技术和执行深海海底搜索作业使命。近年来我国政府十分重视 AUV的发

10、展，投入了较大的人力和财力，在20世纪90年 代我国AUV的研制取得了重大突破，典型代表有“探索者”号1000米无人无缆遥控潜水器和“ CR-01A 6000米无人无缆遥控潜水器。自主式无缆海底机器人是国家863高技术发展计划支持项目，由中国科学院沈阳自动化研究所、声学研究所、中国船舶工业总公司702所、哈尔滨工程大学、上海交通大学等单位联合设计研制的“探索者”号自 主式无缆海底机器人，是以大范围搜索、观察水下1000m或6000m失事目标为主要使命 特征的。我国在军用AUV的研究中也取得了很好的成绩，研制出三型军用智能水下航行器。 其试验平台“智水II号”于1995年夏进行了湖试，在自主导航

11、、自主避障和自主简单 作业等方面取得了成功的试验验证和宝贵的试验数据。此外，北京航空航天大学目前正致力于仿生机器鱼的研究，已研制出多种形式的仿生机器鱼，为水下航行器新型推进器技术和新型结构的研究，奠定了良好的基础2。1.4动态仿真研究现状机械系统动力学仿真分析技术首次出现于1980年前后，作为一门新兴技术，它最初被应用在汽车、铁路等领域中。Woongsan命以提高汽车的稳定性和控制能力为目标， 进行汽车四轮定位系统研究。由于此前研究都是采用简化模型，使用二维结构或以自行车代替，测得的数据很不可靠。后来随着分析手段的提高，动力学仿真分析技术开始大 量地应用于空间科学、石油、机器人等领域，NOEL

12、通过在动力学分析软件中建立、分析和优化模型，得到了飞机起落架的动态性能。Arenz等针对goliath移动机器人模型， 利用ADAMBANSY别MATLA三者联合进行了动力学分析，并针对goliath移动机器人 控制算法进行研究，然后在动力学仿真软件中加以检验。随着水下机器人应用领域的不断的扩展，水下机器人正朝着复杂化、多功能化方向 发展。AUV正是一个由多个系统组成的复杂大系统，它包括载体系统、通讯系统、导航 系统、传感器系统、声纳系统、能源系统、推进系统和智能控制系统。这不仅增加了AUV各系统之间集成的难度，也使得 AUV整体性能验证和评价成为一件困难的事情。由 于AUV勺复杂性、智能控制

13、的局限性和环境感知的低能，因此AUV可靠性、稳定性和自主能力在没有得到确切的验证之前，对AUV进行海上实验是有很大的危险性。另一方面, 由于海洋环境复杂、危险，在海上进行水下实验，除了没有安全性之外，还有实验成本 高，周期长的缺点。因此，AUV整体性能验证和评价已经成为 AUV发展的瓶颈，而建立 AUV实时仿真系统是解决这个瓶颈问题的有效途径。控制系统仿真是控制技术人员在开发控制系统时所经常使用的。通常基于简化的被控系统模型，构建控制系统，设计控制算法，然后利用仿真分析被控系统的性能（如导 弹的反应时间、跟踪精度等）412。1.5本文主要研究内容借鉴国内外较为先进的AUV采用母型设计法，对其外

14、形和部分机构进行改进，使其在达到总体性能要求的同时，尽量减小其体积和降低生产成本。在用SolidWorks完成三维造型设计之后，在ADAM中进行仿真分析，建立简单的控制系统，然后在利用流 体仿真软件FLUENT寸推进器进行流体分析，计算不同转速下推力情况，具体内容如下： 第一章绪论本章论述了研究、开发、利用和保护海洋的重要性和迫切性，概述了水下航行器的 发展和应用状况，分析了国内外 AUV勺研究现状与趋势以及现阶段对动态仿真的研究， 明确了本文研究的任务、意义和主要内容。第二章AUV总体设计和三维建模本章主要是论述了 AUV的总体设计及形体选择原则，根据二维图纸进行三维建模，对主要零件进行功能

15、分析，并简要阐述了所研究对象的主要参数及工作环境。第三章AUV的动力学分析本章为研究机器人总体的运动分别建立了地面坐标系和体坐标系，定义了各运动参数。对海底机器人的运动过程，给出了其各方向的受到的合力及力矩方程，并给出了其 轨迹方程，由于部分参数需经过海试得出，故此只给出计算公式。第四章AUV动态仿真本章是本文的重点，主要是从两个方面对推进器进行动态仿真，一方面是基于 ADAMS/MATLAB行动态仿真，建立控制模型，另一方面则是利用 FLUENT寸推进器进行 流体仿真，分析不同转速时叶片所受到的压力情况以及所产生的推力变化。第五章总结2 AUV总体设计和三维建模水下机器人是一种具有人的一部分

16、或大部分功能，能够在海洋环境下代替人进行某 种作业的自动控制装置。通常是依靠电缆提供的动力(有缆水下机器人)或自身携带的能 源(无缆水下机器人)驱动载体上装有的推进器在水下作三维运动， 并可装设照明、摄像、 声纳等观测设备，多种传感器及用来完成某种作业的机械手和作业工具。由于海洋开发所需从事的项目内容非常多，实施的作业也极富多样性，因此，水下机器人的系统构成、形体、总体布置都不尽相同。虽然有许多国家都在研制和生产水下 机器人，但每种水下机器人都没有进行批量生产，因此，目前还没有一个完善的设计准 则。但在水下机器人结构方面己有一些成功的经验， 可供设计水下机器人结构时作参考。水下机器人外形结构设

17、计所涉及的因素很多，主要包括以下几个方面：水下机器人形体的选择、耐压壳体的设计、推进器的布置、水密接插件、重心和浮心的计算和整体 平衡的考虑。2.1形体的选择AUV外形设计是AUV总体综合优化设计的重要组成部分，并直接关系到 AUV勺诸多 战术技术性能。AUV外形设计技术是研究开发新型 AUV勺必不可少的一项基础技术5。 水下机器人根据使用目的和技术要求的不同，其外形尺寸、结构形式都有很大差异。 一 般情况下，形体的选择要考虑以下原则和要求：(1) 阻力小、航行性能好；(2) 足够的强度；(3) 便于总体布置；(4) 良好的工艺性。2.2设计内容本文所研究的水下机器人驱动装置来源于 5个推进器

18、，推进器是由24V、400W的直流电机带动，由控制系统对五个推进器进行转速调整以控制机器人的上升、下潜、旋转等功能，实现机器人的姿态调整。高压舱内安装有探测及发射装置，控制系统等设备。本课题所设计的对象包括两个高压舱、五个推进器、两个声纳(一个大声纳、一个小声纳）以及起支撑作用的框架该AUV在海水中航行时为中性，即重力与浮力相等且主体的轴线呈水平， 海底机器 人航行到达预定水域后，推进器停止工作，在完成工作任务后，又通过推进器的作用完 成返程。2.3三维建模三维建模的基础就是二维图纸的设计，三维建模成功与否在很大程度上决定于零件 和装配的设计精度，所以在完成设计工作时，应该在满足强度要求的前提

19、下，根据装配 关系进行零件设计。完成整体分析后，先由AutoCAD完成了设计的装配图、零件图，确定各部分尺寸， 然后开始用SolidWorks建模。所设计的三维图主要包括：（1）高压舱（如图2.1）主要由咼压舱盖、舱体、咼压球舱盖三部分组成，共有2个，每个咼压舱由两个卡箍固定于框架之上。材料选用 6063合金，该材料主要合金元素是镁和硅，具有极佳的 加工性能，优良的可焊性、电镀性，良好的抗腐蚀性、韧性，阳极氧化效果优良。图2.1高压舱（2）框架（如图2.2 ）支撑作用，上面固定有高压舱、推进器、声纳等，是由35X 2的钢管焊接而成,上面有8处配焊，作用是固定卡箍。框架要先进行密封性试验，在海底

20、工作时，钢管 内部充满氮气，抵消部分海水压力，防止其在海底工作时发生变形。图22 AUV框架(3) 推进器推进器作为整个AUV动力来源，共有两种，一种是水平推进器，共两个，另一种是 竖直推进器，共三个，通过推进器固定架将其固定于框架之上，各个推进器独立工作， 相互协作实现对整个AUV姿态的控制。推进器上固定有导流罩，主要作用就是在运动中 对迎面而来的海水阻力进行分割化解，从而达到降低阻力的效果。(4) 声纳声纳主要是利用水中声波对水下目标进行探测、定位和通信的电子设备，是水声学中应用最广泛、最重要的一种装置。 它是该海底机器人主要的探测装备，包括大小声纳 各一个，分别安装在AUV前端的上下部位

21、。所设计的三维模型如图2.3所示：图2.3 AUV总装图大声纳；框架；竖直推进器；高压舱；小声纳；一一卡箍；一一水平推进器此AUV重 208kg，总排水量206L,长1500mm宽1240mm高920mm工作深度为800m压强约为8MP温度约为2C,正常工作时的速度2m/s，主要用于深海探测。3 AUV的动力学分析对5自由度AUV勺动力学与运动控制进行研究，应该考虑重力、 浮力、推力以及水 动力的影响，建立水下机器人的动力学模型，对机器人的复杂水下动力学行为进行描述。研究水动力的意义：一是从操纵性的角度研究水下机器人载体的稳定性和快速型， 另一方面在设计控制系统时需要考虑水动力的影响，以便建立

22、AUV勺数学模型。3.1坐标系为便于研究AUV勺运动、数学建模与控制，需建立两个正交坐标系。先建立地面坐 标系E- EnZ作为静坐标系，E为海平面一点，以海平面出发点为原点，EE、En为水平面如图3.1所示：起，Gx与其轴向一致，Gy指向右舷方向，如图3.1所示。本课题所研究的对象在海底的运动是 5自由度的空间运动，包括沿x、z轴的平动 和绕x、y、z轴的转动，由5个推进器来实现。3.2定义运动参数V 机器人重心G的空间速度矢量V在水平面内的投影；V 机器人重心G的空间速度矢量V在纵垂直面内的投影；u、v V在水平面内沿Gx Gy方向的分量；u、wV在Gx、Gz方向的分量；p、q、r沿Gx G

23、y、Gz轴的角速度；© 横倾角，Gy轴和En轴间的夹角；9纵倾角，Gx轴和EE轴在纵垂直面的夹角; 书艏向角，Gx轴和EE轴在水平面的夹角。3.3受力分析为了设计AUV的控制系统，首先需要建立机器人的动力学模型。 在水下运动的机器 人系统是一个非线性的动力学系统，需要确定的参数较多。由于技术和测试条件的限制， 有些参数无法准确测定或者无法测定。为了控制的需要，有必要对系统进行必要的简化， 而只考虑对系统性能起主要作用的影响因素，这里主要考虑重力、浮力、推力和水动力对机器人的影响。1.AUV在水中所受的合外力可用下式表示：nF Ff B PTii 1（注：该方程为矢量方程，为矢量和）（

24、3-1 ）其中：Ff 作用在AUV上的水动力；B AUV的 浮力；P AUV的 重力；nTii1 所有推进器推力之和；斤一第i个推进器的推力；n推力器的总数，该处为n=5。AUV所受的合外力矩方程：nM M F M B MP MTii 1（注：方程为矢量方程）（3-2 ）nMf、Mb、Mp、i /T'分别为水动力、浮力、重力以及推进器推力所产生的力矩。2.力学分析(1) 推力：Tn2D垂向力方程：KT(3-3)n转速 rad / sD螺旋桨直径(113mmKt推力系数(经验值)海水密度(1030 kg/m3)(2) 轴向力方程：2 2m uvrwq Xg qryg pq rZg pr

25、q12L42 2X qq qX rr rXrprp1L3X. u XvrvrXwqwq(3-4)2u12L222XuuUXvvvXwwW2TxPB sin等号左边描述了水下机器人合力作用下的运动规律，右边第一项和第三项分别是角 速度和速度引起的非线性水动力， 第二项是惯性水动力，第四项是推力，第五项是重力 和浮力。(3) 侧向力方程：2 2m v wp ur yg r pZg qr p Xg qp r1 L4(Y'rY'pYp|p|PpYpqPqYqr)2 r p1 3'.'-L3(Y v Yvqvq YwpWP Ywwr)2 vM41( 3-5)-L3(Yur

26、 Ypup &打 v2 w2 2 |r)1 1-L2(Yu2 乂uv Yv|vV v2 w2 2) Ty1 2 'L Y<wvw (P B)cos sin2uq vpzg pXgrpqYg rq p4'L (Z.qq' 2ZppP乙rr2Zrqrq)L3(Z'. wwZwvrZ；pvp)L3ZquqwZwq w(3-6)L2Z*u2w21iiUw1Zw w|1w v2w2 2|w|ZwuwB)cossinZv2L2Zw|w w1辽加2 (PTz(5)摇力矩方程、纵力矩方程、偏航力矩方程I XX p(Ixx Iyy)qr m Yg (w vp uq)

27、 Zg (v ur wp)L5K- r K. p Kp|p|P prpKpq pq Kqqr1212L4L4K. v Kpup KrurvKVqVqKwpwpKwrwr(3-7)L3'2K*uKvuv KvL3K；wVwphcos sinMTxlyyq(I xx I zz)rp m Zg (u wpvr)Xg(w vp uq)L5K. rK' p Kpp|P ppKpq pqKqrqrL4Mv,rMvpVpM quqM |w p(3-8)2 2 2(v w )2 qL3M*u21|wu wMwwWv2 w"M2L3MTyM wuw Mw|w|Wph sinIzz(Iy

28、y I xx) pq m Xg(v urwp) yg(u wq vr)1212L5L4N. rN. vvN. pNr|r|r rpNpq pqNqqrNvqVqNwpWPNwrwr(3-9)L4NrurNpupv2 w2L3N*u2N；uvNvv2w2L3NYwVWMTz式中TxTy Tz为推力在坐标轴y、z上的投影；Mtx MTy Mtz为x、y、z轴的推力矩。（6）根据运动学知识可以将水下机器人姿态角与运动坐标系中角速度的关系可以表示为:p qtgsinrtg cos（3-10）qcosr si n（3-11）(qsinr cos)/cos（3-12）水下机器人重心与动坐标系原点重合，水下

29、机器人的运动轨迹o ucos cos v(cos sin sin sin cos ) w(cos sin sin sin sin )o usin cos v(sin sin sin cos cos ) w(sin sin cos cos sin )usin vcos sin wcos cos(3-13)(3-14)(3-15)(3-16)(8)有海流（海流速度为U （ Ux Uy Uz ）时则水下机器人相对于海流的速度为:Vrx = u-uxVRy=v-uyVrz=W-uz若考虑海流作用时，用相对海流的速度代替前面公式中的速度即可。本章为研究机器人总体的运动分别建立了地面坐标系和体坐标系，定义

30、了各运动参数。对海底机器人的运动过程，给出了其各方向的受到的合力及力矩方程，并给出了其 轨迹方程，由于部分水动力参数需经过海试或湖试才可以得出，故此只给出计算公式。4推进器动态仿真AUV按使命要求在水下平稳地航行是其完成各种任务的基础，要想保证机器人平稳地航行，选择正确的控制方法是至关重要的。尽管现在有许多控制方法可以实现对潜器 的航行控制，但对于实际应用来说，在满足项目要求的情况下，应选取既可行又实用的本文在完成AUV勺三维设计后，用动态仿真软件ADAM对单个推进器进行动态仿真, 建立正确的控制系统。然后再利用 FLUENT对单个推进器进行仿真，实现对螺旋桨叶片 转速与推力的分析。推进器模型

31、如图4.1所示：图4.1推进器模型4.1 ADAMS仿真仿真的基本思路主要是针对单个推进器进行仿真， 验证控制系统的准确性，在ADAMS 中建立底层控制：AUV勺底层控制AUV勺控制系统包括底层控制和高层控制，高层控制是由计算机软件来完成，底层 控制主要由硬件组成，包括执行机构和传感器，这部分是其基础。水下机器人的底层控制主要包括13:(1) 水下机器人的航行控制；(2) 作业工具控制；(3) 传感器控制和信号控制。本文主要完成第一个问题即航行控制。其主要参数是深度、高度、航行速度、航向 角和位置等，为了简化计算和便于仿真，本课题将采用典型的PID调节器作为闭环系统 的控制器进行深度控制。该A

32、UV采用推进器作为动力和控制部件。1. 推进器布置如图:T12. 自动定深回路图体现了定深回路的控制框图:Di9Do图4.1.2深度控制框图图中Do为深度输出，Di为深度输入 典型的PID调节器闭环控制方程：plKd 評Do)(4-1 )tMyc Kp(Di Do) Ki(Di Do)d0Kp Ki Kd为调节器系数。3. 对单个推进器进行单自由度仿真(沿 y向的平动)(1)施加约束：重力，螺旋桨产生的推力(根据推进器的运动由 ADAM自动调节)，阻力(10vv );输入：推力（tl ）;输出：位移（weiyi ）、速度（sudu）;推进器的设定位置：距推进器位置 3m处（2）在MATLAB建

33、立PID控制图所示为控制图：图 4.1.3 MATLAB/Simulink 控制图时间/s图4.1.4位移曲线设定 P gain=200, I gain=25 ， D gain=100。通过ADAMS/MATLAB合仿真得到仿真曲线：（1）位移曲线（如图）稳定振幅0.02m :位 移/m（2）速度曲线（如图）稳定振幅0.3m/s :' 1 1 1 1速 度m/s1 1时间/S图速度曲线图此控制系统能使推进器在允许误差范围内到达预定位置，所以该PID控制是正确的，控制方法是实用的和有效的。 为了使控制效果更为理想，可以通过改变 PID系数的 值得到。采用单一变量法进行分析，如表:表不同P

34、ID参数对比第一组PID位移振幅（m）速度振幅（m/s）200252000.020.35200254000.020.3第二组2002002000.0360.582001002000.0280.4第三组400252000.0250.55100252000.020.2最初200251000.020.3通过减小微分参数可以减小位移和速度的振幅，增加比例环节可以适当的减小到达预定位置的位移，但通过曲线可以看出影响不是非常明显，所以通过比较选定PID的参数分别为100、25、200，生成的仿真曲线速度振幅 0.2m/s,位移振幅0.02m。4.2 FLUENT 仿真对单个推进器进行FLUENTS真，主要

35、是研究其叶片在给定转速时叶片上的压力分 布、叶片表面的速度矢量分布以及分析叶片在旋转时所产生的阻力等，对进一步研究 AUV勺海底控制有一定的指导意义。理论基础1. 流体流动要受物理守恒定律的支配，基本的守恒定律包括：质量守恒定律、动量守恒定律、能量守恒定律，如果流动处于湍流状态，系统还要遵守附加的湍流输运方程8。质量守恒方程:u 0t式中P是流体密度，t是时间，U速度矢量动量方程 du f pdt式中 微元表面粘性应力，p为流体微元体上的压力能量守恒方程：deP u q qdt式中e表示内能，q交换热量状态方程：f(P,q,T) 0内能公式：e e(p,T)傅立叶热传导公式：q K T ,式中

36、K为导热系数2. 静压、动压、总压1 2pgz gH(4-2)(4-3)(4-4)(4-5)(4-6)2静止的流体所具有的压强即为静压强，简称静压，-2是因流体运动而产生的压2强，称为动压强，即动压，gH为总压强，即总压9。关于FLUENT中所涉及的边界压力说明：PabsoluteP gaugeP operat ing等号左侧为绝对压强，右侧第一项为表压强，第二项为操作压强。本文的操作压强为8MP3. k_epsilo n（2eq n）模型标准的k_模型是最基本的二方程模型，k-湍动能，epsilon-耗散率动力学湍流粘度定义为:k2式中C为无量纲常数。标准k_模型k、 的运输方程8:kkUi

37、XjikGkkXjGbYmSk（4-7）tXiUiC1 GxjkC3 Gb )C22S k（4-8）tXiXj其中，Gk是由于平均速度梯度引起的湍动能k的产生项，Gb是由于浮力引起的湍 动能k的产生项，Ym代表可压湍流中脉动扩张的贡献，G C2 C3为经验常数，k 分 别是与湍动能k和耗散率 对应的普朗特数，SkS是用户定义的源项。标准k_模型中的五个可调常数值为：C =0.09， k=1.00，=1.30，G =1.44，C2 =1.92。对于可压流体的流动计算中与浮力有关的系数 C3，当主流方向与重力方向 平行时，取1,当主流方向与重力方向垂直时，取 0。仿真前期准备GAMBIT前处理一一

38、网格划分用CFD（计算流体动力学 Computational Fluid Dynamics，简称CFD方法进行流场计算时，首先要将计算区域离散化，即划分网格。网格是CFD模型的几何表达形式，也是模拟与分析的载体。计算网格的好坏直接影响到数值计算的可行性、收敛性以及计算精度。主要过程：首先将推进器三维模型导入到 GAMBI冲，建立计算域，选择FLUENT5/6 求解器进行网格划分，设置边界条件，输出网格。1. 建立流体区域：首先将SolidWorks生成的推进器模型导入GAMBIT中，根据推进器的高度和叶片的 直径建立半径为120mm高为60mnm勺圆柱体区域，通过布尔差运算建立最终的流体区 域

39、。2. 网格划分：设定网格单元为Tet/Hybrid (网格主要由四面体组成，个别位置可以有六面体、 锥体或楔形体)，选择TGrid混合网格类型。制定在边界上分点时所用的间隔长度为5,即 Interval size=5 。图所示为网格划分图：图网格划分图(spacing=5 )3. 设定边界类型(1) 选择FLUENT5/6求解器(2) 指定边界类型。叶片类型为WALL包括叶片前面，叶片后面和端面，名称分别为ypqian，yphou和ypduan，轮毂的类型为 WALL名称为lungu，进口面的类型为 pressurenlet，名称为in let ，出口面的类型为 pressure_outle

40、t ，名称为pressure_outlet 。423 FLUENT数值模拟计算1. FLUENT求解步骤10:（1）创建几何模型和网格模型（在 GAMBIT前处理软件中完成）；（2）启动FLUENT求解器；（3）导入网格模型；（4）检查网格模型是否存在问题；（5）选择求解器及运行环境；（6）决定计算模型，即是否考虑热交换，是否考虑粘性，是否存在多相等；（7）设置材料特性；（8）设置边界条件；（9）调整用于控制求解的有关参数；（10）初始化流场；（11）开始求解；（12）显示求解结果；（13）保存求解结果。2. 数值计算数值计算是将描述物理现象的偏微分方程在一定的网格系统内离散，用网格节点处的场

41、变量值近似描述微分方程中各项所表示的数学关系，按一定的物理定律或数学原理 构造与微分方程相关的离散代数方程组。引入边界条件后求解离散代数方程组，得到各 网格节点的场变量分布，用这一离散的场变量分布近似代替原微分方程的解析解。（1）计算模型选择分离式求解器（segregated solver ），该求解器主要顺序地、逐一地求解个方 程。（2）粘性模型选择k_epsilon （2eqn）模型，使用标准k- £双方程模型进行湍流计算，使用标 准的近壁函数。模型参数：Cmu=0.09， C1 epsilon 1.44， C2 epsilon 1.92，TKE prandtl number=1

42、 ，TDR prandtl number=1.3 。(3) 工作环境海底800n，工作压力8MP海流静止，温度2C。(4) 边界条件设置设定进出口流体旋转轴为y轴，进口流体的相对压力边界条件为 0MP速度初始条 件为海流静止；出口流体的初始条件与进口流体一致； 叶片壁面(包括ypqian和yphou) 为静止壁面(stational ),计算域流体设为动网格(movingmesh),旋转轴为y轴。(5) 残差收敛条件(con verge nee criteri on)continuity : 0.001 ; x-velocity : 0.001 ; y-velocity : 0.001 ;z-

43、velocity : 0.001 ; k: 0.001 ； epsilon : 0.001。(6) 设置3个监视器mon itor-1 :监测前叶片静压变化情况，报告类型为area-weighted average ，监测表面为ypqian ；monitor-2 :监测后叶片静压变化情况，报告类型为area-weighted average ，监测表面为yphou；monitor-3 :监测进出口质量流量平衡，报告类型(report type )为massflow rate， 检测表面为 pressure_inlet和 pressure_outlet 。3. 仿真运算(1)叶片转速为300rp

44、m,迭代次数iteraion=100 。通过迭代计算可以得到：(a)残差历史contin uity=1.9033e-04；x-velocity=9.2473e-04；y-velocity=9.8468e-04；z-velocity=8.6589e-04；k=5.0959e-04 ；e=9.8342e-04 ；满足残差收敛条件。残差历史曲线如图422所示：图422 残差历史(b)质量流量平衡曲线Mass Flow Rate刚开始时急剧上升最大时到达 20kg/s左右,主要由于刚开始旋转 时海流波动较大，经过100次迭代后，质量流量逐步降为-2.0131e-03kg/s，并在其附 近波动，基本保持

45、平衡，小于设定的收敛值，即达到了预期效果。通过后处理过程可以 得到：进口流量 pressurenlet: 7.3810244kg/s ；出 口流量 pressure_outlet : 7.379024kg/s ； 净流量 net mass-flow : -0.0020003319kg/s。质量流量曲线如图:图4.2.3质量流量平衡(C)阻力系数变化曲线如图424，由曲线可以看出，在转速为 300rpm时，开始时由于叶片处于静止的 海流中，所以阻力系数Cd波动较大，后来逐步趋向于稳定，稳定值为Cd =0.15，阻力系数的变化直接反映了阻力的变化，稳定时阻力并不是很大。图424阻力变化曲线（d）监

46、测叶片前后面静压变化通过监测叶片前后面可以得到此时前后面的静压情况，进而可以分析产生的推力ypqian 静压稳定值为-1.5872e+03pa，yphou 静压稳定值为 9.2495e+02pa。Area（e）利用FLUENT-3D进行后处理在后处理阶段，可以得到静压等高线图（如图和速度矢量图（如图 428，max=4.201029m/s,4.2.7 ，max=6150pa min=-2380pa)min=0.00999315m/s)。压差：p=2512.15pa;推力：F ps 2512.15 3.14 0.122113.589(N)。静压变化曲线如图425和图426 :427 叶片静压等高

47、线分布图428速度矢量分布图通过图像可以看出叶片的前后面有明显的压差，这也是推进器推力的来源，也可以看出在叶片的边缘速度比较大，叶片的连接处压力较大在高速旋转时会降低其使用寿 命，也容易造成损坏，从而会对整个 AUV的定位等产生影响，加大定位误差，也会在一 定程度上造成测量结果的精度降低。(2)下面分别将转速设为 400rpm, 500rpm, 600rpm, 800rpm和lOOOrpm,对比其 阻力系数及叶片前后面静压、速度、质量流量的变化(静压和速度矢量变化趋势与转速 为300rpm时类似，这里不再给出图像，只给出关键值，如表421所示。图429给出了转速和推力之间的关系曲线，可以直观的

48、看出二者的关系。450400350300260200卜50*or0Dnu3叶片推力与转速姜系曲线400500600 、700800900 WOO转 ij/rpm图429叶片推力与转速关系曲线表阻力系数、流体速度、静压对比叶片转速（rpm）3004005006008001000阻力系数cd50.40.751叶片 流体 速度（m/s）min9.993e-32.419e-24.484e-23.422e-23.115e-21.184e-1max4.2015.6867.1748.65911.74815.138质量流量（kg/s）in let-7.381-12.203-17.178-2

49、7.148-39.227-62.727outlet7.37912.20117.17427.19339.22562.736net-2e-3-2.348e-3-4.387e-3-1.846e-3-2.155e-3-9.789e-3叶片 静压（pa）ypqia n-1.587e+3-2.603e+3-3.952e+3-4.627e+3-7.320e+3-9.695e+3yephou9.250e+21.474e+32.164e+32.107e+32.705e+31.897e+3p2512.154077611667341002511592推力（Nps113.589184.346276.541304.48

50、5453.290524.144由表421可以从质量流量的结果看出，随着转速的增加，质量流量损失不断增加, 但在误差允许范围以内，总体而言，所建计算域是符合质量守恒定律的；根据阻力和推 力值，可以明显地看出，随着转速的不断增大，阻力系数不断增大，这意味着阻力不断 增大；同时，叶片旋转时，流体各处的速度大小不一，造成海水压强的变化，叶片前后 面间产生压差，利用压差推进器得以向前运动。 随着转速的增加，叶片前后面之间的压 差不断增加，这就导致了推力的增加。 当AUV在海底运动时，因推力的作用向前加速运 动，同时由于海水具有一定的粘性而产生阻力，速度越大阻力也越大，最后达到平衡。5总结从运动控制的角度出发，对AUV进行了力学研究和水动力分析，综合考虑其所受的 重力、浮力等，在三维模型的基础上，运用 ADAMS/MATLAB_Simulin进行运动学仿真, 对单个推进器运动控制进行研究与分析，通过仿真曲线较好地说明了所建立的控制系统 的正确性。通过FLUENT寸推进器进行流体仿真，得到