CN104133374B - 一种船舶及石油平台动力定位系统的推力分配方法 - Google Patents
一种船舶及石油平台动力定位系统的推力分配方法 Download PDFInfo
- Publication number
- CN104133374B CN104133374B CN201410348488.7A CN201410348488A CN104133374B CN 104133374 B CN104133374 B CN 104133374B CN 201410348488 A CN201410348488 A CN 201410348488A CN 104133374 B CN104133374 B CN 104133374B
- Authority
- CN
- China
- Prior art keywords
- alpha
- thrust
- propeller
- angle
- ship
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 230000002153 concerted effect Effects 0.000 claims description 4
- 238000000205 computational method Methods 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 239000000523 sample Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Landscapes
- Earth Drilling (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明是一种船舶及石油平台动力定位系统的推力分配方法。包括如下步骤:步骤一:推进器推进角的确定;步骤二:求取结构矩阵;步骤三:伪逆算法求推力分配;步骤四:推力限幅确定。本发明实现了船舶及石油平台动力定位推力分配的快速性,稳定性,高效性,是一种方便实用的船舶及石油平台动力定位系统的推力分配方法。
Description
技术领域
本发明是一种船舶及石油平台动力定位系统的推力分配方法,特别是一种水面船舶及石油平台动力定位系统的多螺旋桨推力分配方法,属于船舶及石油平台动力定位系统的推力分配方法的创新技术。
背景技术
针对船舶/石油平台平面定位的需要,采用多组螺旋桨推力以消除船舶/石油平台上受到的海风、海浪、海流的影响,使船舶/石油平台等按预定的位置、航向、轨迹、速度保持姿态。船舶/石油平台动力定位系统是一种闭环控制系统,它是通过推进器提供的推力,使海洋结构物在海上保持一定的位置和艏向角的一种先进的定位系统。它的基本原理是通过获得海洋结构物的位置误差数据,由控制器计算出抗拒环境干扰所需的推力和转矩,并通过推力系统的实现,以此达到消除位置误差的目的。此动力定位系统一般由控制系统、测量系统和推力系统组成。推力系统是船舶动力定位系统的一个重要组成部分,作为船舶动力定位系统的执行机构,其作用是提供控制系统所需要的推力和力矩,用来满足定位的需要,推力系统主要由发电机组、输电设备、推进器等组成。推进器一般来说是指螺旋桨,应用比较广泛的是槽道式推进器,全回转推进器和吊舱推进器。
伪逆算法是动力定位工程应用中普遍采用的方法,这种算法无需多次进行迭代计算,方法简单,但无法解决推进角及禁区角问题。现有的技术当中,有利用蚁群算法、遗传算法、二次序列规划算法等来解决推力分配最优解问题,但这些算法需要多次迭代,计算时间长,不易用于船舶/石油平台的实时操作控制中。
发明内容
本发明的目的在于考虑上述问题而提供一种实现了船舶及石油平台动力定位推力分配的快速性,稳定性,高效性的船舶及石油平台动力定位系统的推力分配方法。
本发明的技术方案是:本发明的船舶及石油平台动力定位系统的推力分配方法,包括如下步骤:步骤一:推进器推进角的确定;步骤二:求取结构矩阵;步骤三:伪逆算法求推力分配;步骤四:推力限幅确定。
本发明针对可调节的带方向角约束的多螺旋桨的推进角度计算问题和各螺旋桨的推力分配问题,即带约束的多参数的快速最优解的求解问题,本发明提出先分配到推进角方法,实现推力方向角的快速跟踪和推力的高效性;再分配到推力的方法,采用伪逆算法实现不同推进能力的螺旋桨的协调工作,有效的避免了相邻近的推进器之间推力干扰引起的推力衰减问题,解决了其他算法中的多奇异值问题,多次迭代计算问题,克服了收敛时间长的缺点,实现了船舶/石油平台动力定位推力分配计算的快速性,稳定性,高效性,易于在船舶、石油平台上实现。本发明是一种方便实用的船舶及石油平台动力定位系统的推力分配方法。
附图说明
图1为船舶/石油平台推进器分布与推进角度示意图;
图2为船舶/石油平台推进器间距与偏差角关系图;
图3为船舶/石油平台动力定位系统推力分配程序框图;
图4为推力分配结果对比仿真图;
图5为各螺旋桨推力分配结果仿真图。
具体实施方式
实施例:
本发明船舶及石油平台动力定位系统的推力分配方法,包括如下步骤:
步骤一:推进器推进角的确定;
步骤二:求取结构矩阵;
步骤三:伪逆算法求推力分配;
步骤四:推力限幅确定。
上述步骤一推进器推进角的确定方法如下:
1)计算全回转推进器的期望推进方向的方向角α
取船舶所受合力的方向角作为全回转推进器期望推进方向的方向角α,以下称推进器推进角,则全回转推进器推进角α为:
定义τ=[Fu,Fv,Nz]T,Fu,Fv,Nz分别为船舶在纵荡、横荡和艏摇三个方向上的所受合力(矩),α以船头正前方为零度,从俯视看,顺时针增加,一周为360°;Fu,Fv,Nz的正方向分别为船头正前方,船体正右侧,从俯视看,顺时针方向;
2)全回转推进器推进角α的优化
本发明中的推进器推进角度的计算方法和伪逆算法,能实现对连续变化的给定信号的快速跟踪,如果期望推进器推进角α和当前时刻α'之差满足下式:
取
代替α;
3)确定两个相邻的全回转推进器推进角
如果两个相邻的全回转推进器i和j不满足合理间距di,j>2(Di+Dj),引入推进器偏差角Δαi,j,用于解决两个相近的全回转螺旋桨之间产生的推力干扰问题,Di、Dj分别为第i、第j个推进器的直径,则推进器推进角满足:
其中αi、αj分别为第i、第j个推进器的推进角;推进器偏差角Δαi,j满足:
上述步骤二求取结构矩阵的方法是:
定义τ=A(α)f,f=[f1,…,fp,fp-1,…,fn]T,fi表示第i个推进器的推力;A为结构矩阵,由推进器的类型、位置分布以及推进器角度得到:
其中,Lix,Liy分别代表第i个推进器在纵荡方向和横荡方向的力臂。
上述步骤三伪逆算法求推力分配的方法是:
根据优化的目标函数:以及约束条件:S.t.τ=A(α)·f,求解推力f,其中W为权重矩阵,定义推力的拉格朗日函数为:
L(f,λ)=fTWf+λT(τ-A(α)f) (1.6)
推力的拉格朗日函数对变量f求导:
由局部最优定理可知:
f=(W+WT)-1A(α)Tλ (1.9)
(W+WT)T>0时,u取最小值,其中W+WT为奇异矩阵,且Wii>0(i=1…n),由边界条件可知:
τ=A(α)f=A(α)(W+WT)-1A(α)Tλ
(2.0) 从而得到:
λ=(A(α)(W+WT)-1A(α)T)-1τ (2.1)
代入上式,可得推力:
f=(W+WT)-1A(α)T(A(α)(W+WT)-1A(α)T)-1τ (2.2)如果A(α)-1存在,由式(2.2),得
f=(W+WT)-1A(α)T(A(α)T)-1(W+WT)A(α)-1τ=A(α)-1τ
(2.3)
即,当A(α)奇异矩阵时,(2.2)式可实现推力的逆分配,然而结构矩阵A(α)一般为非奇异矩阵,因此我们利用W+WT实现推力的伪逆分配;
当W≠I(I为单位矩阵),(2.2)式为带权重的伪逆推力分配算法;当W=I,为等价或者无权重的伪逆推力分配算法。
上述步骤四推力限幅的确定方法是:判断推进器推力f是否超过边界值,如果有推进器推力超出此推进器推力的边界值,则取此边界值作为此推进器的输出,实现推力分配。
图1为船舶/石油平台推进器分布示意图。此示意图是以一艘按26:1缩小的2.8m供给船,以此模型为计算对象,说明本发明方法。1#和2#推进器均位于船头,分布在船舶的纵向中轴线上;3#和4#推进器均位于船尾,对称分布在船舶的纵向中轴线的两侧。1#和2#推进器均采用槽道式推进器,安装有转速传感器;3#和4#推进器采用全回转推进器,安装有转速传感器和全回转角度传感器。实时控制计算机、输入输出电路、推进器、传感器、采样电路等共同组成整个推力分配系统。推进器1#~2#可分别产生最大4.5N推力,3#~4#可分别产生最大10N推力。全回转推进器3#和4#的推进角最大变化速度max|Δα|=30°/s。
各个推进器参数如下表所示:
推进器个数大于方向度3,因此可以产生多个不同的推力和方向组合,满足推力和力矩要求,本发明可以从中找到最优解。
图3是船舶动力定位系统推力分配程序框图,程序开始,其步骤如下:
1)计算全回转推进器的期望推进角α:
2)优化α
3)求各全回转推进器推进角
本实例中两个相邻的全回转推进器3和4不满足合理间距d3,4>2(D3+D4),引入推进器偏差角取Δα3,4=30°,则:
4)求取结构矩阵A(α)
由附图1得到本实例的结构矩阵A(α):
5)伪逆算法求推力分配
权值矩阵W采用单位矩阵,由
f=(W+WT)-1A(α)T(A(α)(W+WT)-1A(α)T)-1τ
6)推力限幅
7)结束
在VC编译环境下,可以利用C语言对上述方法在船舶中实现。
在Matlab环境下对算法进行仿真,给出期望的推力和力矩τ=[Fu,Fv,Nz]T,通过本发明所采用的方法,将力和力矩有效的分配到各个推进器。从图4图5可以看出,①当期望推进器推力(矩)平稳变化时,各推进器推进角和推力均可以在能力范围内平稳快速跟踪给定信号,能够有效避免奇异结构的发生。②在t=240附近,虽然期望推力和力矩突变较大,各推进器推进角变化在推进器能力(30o/s)范围内,依然可以快速高效的实现推力期望τ=[Fu,Fv,Nz]T。事实证明本发明中的推力分配算法可以实现1#~4#各个推进器的推力f和3#~4#全回转推进器推进角α的最佳分配,达到期望τ。
本发明提出了以伪逆算法为基础,解决推力禁区角、推力限制、推进器能力分配、推进器推进角变化率有限等问题的,一种高效率高速度的推力分配方法。从计算过程看,本发明所采用的方法避免了大量的迭代计算,节约计算时间,有效提高船舶动力定位系统的响应速度。另外,从计算结果可以看出,本发明所采用的方法能够高效的将纵荡、横荡和回转力矩分配给各个推进器,推进器推力变化平稳,推进角位置合理,能够有效避免奇异结构的发生。
以上仅为本发明的优选实施例,但本发明的设计构思并不局限于此,凡利用此构思对本发明做出的非实质性修改,也均落入本发明的保护范围之内。
Claims (4)
1.一种船舶及石油平台动力定位系统的推力分配方法,其特征在于包括如下步骤:
步骤一:推进器推进角的确定;
步骤二:求取结构矩阵;
步骤三:伪逆算法求推力分配;
步骤四:推力限幅确定;
上述步骤一推进器推进角的确定方法如下:
1)计算全回转推进器的期望推进方向的方向角α
取船舶所受合力的方向角作为全回转推进器期望推进方向的方向角α,以下称推进器推进角,则全回转推进器推进角α为:
定义τ=[Fu,Fv,Nz]T,Fu,Fv,Nz分别为船舶在纵荡、横荡和艏摇三个方向上的所受合力(矩),α以船头正前方为零度,顺时针增加,一周为360°;Fu,Fv,Nz的正方向分别为船头正前方,船体正右侧,顺时针方向;
2)全回转推进器推进角α的优化
本发明中的推进器推进角度的计算方法和伪逆算法,能实现对连续变化的给定信号的快速跟踪,如果期望推进器推进角α和当前时刻α'之差满足下式:
取
代替α;Δαmax为推进器偏差角最大值;
3)确定两个相邻的全回转推进器推进角
如果两个相邻的全回转推进器i和j不满足合理间距di,j>2(Di+Dj),引入推进器偏差角Δαi,j,用于解决两个相近的全回转螺旋桨之间产生的推力干扰问题,Di、Dj分别为第i、第j个推进器的直径,则推进器推进角满足:
其中αi、αj分别为第i、第j个推进器的推进角;推进器偏差角Δαi,j满足:
2.根据权利要求1所述的船舶及石油平台动力定位系统的推力分配方法,其特征在于上述步骤二求取结构矩阵的方法是:
定义τ=A(α)f,f=[f1,...,fp,fp-1,...,fn]T,fi表示第i个推进器的推力;A为结构矩阵,由推进器的类型、位置分布以及推进器角度得到:
其中,Lix,Liy分别代表第i个推进器在纵荡方向和横荡方向的力臂。
3.根据权利要求1所述的船舶及石油平台动力定位系统的推力分配方法,其特征在于上述步骤三伪逆算法求推力分配的方法是:
根据优化的目标函数:以及约束条件:S.t.τ=A(α)·f,求解推力f,其中W为权重矩阵,定义推力的拉格朗日函数为:
L(f,λ)=fTWf+λT(τ-A(α)f) (1.6)
推力的拉格朗日函数对变量f求导:
由局部最优定理可知:
f=(W+WT)-1A(α)Tλ (1.9)
(W+WT)T>0时,u取最小值,其中W+WT为奇异矩阵,且
Wii>0(i=1...n),由边界条件可知:
τ=A(α)f=A(α)(W+WT)-1A(α)Tλ
(2.0)从而得到:
λ=(A(α)(W+WT)-1A(α)T)-1τ
(2.1)
代入上式,可得推力:
f=(W+WT)-1A(α)T(A(α)(W+WT)-1A(α)T)-1τ
(2.2)
如果A(α)-1存在,由式(2.2),得
f=(W+WT)-1A(α)T(A(α)T)-1(W+WT)A(α)-1τ=A(α)-1τ
(2.3)
即,当A(α)奇异矩阵时,(2.2)式可实现推力的逆分配,然而结构矩阵A(α)一般为非奇异矩阵,因此我们利用W+WT实现推力的伪逆分配;
当W≠I(I为单位矩阵),(2.2)式为带权重的伪逆推力分配算法;当W=I,为等价或者无权重的伪逆推力分配算法。
4.根据权利要求1所述的船舶及石油平台动力定位系统的推力分配方法,其特征在于上述步骤四推力限幅的确定方法是:判断推进器推力f是否超过边界值,如果有推进器推力超出此推进器推力的边界值,则取此边界值作为此推进器的输出,实现推力分配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410348488.7A CN104133374B (zh) | 2014-05-23 | 2014-07-21 | 一种船舶及石油平台动力定位系统的推力分配方法 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2014102231744 | 2014-05-23 | ||
CN201410223174.4 | 2014-05-23 | ||
CN201410223174 | 2014-05-23 | ||
CN201410348488.7A CN104133374B (zh) | 2014-05-23 | 2014-07-21 | 一种船舶及石油平台动力定位系统的推力分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104133374A CN104133374A (zh) | 2014-11-05 |
CN104133374B true CN104133374B (zh) | 2017-01-25 |
Family
ID=51806095
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410348488.7A Expired - Fee Related CN104133374B (zh) | 2014-05-23 | 2014-07-21 | 一种船舶及石油平台动力定位系统的推力分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104133374B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301963B (zh) * | 2015-11-17 | 2018-04-06 | 江苏科技大学 | 一种基于船舶功率管理系统的推力优化分配方法 |
CN105808837B (zh) * | 2016-03-04 | 2019-06-07 | 武汉理工大学 | 面向船舶动力定位系统控制力平缓变化时推进器水动力干扰处理方法 |
CN108333926B (zh) * | 2018-01-19 | 2020-09-04 | 武汉理工大学 | 动力定位船舶的固定角度推力分配方法 |
CN108490770A (zh) * | 2018-02-28 | 2018-09-04 | 哈尔滨工程大学 | 一种基于混合算法的船舶动力定位系统推力分配方法 |
CN108628171B (zh) * | 2018-06-15 | 2021-03-30 | 哈尔滨工程大学 | 一种用于船舶动力定位系统的改进推力分配方法 |
CN109085756B (zh) * | 2018-08-27 | 2020-11-10 | 西安交通大学 | 一种基于遗传算法优化的水下机器人推力分配方法及系统 |
WO2020252722A1 (zh) * | 2019-06-20 | 2020-12-24 | 陈勇贵 | 一种新型的水翼船及其控制方法 |
CN111061285B (zh) * | 2019-12-12 | 2022-08-02 | 哈尔滨工程大学 | 一种船舶动力定位推力分配方法 |
CN113075884B (zh) * | 2021-03-29 | 2022-07-15 | 哈尔滨工程大学 | 基于自适应遗传-最小二乘互联预测系统的推力分配方法 |
CN115994456B (zh) * | 2023-03-21 | 2023-06-06 | 中国船舶集团有限公司第七〇七研究所 | 一种用于船舶推力分配的全回转推进器推力衰减处理方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102385665A (zh) * | 2011-08-30 | 2012-03-21 | 无锡中讯科技有限公司 | 采用遗传算法的船舶动力定位系统推力分配方法 |
CN102508431A (zh) * | 2011-11-04 | 2012-06-20 | 江苏科技大学 | 一种海洋钻井平台动力定位系统推力分配方法 |
-
2014
- 2014-07-21 CN CN201410348488.7A patent/CN104133374B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102385665A (zh) * | 2011-08-30 | 2012-03-21 | 无锡中讯科技有限公司 | 采用遗传算法的船舶动力定位系统推力分配方法 |
CN102508431A (zh) * | 2011-11-04 | 2012-06-20 | 江苏科技大学 | 一种海洋钻井平台动力定位系统推力分配方法 |
Non-Patent Citations (2)
Title |
---|
A Four-Quadrant Thrust Estimation Scheme Based on Chebyshev Fit and Experiment of Ship Model;Baoyu Ye,等;《The Open Mechanical Engineering Journal》;20121231;第6卷;第148-154页 * |
船舶动力定位系统推力分配策略研究;王钦若,等;《控制工程》;20130131;第20卷(第1期);第30-33页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104133374A (zh) | 2014-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104133374B (zh) | 一种船舶及石油平台动力定位系统的推力分配方法 | |
CN105446821B (zh) | 一种基于改进神经网络的智能水下机器人推进器故障诊断方法 | |
CN101408772B (zh) | Auv智能避碰方法 | |
Sutulo et al. | Mathematical models for simulation of manoeuvring performance of ships | |
CN106773722B (zh) | 一种基于人工鱼群算法的船舶动力定位系统推力分配方法 | |
CN103092077B (zh) | 动力定位系统的推力分配方法 | |
CN106527133A (zh) | 一种船舶多桨协调控制分配方法 | |
Li et al. | A study on the performance of cascade hard sails and sail-equipped vessels | |
Sanada et al. | ONR Tumblehome course keeping and maneuvering in calm water and waves | |
CN108845576A (zh) | 一种基于粒子群与序列二次规划结合的推力分配方法 | |
CN101859147A (zh) | 船舶航向智能协调控制方法 | |
US20210223747A1 (en) | Dynamic positioning and thrust distribution device and method based on artificial neural network | |
CN102841970B (zh) | 采用动态禁止角的动力定位推力分配装置及其分配方法 | |
Lim et al. | Electric power consumption predictive modeling of an electric propulsion ship considering the marine environment | |
Ueno et al. | Rudder effectiveness and speed correction in practice at tank test | |
Kosleck | Prediction of wave-structure interaction by advanced wave field forecast | |
Bańka et al. | Design of a multivariable neural controller for control of a nonlinear MIMO plant | |
Koenhardono et al. | Neural network-based engine propeller matching (NN-EPM) for trimaran patrol ship | |
Park et al. | A study on the key performance indicator of the dynamic positioning system | |
Jürgens et al. | Influence of thruster response time on dp capability by time-domain simulations | |
CN113126503B (zh) | 一种考虑推进器故障的智能船舶推力分配方法 | |
Luo et al. | Reviews on the power management for shipboard energy storage systems | |
Ruth et al. | Thrust allocation with linear constrained quadratic cost function | |
Cavallo et al. | A Model for Simulating Unsteady Wake Interference Effects in Aerial Refueling | |
Zaghi et al. | Rans hydrodynamic characterization of a usv swath configuration including design optimization |
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: 20170125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |