CN103895814A - 一种船舶舵减横摇的控制方法 - Google Patents
一种船舶舵减横摇的控制方法 Download PDFInfo
- Publication number
- CN103895814A CN103895814A CN201410122987.4A CN201410122987A CN103895814A CN 103895814 A CN103895814 A CN 103895814A CN 201410122987 A CN201410122987 A CN 201410122987A CN 103895814 A CN103895814 A CN 103895814A
- Authority
- CN
- China
- Prior art keywords
- overbar
- delta
- plusminus
- rudder
- centerdot
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明涉及一种船舶舵减横摇的控制方法,其特征在于:对船舶舵减横摇运动进行数学描述;船舶横摇运动模型所输出的横摇角信号输入至高增益观测器,高增益观测器输出横摇角速度信号反馈输入至解析模型预测控制器;解析模型预测控制器输出指令舵角信号,实现船舶舵减横摇。
Description
技术领域
本发明涉及船舶工程、控制科学与控制工程领域,尤其涉及一种船舶舵减横摇的控制方法。
背景技术
由于船舶横摇运动阻尼很小,使得船舶在风浪中会产生剧烈的横摇,过大的横摇会对船舶航行性能和安全性带来很大的影响,为保证船舶在复杂海况安全航行,舵减横摇技术作为一种新型控制思想,近年来受到了极大的关注,而且用舵来进行减摇装置的设计,简单方便,价格低廉,因此,舵减摇技术是船舶运动控制领域中一个重要研究课题。船舶运动本质上是非线性的,从船舶运动环境和自身运动特征可知,其精确的数学模型难以得到,造成模型参数具有不确定性,因而,船舶舵减横摇控制器应该基于非线性控制理论来设计,且必须具有对模型参数摄动的鲁棒性。船舶舵减横摇控制要求在非线性模型具有参数不确定前提下对海浪引起的横摇运动抑制,舵减横摇摇控制器设计必须满足对模型参数摄动的鲁棒性。但传统PID,LQG控制算法,不能有效处理模型非线性和参数不确定性问题。
发明内容
本发明目的在于提供一种船舶舵减横摇的控制方法,能够有效解决非线性系统模型参数不确定性问题易于建模、响应迅速、控制性能好、鲁棒性强。
实现本发明目的技术方案:
一种船舶舵减横摇的控制方法,其特征在于:对船舶舵减横摇运动进行数学描述;船舶横摇运动模型所输出的横摇角信号输入至高增益观测器,高增益观测器输出横摇角速度信号反馈输入至解析模型预测控制器;解析模型预测控制器输出指令舵角信号,实现船舶舵减横摇。
船舶横摇运动模型所输出的横摇角信号反馈输入至解析模型预测控制器。
实时记录船舶的横摇角和控制舵角,进行控制效果监控,判断船舶是否能够抑制船舶横摇运动到期望程度,如果横摇抑制效果不好,则调整解析模型预测控制器。
对船舶舵减横摇运动进行数学描述具体包括以下内容,
船舶横向运动的非线性数学模型描述为
对照式(2),则得出船舶舵减横摇的不确定非线性系统为:
其中,
高增益观测器通过如下方法实现,
取高增益观测器状态方程如下:
0=-θp(θ)-(ATp(θ)+p(θ)A)+CTC (5)
其中,
p-1(θ)C0 T=[Cn 1θ,Cn 2θ2,…,Cn nθn]T (6)
将式(3)写成形如式(4)的状态方程为:
其中,
选取高增益观测器增益p-1(θ)C0 T=[2θ,θ2]T,式(7)表示为:
解析模型预测控制器通过如下方法实现,
1)含有不确定项的解析模型预测控制规律;
将带有不确定项的式(3)转换为:
其中,f*(x)=f(x)+Δf(x),g*(x)=g(x)+Δg(x),不失一般性,假设系统(11)的平衡点xo,有f*(xo)=0,g*(xo)≠0,h(xo)=0,称式(11)为重定义下的标称模型;
式(11)的滚动时域的性能函数为:
系统(11)在t时刻的预测控制问题描述为:
状态变量的初始值给为:
基于解析模型预测控制思想得到:
由式(16),可以得到带有不确定项的解析模型预测控制的最优控制律:
2)消除控制规律中的不确定项;
由式(11)可知:
则知g*(x)有界,又有:
其中,式(25)处的±符号选取与f(x)对应项相同;
边界值:
其中,式(26)b1项前的±符号与g(x)相同,含fd项前的符号与f(x)对应项相同;
取u1为其边界值,得到解析模型预测控制器:
其中,±符号选取与其对应前一项符号一致。
结合高增益观测器反馈信号,解析模型预测控制器通过如下方法实现控制,
其中,ur为解析模型预测控制器输出经过舵机执行器后的实际舵角。
本发明具有的有益效果:
针对非线性系统模型参数不确定性问题,本发明提出了一种改进解析模型预测控制方法。先将其重定义成带有不确定项标称系统模型,基于解析模型预测控制理论对重定义标称模型进行推导,得到含有不确定项的控制律。由于控制规律有界,取其边界值,从理论推导上消除了控制器中的不确定项,有效的解决了模型中的参数不确定性问题。船舶舵减横摇控制系统仅能测量船舶横摇角,本发明控制方法中引入状态观测器以获得控制器需要的横摇角速率信息。相比其他观测器,高增益观测器的稳定性和精度有理论保证,计算简便参数少,参数一旦选定,无需调整。由于高增益观测器设计过程中部分函数需满足局部李普希兹条件,这样舵减横摇控制最终需给出其满足设计高增益观测器所需条件,本发明高增益观测器设计通过构建满足局部李普希兹条件条件的观测器系数函数,实现对横摇角速度信息的观测。本发明通过状态观测器直接获得观测值,解决了横摇角速度不可测问题。本发明能够有效抑制模型的参数摄动,横摇运动可得到有效抑制,本发明具有易于建模、响应迅速、控制性能较好、鲁棒性强和逻辑结构简单等优点。
附图说明
图1为本发明控制方法流程图;
图2为基于本发明控制方法的控制系统组成框图。
具体实施方式
如图2所示,船舶横摇运动模型所输出的横摇角信号输入至高增益观测器,高增益观测器的输入信号还包括舵机实际输出的舵角信号。高增益观测器输出横摇角速度信号反馈输入至解析模型预测控制器;解析模型预测控制器输出指令舵角信号,实现船舶舵减横摇。船舶横摇运动模型所输出的横摇角信号反馈输入至解析模型预测控制器,解析模型预测控制器的输入信号还包括期望横摇角指令。
本发明是由船舶舵减横摇系统的改进解析模型预测控制器接受期望横摇角信息、高增益观测器信息及罗经获得的系统输出横摇角信息,经过运算输出指令舵角,舵机执行器根据舵角指令输出实际舵角给船体,船舶抑制产生的横摇角,输出的横摇角经由罗经测得并输入给高增益观测器和控制器。实时记录船舶的横摇角和控制舵角,进行控制效果监控,判断船舶是否能够抑制船舶横摇运动到期望程度,如果横摇抑制效果不好,则调整解析模型预测控制器,如此形成闭环控制系统。
一、对船舶舵减横摇运动进行数学描述,即完成船舶舵减横摇运动的数学模型描述,海浪扰动对船舶横摇运动产生的扰动力矩的数学描述,由于船舶航行速度及初稳性高度变化对船舶舵减横摇系统水动力系数影响的描述三项工作。具体包括以下内容:
船舶横向运动的非线性数学模型可以描述为
其中, g(x)=b,h(x)=x1,w为海浪扰动,a1,a2,a3,a4为计算出的已知系数。而由于速度和初稳性高度变化引起系数的变化量为:Δa1,Δa2,Δa3,Δa4,Δb,这些量的变化一般不能准确求出,但可以预估计其最大绝对值,记为:mi(i=1,2,3,4)及b1。
对照式(2),则得出船舶舵减横摇的不确定非线性系统为:
其中,
二、高增益观测器的设计和舵减横摇系统满足设计高增益观测器所需条件的证明
1)高增益观测器的设计。
取高增益观测器状态方程如下:
0=-θp(θ)-(ATp(θ)+p(θ)A)+CTC (10)
其中,
p-1(θ)C0 T=[Cn 1θ,Cn 2θ2,…,Cn nθn]T (6)
将式(3)写成形如式(4)的状态方程为:
其中,
选取高增益观测器增益p-1(θ)C0 T=[2θ,θ2]T,式(7)表示为:
2)舵减横摇系统满足设计高增益观测器所需条件的证明。
所设计的高增益观测器必须满足下面两个假设条件:
假设1:b(·),φ(·)满足局部李普希兹条件。
假设2:控制输入为有界的状态估计。
令m=x1,n=x2。假设x1,x2∈[-1.57,1.57],要证明存在L1,L2>0,使得,|φ(m,n1)-φ(m,n2)|≤L1|n1-n2|和|φ(m1,n)-φ(m2,n)|≤L2|m1-m2|同时成立。
取L1=3.14|a2+Δa2|,|φ(m,n1)-φ(m,n2)|≤L1|n1-n2|成立;
取L2=7.40|a2+Δa2|,|φ(m1,n)-φ(m2,n)|≤L2|m1-m2|成立。
由于|u|≤35°,输入有界,满足假设2)。
三、非线性解析模型预测控制控制器设计。从理论推导上消除了控制器中的不确定项,完成船舶舵减横摇改进非线性解析模型预测控制器设计,具体为
1)含有不确定项的解析模型预测控制规律。
将带有不确定项的式(3)转换为:
其中,f*(x)=f(x)+Δf(x),g*(x)=g(x)+Δg(x)。不失一般性,假设系统(11)的平衡点xo,有f*(xo)=0,g*(xo)≠0,h(xo)=0,称式(11)为重定义下的标称模型。
系统(11)的滚动时域的性能函数为:
系统(11)在t时刻的预测控制问题描述为:
基于解析模型预测控制思想得到:
由式(16),可以得到带有不确定项的解析模型预测控制的最优控制律:
2)消除控制规律中的不确定项。
由式(11)可知:
则知f*(x)有界,并记
则知g*(x)有界。又有:
其中,式(25)处的±符号选取与f(x)对应项相同。
边界值:
其中,式(26)b1项前的±符号与g(x)相同,含fd项前的符号与f(x)对应项相同。
综上,可知u(t)有界,取u1为其边界值,得到改进的舵减横摇解析模型预测控制器:
其中,±符号选取与其对应前一项符号一致。
结合高增益观测器和改进的解析模型预测控制律,船舶舵减横摇控制系统为:
其中,ur为控制器输出经过舵机执行器后的实际舵角。
五、船舶舵减横摇系统的控制效果监控。实时记录船舶的横摇角和控制舵角,判断船舶是否能够抑制船舶横摇运动到期望程度,如果能够满足程序结束,如果横摇抑制效果不好,调整控制器的设计,如图1所示。
Claims (7)
1.一种船舶舵减横摇的控制方法,其特征在于:对船舶舵减横摇运动进行数学描述;船舶横摇运动模型所输出的横摇角信号输入至高增益观测器,高增益观测器输出横摇角速度信号反馈输入至解析模型预测控制器;解析模型预测控制器输出指令舵角信号,实现船舶舵减横摇。
2.根据权利要求1所述的船舶舵减横摇的控制方法,其特征在于:船舶横摇运动模型所输出的横摇角信号反馈输入至解析模型预测控制器。
3.根据权利要求2所述的船舶舵减横摇的控制方法,其特征在于:实时记录船舶的横摇角和控制舵角,进行控制效果监控,判断船舶是否能够抑制船舶横摇运动到期望程度,如果横摇抑制效果不好,则调整解析模型预测控制器。
5.根据权利要求4所述的船舶舵减横摇的控制方法,其特征在于:高增益观测器通过如下方法实现,
取高增益观测器状态方程如下:
其中,“-”代表观测,为非线性项,p(θ)定义为如下方程的解:
0=-θp(θ)-(ATp(θ)+p(θ)A)+CTC (5)
其中,
p-1(θ)C0 T=[Cn 1θ,Cn 2θ2,…,Cn nθn]T (6)
将式(3)写成形如式(4)的状态方程为:
其中,
选取高增益观测器增益p-1(θ)C0 T=[2θ,θ2]T,式(7)表示为:
6.根据权利要求5所述的船舶舵减横摇的控制方法,其特征在于:解析模型预测控制器通过如下方法实现,
1)含有不确定项的解析模型预测控制规律;
将带有不确定项的式(3)转换为:
其中,f*(x)=f(x)+Δf(x),g*(x)=g(x)+Δg(x),不失一般性,假设系统(11)的平衡点xo,有f*(xo)=0,g*(xo)≠0,h(xo)=0,称式(11)为重定义下的标称模型;
式(11)的滚动时域的性能函数为:
系统(11)在t时刻的预测控制问题描述为:
基于解析模型预测控制思想得到:
其中,K=[k0,k1,…,kρ-1]代表矩阵Γll -1Γρl T的第一行元素,由预测周期T1、控制阶次l及相关度ρ决定,为最优控制律
由式(16),可以得到带有不确定项的解析模型预测控制的最优控制律:
2)消除控制规律中的不确定项;
由式(11)可知:
则知g*(x)有界,又有:
其中,式(25)处的±符号选取与f(x)对应项相同;
其中,式(26)b1项前的±符号与g(x)相同,含fd项前的符号与f(x)对应项相同;
取u1为其边界值,得到解析模型预测控制器:
其中,±符号选取与其对应前一项符号一致。
7.根据权利要求6所述的船舶舵减横摇的控制方法,其特征在于:结合高增益观测器反馈信号,解析模型预测控制器通过如下方法实现控制,
其中,ur为解析模型预测控制器输出经过舵机执行器后的实际舵角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410122987.4A CN103895814B (zh) | 2014-03-28 | 2014-03-28 | 一种船舶舵减横摇的控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410122987.4A CN103895814B (zh) | 2014-03-28 | 2014-03-28 | 一种船舶舵减横摇的控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103895814A true CN103895814A (zh) | 2014-07-02 |
CN103895814B CN103895814B (zh) | 2016-09-14 |
Family
ID=50987472
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410122987.4A Expired - Fee Related CN103895814B (zh) | 2014-03-28 | 2014-03-28 | 一种船舶舵减横摇的控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103895814B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104808662A (zh) * | 2015-03-13 | 2015-07-29 | 哈尔滨工程大学 | 一种基于数据驱动的抑制船舶航向扰动的控制方法 |
CN106842910A (zh) * | 2016-11-17 | 2017-06-13 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 一种基于干扰观测器的船舶减摇鳍滑模控制方法 |
CN107145072A (zh) * | 2017-05-31 | 2017-09-08 | 威海海洋职业学院 | 船舶舵减横摇自适应控制系统 |
CN107140110A (zh) * | 2017-03-21 | 2017-09-08 | 山东省科学院海洋仪器仪表研究所 | 一种船舶大幅横摇运动非线性阻尼系数识别方法 |
CN108508873A (zh) * | 2018-05-03 | 2018-09-07 | 厦门大学 | 超声波辅助检测的潜艇舵系统故障诊断与容错控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1429736A (zh) * | 2001-12-31 | 2003-07-16 | 松村纪孝 | 船舶减摇装置及其控制方法 |
CN1521083A (zh) * | 2003-02-12 | 2004-08-18 | 松村纪孝 | 船舶减摇装置的控制方法 |
CN102298348A (zh) * | 2011-07-05 | 2011-12-28 | 哈尔滨工程大学 | 船舶舵与翼舵或鳍与翼鳍联合控制系统数字控制器及方法 |
CN102910265A (zh) * | 2012-10-31 | 2013-02-06 | 哈尔滨工程大学 | 一种基于双重控制h∞设计的舵鳍联合减摇方法 |
CN102951272A (zh) * | 2012-10-30 | 2013-03-06 | 哈尔滨工程大学 | 一种船舶舵减摇的综合再置控制方法 |
EP2312271B1 (de) * | 2009-09-10 | 2013-10-23 | ThyssenKrupp Marine Systems GmbH | Verfahren zur Bestimmung von Seegangsgrößen |
KR20140037690A (ko) * | 2012-09-19 | 2014-03-27 | 삼성중공업 주식회사 | 선박의 운항안정성 제어장치 및 방법 |
-
2014
- 2014-03-28 CN CN201410122987.4A patent/CN103895814B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1429736A (zh) * | 2001-12-31 | 2003-07-16 | 松村纪孝 | 船舶减摇装置及其控制方法 |
CN1521083A (zh) * | 2003-02-12 | 2004-08-18 | 松村纪孝 | 船舶减摇装置的控制方法 |
EP2312271B1 (de) * | 2009-09-10 | 2013-10-23 | ThyssenKrupp Marine Systems GmbH | Verfahren zur Bestimmung von Seegangsgrößen |
CN102298348A (zh) * | 2011-07-05 | 2011-12-28 | 哈尔滨工程大学 | 船舶舵与翼舵或鳍与翼鳍联合控制系统数字控制器及方法 |
KR20140037690A (ko) * | 2012-09-19 | 2014-03-27 | 삼성중공업 주식회사 | 선박의 운항안정성 제어장치 및 방법 |
CN102951272A (zh) * | 2012-10-30 | 2013-03-06 | 哈尔滨工程大学 | 一种船舶舵减摇的综合再置控制方法 |
CN102910265A (zh) * | 2012-10-31 | 2013-02-06 | 哈尔滨工程大学 | 一种基于双重控制h∞设计的舵鳍联合减摇方法 |
Non-Patent Citations (3)
Title |
---|
杜佳璐等: "基于高增益观测器的船舶动力定位系统的输出反馈控制", 《控制理论与应用》, vol. 30, no. 11, 30 November 2013 (2013-11-30), pages 1486 - 1491 * |
王晓飞: "基于解析模型预测控制的欠驱动船舶路径跟踪控制研究", 《万方学位论文》, 3 August 2011 (2011-08-03), pages 28 - 57 * |
赖延辉等: "舵减横摇研究综述", 《海军工程学院学报》, no. 3, 30 September 1992 (1992-09-30), pages 6 - 15 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104808662A (zh) * | 2015-03-13 | 2015-07-29 | 哈尔滨工程大学 | 一种基于数据驱动的抑制船舶航向扰动的控制方法 |
CN106842910A (zh) * | 2016-11-17 | 2017-06-13 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 一种基于干扰观测器的船舶减摇鳍滑模控制方法 |
CN107140110A (zh) * | 2017-03-21 | 2017-09-08 | 山东省科学院海洋仪器仪表研究所 | 一种船舶大幅横摇运动非线性阻尼系数识别方法 |
CN107140110B (zh) * | 2017-03-21 | 2019-07-09 | 山东省科学院海洋仪器仪表研究所 | 一种船舶大幅横摇运动非线性阻尼系数识别方法 |
CN107145072A (zh) * | 2017-05-31 | 2017-09-08 | 威海海洋职业学院 | 船舶舵减横摇自适应控制系统 |
CN108508873A (zh) * | 2018-05-03 | 2018-09-07 | 厦门大学 | 超声波辅助检测的潜艇舵系统故障诊断与容错控制方法 |
CN108508873B (zh) * | 2018-05-03 | 2020-01-03 | 厦门大学 | 超声波辅助检测的潜艇舵系统故障诊断与容错控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103895814B (zh) | 2016-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104076821B (zh) | 基于模糊自适应观测器的欠驱动水面艇轨迹跟踪控制系统 | |
CN103895814A (zh) | 一种船舶舵减横摇的控制方法 | |
CN103970021B (zh) | 一种基于模型预测控制的松弛动力定位控制系统 | |
Banazadeh et al. | Frequency domain identification of the Nomoto model to facilitate Kalman filter estimation and PID heading control of a patrol vessel | |
CN103676654B (zh) | 基于干扰补偿的船舶动力定位速度估计系统及方法 | |
CN101707016A (zh) | 航海模拟器用船舶六自由度运动数学模型的建立方法 | |
CN103217160A (zh) | 一种船舶动力定位卡尔曼滤波方法 | |
CN104635773A (zh) | 一种基于改进强跟踪滤波状态观测器的船舶动力定位方法 | |
Kjerstad et al. | Disturbance rejection by acceleration feedforward for marine surface vessels | |
CN104316025A (zh) | 一种基于船体姿态信息估计海浪浪高的系统 | |
CN104142626B (zh) | 一种基于逆系统和内模控制的船舶动力定位控制方法 | |
CN104808662B (zh) | 一种基于数据驱动的抑制船舶航向扰动的控制方法 | |
CN104155043B (zh) | 一种动力定位系统外界环境力测量方法 | |
CN103935480A (zh) | 一种基于解析模型预测控制设计的船舶舵减摇方法 | |
Tomera | Nonlinear observers design for multivariable ship motion control | |
Battista et al. | An energy-casimir approach to underwater vehicle depth and heading regulation in short crested waves | |
Saoud et al. | Modeling and control design of a robotic sailboat | |
CN112965363A (zh) | 一种基于外源卡尔曼滤波的船舶动力定位系统 | |
Brodtkorb et al. | Sea state estimation using model-scale DP measurements | |
Auestad et al. | Heave motion estimation on a craft using a strapdown inertial measurement unit | |
Valentinis et al. | A Maneuvering Model for an Underwater Vehicle Near a Free Surface—Part III: Simulation and Control Under Waves | |
CN105180944A (zh) | 一种船体侧滑速度误差判定补偿方法 | |
CN111651943B (zh) | 一种船舶运动控制方法、终端设备及存储介质 | |
Moe et al. | Path following of underactuated marine underwater vehicles in the presence of unknown ocean currents | |
Wang et al. | Nonsingular fast terminal sliding mode-based robust adaptive structural reliable position-mooring control with uncertainty estimation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160914 |
|
CF01 | Termination of patent right due to non-payment of annual fee |