发明名称 一种DVL测速失效时的UUV反步滑模动力定位控制方法
摘要 一种DVL测速失效时的UUV反步滑模动力定位控制方法,涉及一种UUV动力定位控制方法。为了解决在UUV反步滑模动力定位控制中测速传感器DVL测量失效的问题。包括:在DVL测速失效时,测量UUV北向、东向位置和艏向角,速度估计器在线估计出北向、东向速度和艏向角。根据北向、东向、艏向角测量值与期望值的误差和估计的北向、东向速度、艏向角,构造使UUV渐进稳定的反步滑模控制律,解算出UUV在纵荡、横荡和艏摇三个自由度上的控制向量。控制向量经过推力分配,得到UUV主推进器、水平辅助推进器的推力,从而使UUV达到期望的北向、东向位置和艏向角,实现在DVL测速失效时UUV在水平面内的反步滑模动力定位控制。
申请公布号 CN105867417A 申请公布日期 2016.08.17
申请号 CN201610252029.8 申请日期 2016.04.21
申请人 哈尔滨工程大学 发明人 徐健;康晓峰;严浙平;张耕实;汪慢
分类号 G05D1/10(2006.01)I;G05B13/04(2006.01)I 主分类号 G05D1/10(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 岳昕
主权项 一种DVL测速失效时的UUV反步滑模动力定位控制方法,其特征在于,所述方法包括如下步骤:步骤1:UUV在水面进行纵荡、横荡和艏摇的三自由度反步滑模动力定位控制,获得期望向量η<sub>d</sub>=[x<sub>d</sub> y<sub>d</sub> ψ<sub>d</sub>]<sup>T</sup>,x<sub>d</sub>、y<sub>d</sub>和ψ<sub>d</sub>分别表示UUV在大地坐标下的北向位置、东向位置和艏向角的期望值;步骤2:对UUV的经纬度和艏向角进行采集,同时采用速度传感器DVL对UUV进行速度采集,确定向量η=[x y ψ]<sup>T</sup>,x、y和ψ分别表示采集的UUV在大地坐标下的北向位置、东向位置和艏向角的值;步骤3:故障检测模块实时检测,当检测到速度传感器DVL测速无效时,产生触发信息启动速度估计器对UUV的速度信息进行估计,速度估计器是以UUV的水平面三自由度动力学模型为基础,以采集的北向位置x、东向位置y和艏向角ψ作为速度估计器的输入,并且引入了北向位置、东向位置和艏向角估计误差<img file="FDA0000971277980000011.GIF" wi="187" he="54" />的滑模面函数以消除速度估计器模型的不确定性:<maths num="0001"><math><![CDATA[<mrow><mover><mover><mi>&eta;</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mo>=</mo><mi>J</mi><mover><mi>v</mi><mo>^</mo></mover><mo>+</mo><msub><mi>L</mi><mn>1</mn></msub><mi>s</mi><mi>a</mi><mi>t</mi><mrow><mo>(</mo><msub><mi>s</mi><mn>1</mn></msub><mo>/</mo><msub><mi>&phi;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971277980000012.GIF" wi="1198" he="71" /></maths><maths num="0002"><math><![CDATA[<mrow><mover><mover><mi>v</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mo>=</mo><msup><mi>M</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mo>-</mo><mi>C</mi><mover><mi>v</mi><mo>^</mo></mover><mo>-</mo><mi>D</mi><mover><mi>v</mi><mo>^</mo></mover><mo>-</mo><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mi>d</mi></msub><mo>+</mo><mi>&tau;</mi><mo>+</mo><msub><mi>L</mi><mn>2</mn></msub><mi>s</mi><mi>a</mi><mi>t</mi><mo>(</mo><mrow><msub><mi>s</mi><mn>1</mn></msub><mo>/</mo><msub><mi>&phi;</mi><mn>1</mn></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971277980000013.GIF" wi="1202" he="78" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mover><mover><mi>&tau;</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mi>d</mi></msub><mo>=</mo><mo>-</mo><msup><mi>T</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mi>d</mi></msub><mo>+</mo><msub><mi>L</mi><mn>3</mn></msub><mi>s</mi><mi>a</mi><mi>t</mi><mrow><mo>(</mo><msub><mi>s</mi><mn>1</mn></msub><mo>/</mo><msub><mi>&phi;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971277980000014.GIF" wi="1213" he="77" /></maths>其中,<img file="FDA0000971277980000015.GIF" wi="317" he="71" />为UUV在大地坐标系下的北向位置、东向位置和艏向角的估计向量,J∈R<sup>3×3</sup>是大地坐标系与UUV坐标系之间的变换矩阵,<img file="FDA0000971277980000016.GIF" wi="301" he="78" />为UUV坐标下的北向速度、东向速度和艏向角速度估计向量,sat(*)为饱和函数,<img file="FDA0000971277980000017.GIF" wi="190" he="61" />为速度估计器的滑模面函数,φ<sub>1</sub>为边界层厚度;M∈R<sup>3×3</sup>为UUV的惯性质量矩阵;C∈R<sup>3×3</sup>为UUV离心力和哥式力矩阵;D∈R<sup>3×3</sup>为粘性矩阵,<img file="FDA0000971277980000018.GIF" wi="173" he="70" />为外界环境的扰动估计向量,τ∈R<sup>3×1</sup>为控制力向量;L<sub>1</sub>、L<sub>2</sub>、L<sub>3</sub>∈R<sup>3×3</sup>为速度估计器增益对角矩阵;饱和函数sat(*)具体形式如下:<maths num="0004"><math><![CDATA[<mrow><mi>s</mi><mi>a</mi><mi>t</mi><mrow><mo>(</mo><msub><mi>s</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>sgn</mi><mrow><mo>(</mo><msub><mi>s</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>/</mo><msub><mi>&phi;</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mo>|</mo><msub><mi>s</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>/</mo><msub><mi>&phi;</mi><mn>1</mn></msub><mo>|</mo><mo>&gt;</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>s</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>/</mo><msub><mi>&phi;</mi><mn>1</mn></msub></mrow></mtd><mtd><mrow><mo>|</mo><msub><mi>s</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>/</mo><msub><mi>&phi;</mi><mn>1</mn></msub><mo>|</mo><mo>&le;</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971277980000019.GIF" wi="1374" he="151" /></maths>其中,φ<sub>1</sub>为正常数,表示边界层的厚度,i=1,2,3;步骤4:根据期望向量η<sub>d</sub>与向量<img file="FDA0000971277980000021.GIF" wi="59" he="54" />获得偏差向量<img file="FDA0000971277980000022.GIF" wi="250" he="62" />反步滑膜动力定位控制器再根据向量<img file="FDA0000971277980000023.GIF" wi="35" he="47" />和偏差向量η<sub>e</sub>,获得控制力向量τ;步骤5:根据获得的控制力向量τ进行推力分配,使UUV的主推进器和水平辅助推进器在分配得到的推力下工作,从而使UUV北向位置、东向位置和艏向角达到期望值。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号