发明名称 一种动力定位船舶循迹导引控制方法
摘要 本发明公开了一种动力定位船舶循迹导引控制方法。该算法首先建立大地坐标系和船体坐标系,以及船舶的数学模型;然后通过设定路径点生成的路径与船舶之间的几何位置关系生成导引策略,计算船舶当前时刻的期望艏向;应用控制算法得到控制船舶达到期望艏向所需的控制力矩;最后计算控制船舶达到期望运动速度所需要的纵向推力。该算法根据航迹信息以及航迹和船舶的相对位置计算高速循迹直线段和转弯段船舶的期望艏向,并通过纵向推力控制船舶的运动速度,计算过程简单易行,适合于工程应用。
申请公布号 CN103576555A 申请公布日期 2014.02.12
申请号 CN201310563586.8 申请日期 2013.11.14
申请人 哈尔滨工程大学 发明人 付明玉;张爱华;余玲玲;谢笑颖;焦建芳
分类号 G05B13/04(2006.01)I;G05D1/02(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 代理人
主权项 一种动力定位船舶循迹导引控制方法,其特征在于包括如下步骤:(1)选取海平面上任意一点为原点,以正东方向为横轴、正北方向为纵轴,建立大地坐标系EOEN;然后,以船舶的几何中心为原点,以右舷方向为横轴、船艏方向为纵轴,建立船体坐标系YObX;(2)在大地坐标系和船体坐标系中建立动力定位船舶三自由度船舶高速运动数学模型: <mfenced open='' close=''> <mtable> <mtr> <mtd> <mover> <mi>&eta;</mi> <mo>&CenterDot;</mo> </mover> <mo>=</mo> <mi>Jv</mi> </mtd> </mtr> <mtr> <mtd> <mi>M</mi> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> <mo>+</mo> <mi>D</mi> <mrow> <mo>(</mo> <mi>v</mi> <mo>)</mo> </mrow> <mi>v</mi> <mo>+</mo> <mi>C</mi> <mrow> <mo>(</mo> <mi>v</mi> <mo>)</mo> </mrow> <mi>v</mi> <mo>=</mo> <mi>&tau;</mi> </mtd> </mtr> </mtable> </mfenced>式中,η表示船舶在大地坐标系中位置和艏向向量[n,e,ψ]Τ;v表示船舶在船体坐标系中线速度和角速度向量[u,v,r]Τ;τ为推进器产生的力和力矩向量[Fx,Fy,Nz]Τ;J为从船体坐标到大地坐标的转换矩阵, <mrow> <mi>J</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>cos</mi> <mi>&psi;</mi> </mtd> <mtd> <mo>-</mo> <mi>sin</mi> <mi>&psi;</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <mi>&psi;</mi> </mtd> <mtd> <mi>cos</mi> <mi>&psi;</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>M为船舶惯性矩阵, <mrow> <mi>M</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>m</mi> <mo>-</mo> <msub> <mi>X</mi> <mover> <mi>u</mi> <mo>.</mo> </mover> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>m</mi> <mo>-</mo> <msub> <mi>Y</mi> <mover> <mi>v</mi> <mo>.</mo> </mover> </msub> </mtd> <mtd> <msub> <mi>mx</mi> <mi>g</mi> </msub> <mo>-</mo> <msub> <mi>Y</mi> <mover> <mi>r</mi> <mo>.</mo> </mover> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>mx</mi> <mi>g</mi> </msub> <mo>-</mo> <msub> <mi>N</mi> <mover> <mi>v</mi> <mo>.</mo> </mover> </msub> </mtd> <mtd> <msub> <mi>I</mi> <mi>z</mi> </msub> <mo>-</mo> <msub> <mi>N</mi> <mover> <mi>r</mi> <mo>.</mo> </mover> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>m为船舶质量;Iz为转动惯量;xg为船体坐标系中船舶质心纵向坐标;其他参数均为一阶水动力导数;船舶水动力阻尼项为D(v)=D+Dn(v); <mrow> <mi>D</mi> <mo>=</mo> <mo>-</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>X</mi> <mi>u</mi> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>Y</mi> <mi>v</mi> </msub> </mtd> <mtd> <msub> <mi>Y</mi> <mi>r</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>N</mi> <mi>v</mi> </msub> </mtd> <mtd> <msub> <mi>N</mi> <mi>r</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>为线性部分; <mrow> <msub> <mi>D</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>X</mi> <mrow> <mi>u</mi> <mo>|</mo> <mi>u</mi> <mo>|</mo> </mrow> </msub> <mo>|</mo> <mi>u</mi> <mo>|</mo> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>Y</mi> <mrow> <mi>v</mi> <mo>|</mo> <mi>v</mi> <mo>|</mo> </mrow> </msub> <mo>|</mo> <mi>v</mi> <mo>|</mo> <mo>+</mo> <msub> <mi>Y</mi> <mrow> <mo>|</mo> <mi>r</mi> <mo>|</mo> <mi>v</mi> </mrow> </msub> <mo>|</mo> <mi>r</mi> <mo>|</mo> </mtd> <mtd> <msub> <mi>Y</mi> <mrow> <mo>|</mo> <mi>v</mi> <mo>|</mo> <mi>r</mi> </mrow> </msub> <mo>|</mo> <mi>v</mi> <mo>|</mo> <mo>+</mo> <msub> <mi>Y</mi> <mrow> <mo>|</mo> <mi>r</mi> <mo>|</mo> <mi>r</mi> </mrow> </msub> <mo>|</mo> <mi>r</mi> <mo>|</mo> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>N</mi> <mrow> <mi>v</mi> <mo>|</mo> <mi>v</mi> <mo>|</mo> </mrow> </msub> <mo>|</mo> <mi>v</mi> <mo>|</mo> <mo>+</mo> <msub> <mi>N</mi> <mrow> <mo>|</mo> <mi>r</mi> <mo>|</mo> <mi>v</mi> </mrow> </msub> <mo>|</mo> <mi>r</mi> <mo>|</mo> </mtd> <mtd> <msub> <mi>N</mi> <mrow> <mo>|</mo> <mi>v</mi> <mo>|</mo> <mi>r</mi> </mrow> </msub> <mo>|</mo> <mi>v</mi> <mo>|</mo> <mo>+</mo> <msub> <mi>N</mi> <mrow> <mo>|</mo> <mi>r</mi> <mo>|</mo> <mi>r</mi> </mrow> </msub> <mo>|</mo> <mi>r</mi> <mo>|</mo> </mtd> </mtr> </mtable> </mfenced> </mrow>为非线性部分; <mrow> <mi>C</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <msub> <mi>Y</mi> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> </msub> <mo>)</mo> </mrow> <mi>v</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>mx</mi> <mi>g</mi> </msub> <mo>-</mo> <msub> <mi>Y</mi> <mover> <mi>r</mi> <mo>&CenterDot;</mo> </mover> </msub> <mo>)</mo> </mrow> <mi>r</mi> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <msub> <mi>X</mi> <mover> <mi>u</mi> <mo>&CenterDot;</mo> </mover> </msub> <mo>)</mo> </mrow> <mi>u</mi> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <msub> <mi>Y</mi> <mover> <mi>v</mi> <mo>&CenterDot;</mo> </mover> </msub> <mo>)</mo> </mrow> <mi>v</mi> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>mx</mi> <mi>g</mi> </msub> <mo>-</mo> <msub> <mi>Y</mi> <mover> <mi>r</mi> <mo>&CenterDot;</mo> </mover> </msub> <mo>)</mo> </mrow> <mi>r</mi> </mtd> <mtd> <mo>-</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <msub> <mi>X</mi> <mover> <mi>u</mi> <mo>&CenterDot;</mo> </mover> </msub> <mo>)</mo> </mrow> <mi>u</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>为科里奥利和向心力矩阵;(3)通过船体设定的路径点和当前时刻船舶的大地坐标位置,获得当前时刻的船舶期望艏向;(4)利用反步法,根据当前时刻期望艏向ψd与船舶的期望速度ud,计算控 制船舶高速循迹所需的纵向力Fx和力矩Nz;(5)将控制力和力矩传递给动力定位船舶的数学模型,解算下一时刻船舶的状态信息,然后由当前时刻更新到下一时刻;(6)重复执行步骤(2)‑‑步骤(5),直到最后一条航迹循迹完成。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室