[发明专利]一种作业型ROV六自由度运动控制方法在审

专利信息
申请号: 201610003418.7 申请日: 2016-01-04
公开(公告)号: CN105607476A 公开(公告)日: 2016-05-25
发明(设计)人: 李新飞;袁利毫;韩端锋;崔洪斌;谢耀国;姜迈 申请(专利权)人: 哈尔滨工程大学
主分类号: G05B13/04 分类号: G05B13/04
代理公司: 暂无信息 代理人: 暂无信息
地址: 150001 黑龙江省哈尔滨市南岗区*** 国省代码: 黑龙江;23
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明属于无人水下机器人运动控制技术领域,具体涉及一种过驱动作业型遥控水下机器人的作业型ROV六自由度运动控制方法。本发明包括:建立ROV水动力学模型;建立作业型ROV在六个自由度方向的控制力和力矩向量;计算本体坐标系{b}中作业型ROV的运动状态;建立ROV的运动学方程;将ROV本体坐标系中的角速度运动状态转化为大地坐标系中来表示;将ROV本体坐标系中的角位移状态转化为大地坐标系中来表示;根据输入的大地坐标系中的ROV期望航向角,计算在本体坐标系中ROV的期望航向角。本发明中可以逼真的模拟通过外部控制输入来操纵作业型ROV运动及姿态变化的过程,具有能够逼真的模拟水下作业过程中ROV的各种运动及控制参数变化过程的优点。
搜索关键词: 一种 作业 rov 自由度 运动 控制 方法
【主权项】:
一种作业型ROV六自由度运动控制方法,其特征在于:(1)建立ROV水动力学模型:建立作业型ROV在本体坐标系{b}中的六自由度的水动力学方程:<mrow><mi>M</mi><mover><mi>&upsi;</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mi>C</mi><mrow><mo>(</mo><mi>&upsi;</mi><mo>)</mo></mrow><mi>&upsi;</mi><mo>+</mo><mi>D</mi><mrow><mo>(</mo><mi>&upsi;</mi><mo>)</mo></mrow><mi>&upsi;</mi><mo>+</mo><mi>g</mi><mrow><mo>(</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>=</mo><mi>&tau;</mi></mrow>式中:M∈(6×6)为系统惯性矩阵;C(υ)∈(6×6)为科氏向心力矩阵;D(υ)∈(6×6)为由水动力引起阻尼力矩阵;g(η)∈(6×1)为重力和浮力引起的力和力矩向量;υ∈(6×1)为ROV的速度和角速度向量;为ROV的加速度和角加速度向量;τ∈(6×1)为ROV的推力及推力矩向量;(2)建立作业型ROV在六个自由度方向的控制力和力矩向量:作业型ROV总共安装8台呈矢量布置的液压螺旋桨推进器,将在ROV本体坐标系{b}中的六个方向产生控制力及力矩向量:τ=[XT YT ZT KT MT NT]T(3)计算本体坐标系{b}中作业型ROV的运动状态:根据作业型ROV模型试验得到的水动力参数结果,计算出式中ROV水动力学模型中的系数矩阵M、C(υ)、D(υ)、g(η),将六自由度的控制力和力矩向量τ输入到ROV水动力学模型,即可得到本体坐标系{b}中的运动状态,包括速度[u v w p q r]T和位移[X Y Z α β γ]T;(4)建立ROV的运动学方程,将ROV本体坐标系{b}中的线速度状态转化为大地坐标系{n}中来表示:ROV在本体坐标系{b}中的速度向量为:υ=[v ω]T式中:v=[u v w]T为ROV在本体坐标系中三个线速度向量;ω=[p q r]T为ROV在本体坐标系中的三个角速度向量;大地坐标系中ROV的线速度向量为:<mrow><mover><mi>p</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>R</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mi>v</mi></mrow>式中:<mrow><mover><mi>p</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>N</mi><mo>&CenterDot;</mo></mover></mtd><mtd><mover><mi>E</mi><mo>&CenterDot;</mo></mover></mtd><mtd><mover><mi>D</mi><mo>&CenterDot;</mo></mover></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>为大地坐标系的线速度向量;Θ=[φ θ ψ]T为大地坐标系{n}和本体坐标系{b}之间的欧拉角向量;R(Θ)为本体坐标系{b}到大地坐标系{n}的线速度转换矩阵;<mrow><mi>R</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>c</mi><mi>&psi;</mi><mi>c</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>s</mi><mi>&psi;</mi><mi>c</mi><mi>&phi;</mi><mo>+</mo><mi>c</mi><mi>&psi;</mi><mi>s</mi><mi>&theta;</mi><mi>s</mi><mi>&phi;</mi></mrow></mtd><mtd><mrow><mi>s</mi><mi>&psi;</mi><mi>s</mi><mi>&phi;</mi><mo>+</mo><mi>c</mi><mi>&psi;</mi><mi>c</mi><mi>&phi;</mi><mi>s</mi><mi>&theta;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>s</mi><mi>&psi;</mi><mi>c</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mi>c</mi><mi>&psi;</mi><mi>c</mi><mi>&phi;</mi><mo>+</mo><mi>s</mi><mi>&phi;</mi><mi>s</mi><mi>&theta;</mi><mi>s</mi><mi>&psi;</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>c</mi><mi>&psi;</mi><mi>s</mi><mi>&phi;</mi><mo>+</mo><mi>s</mi><mi>&theta;</mi><mi>s</mi><mi>&psi;</mi><mi>c</mi><mi>&phi;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>s</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mi>c</mi><mi>&theta;</mi><mi>s</mi><mi>&phi;</mi></mrow></mtd><mtd><mrow><mi>c</mi><mi>&theta;</mi><mi>c</mi><mi>&phi;</mi></mrow></mtd></mtr></mtable></mfenced></mrow>式中:c表示为cos(·),s表示为sin(·);(5)建立ROV的运动学方程,将ROV本体坐标系{b}中的位移状态转化为大地坐标系{n}中来表示,大地坐标系{n}中ROV的位移向量:p=∫R(Θ)vdv式中:p=[N E D]T为大地坐标系{n}中的位移向量(6)将ROV本体坐标系中的角速度运动状态转化为大地坐标系中来表示:在大地坐标系中ROV的角速度向量可由下式计算:<mrow><mover><mi>&Theta;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>T</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mi>&omega;</mi></mrow>式中:<mrow><mover><mi>&Theta;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>&phi;</mi><mo>&CenterDot;</mo></mover></mtd><mtd><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover></mtd><mtd><mover><mi>&psi;</mi><mo>&CenterDot;</mo></mover></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>为大地坐标系中水下机器人的角速度向量;T(Θ)为本体坐标系到大地坐标系的角速度转换矩阵;<mrow><mi>T</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mi>s</mi><mi>&phi;</mi><mi>t</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mi>s</mi><mi>&phi;</mi><mi>t</mi><mi>&theta;</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>c</mi><mi>&phi;</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>s</mi><mi>&phi;</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>s</mi><mi>&phi;</mi><mo>/</mo><mi>c</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mi>c</mi><mi>&phi;</mi><mo>/</mo><mi>c</mi><mi>&theta;</mi></mrow></mtd></mtr></mtable></mfenced></mrow>式中:t表示为tan(·);(7)将ROV本体坐标系中的角位移状态转化为大地坐标系中来表示:大地坐标系中ROV的角位移向量:Θ=∫T(Θ)ωdω(8)计算本体坐标系ROV期望位移值和大地坐标系中ROV期望位移值之间的关系:本体坐标系中ROV三个方向的线速度:<mrow><mi>v</mi><mo>=</mo><msup><mi>R</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mover><mi>p</mi><mo>&CenterDot;</mo></mover></mrow>式中:R‑1(Θ)——R(Θ)的逆矩阵;本体坐标系中ROV三个方向的位移为:<mrow><mi>P</mi><mo>=</mo><munderover><mo>&Integral;</mo><mn>0</mn><mi>t</mi></munderover><mrow><mi>v</mi><mi>d</mi><mi>t</mi></mrow><mo>=</mo><munderover><mo>&Integral;</mo><mn>0</mn><mi>t</mi></munderover><msup><mi>R</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mover><mi>p</mi><mo>&CenterDot;</mo></mover><mi>d</mi><mi>t</mi></mrow>式中:P=[X Y Z]T为本体坐标系中的ROV的三个方向位移向量;(9)根据控制输入的大地坐标系中的ROV期望位移值[N0 E0 D0]T,计算在本体坐标系中ROV的期望位移值[X0 Y0 Z0]T[X0 Y0 Z0]T=R‑1(Θ)[N0 E0 D0]T(10)计算本体坐标系ROV期望姿态角和大地坐标系中ROV期望姿态角之间的关系:本体坐标系中ROV三个方向的角速度为:<mrow><mi>&omega;</mi><mo>=</mo><msup><mi>T</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mover><mi>&Theta;</mi><mo>&CenterDot;</mo></mover></mrow>式中:T‑1(Θ)为T(Θ)的逆矩阵;本体坐标系中ROV三个方向的角位移为:<mrow><mi>Q</mi><mo>=</mo><munderover><mo>&Integral;</mo><mn>0</mn><mi>t</mi></munderover><mi>&omega;</mi><mi>d</mi><mi>t</mi><mo>=</mo><munderover><mo>&Integral;</mo><mn>0</mn><mi>t</mi></munderover><msup><mi>T</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mover><mi>&Theta;</mi><mo>&CenterDot;</mo></mover><mi>d</mi><mi>t</mi></mrow>式中:Q=[α β γ]T为本体坐标系中的ROV的角位移向量;(11)根据输入的大地坐标系中的ROV期望航向角,计算在本体坐标系中ROV的期望航向角:在大地坐标系中ROV的期望姿态角为[φ0 θ0 ψ0]T,则本体坐标系中的期望姿态角值为[α0 β0 γ0]T为:[α0 β0 γ0]T=T‑1(Θ)[φ0 θ0 ψ0]T式中:α0为本体坐标系中的期望横倾角,对作业型ROV来说α0=0;β0为本体坐标系中的期望纵倾角,对作业型ROV来说β0=0;γ0为本体坐标系中的期望航向角;对作业型ROV来说,本体坐标系中的期望航向角为:γ0=ψ0;(12)建立作业型ROV外部控制输入的功能框图:ROV外部控制输入包括:ROV操纵面板输入控制和ROV软件输入控制,ROV操纵面板输入控制包括:操纵手柄输入、微调旋钮、自动按钮控制、ROV软件输入控制包括:键盘输入相对位置和软件自动控制;ROV外部控制命令的优先级是:ROV操纵面板输入控制和ROV软件输入控制优先级相同;对于ROV操纵面板输入控制来说:操纵手柄输入命令优先级高于自动按钮控制,保持高程和保持深度命令优先级相同,操纵手柄优先级高于微调旋钮命令;微调旋钮优先级高于自动按钮控制命令;对于ROV软件输入控制命令来说,其命令优先级的顺序与ROV操纵面板输入控制命令优先级顺序一致;(13)构建ROV六自由度运动控制器;(14)通过ROV外部控制输入来实现ROV的六自由度运动控制:ROV的各种外部控制命令经过控制命令逻辑判断模块,输出有效控制命令至ROV运动引导律计算模块和ROV运动控制器;ROV运动控制器输出六个方向的推力和推力矩,操纵ROV按照期望控制命令运动、定向或悬停定位。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工程大学,未经哈尔滨工程大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201610003418.7/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top