一种船舶及石油平台动力定位系统的推力分配方法.docx

一种船舶及石油平台动力定位系统的推力分配方法一种船舶及石油平台动力定位系统的推力分配方法本发明是一种船舶及石油平台动力定位系统的推力分配方法包括如下步骤：步骤一：推进器推进角的确定；步骤二：求取结构矩阵；步骤三：伪逆算法求推力分配；步骤四：推力限幅确定本发明实现了船舶及石油平台动力定位推力分配的快速性，稳定性，高效性，是一种方便实用的船舶及石油平台动力定位系统的推力分配方法专利说明】一种船舶及石油平台动力定位系统的推力分配方法【技术领域】 [0001]本发明是一种船舶及石油平台动力定位系统的推力分配方法，特别是一种水面船舶及石油平台动力定位系统的多螺旋桨推力分配方法，属于船舶及石油平台动力定位系统的推力分配方法的创新技术背景技术】 [0002]针对船舶/石油平台平面定位的需要，采用多组螺旋桨推力以消除船舶/石油平台上受到的海风、海浪、海流的影响，使船舶/石油平台等按预定的位置、航向、轨迹、速度保持姿态船舶/石油平台动力定位系统是一种闭环控制系统，它是通过推进器提供的推力，使海洋结构物在海上保持一定的位置和艏向角的一种先进的定位系统它的基本原理是通过获得海洋结构物的位置误差数据，由控制器计算出抗拒环境干扰所需的推力和转矩，并通过推力系统的实现，以此达到消除位置误差的目的。

此动力定位系统一般由控制系统、测量系统和推力系统组成推力系统是船舶动力定位系统的一个重要组成部分，作为船舶动力定位系统的执行机构，其作用是提供控制系统所需要的推力和力矩，用来满足定位的需要，推力系统主要由发电机组、输电设备、推进器等组成推进器一般来说是指螺旋桨，应用比较广泛的是槽道式推进器，全回转推进器和吊舱推进器 [0003]伪逆算法是动力定位工程应用中普遍采用的方法，这种算法无需多次进行迭代计算，方法简单，但无法解决推进角及禁区角问题现有的技术当中，有利用蚁群算法、遗传算法、二次序列规划算法等来解决推力分配最优解问题，但这些算法需要多次迭代，计算时间长，不易用于船舶/石油平台的实时操作控制中发明内容】 [0004]本发明的目的在于考虑上述问题而提供一种实现了船舶及石油平台动力定位推力分配的快速性，稳定性，高效性的船舶及石油平台动力定位系统的推力分配方法 [0005]本发明的技术方案是:本发明的船舶及石油平台动力定位系统的推力分配方法，包括如下步骤:步骤一:推进器推进角的确定；步骤二:求取结构矩阵；步骤三:伪逆算法求推力分配；步骤四:推力限幅确定 [0006]本发明针对可调节的带方向角约束的多螺旋桨的推进角度计算问题和各螺旋桨的推力分配问题，即带约束的多参数的快速最优解的求解问题，本发明提出先分配到推进角方法，实现推力方向角的快速跟踪和推力的高效性；再分配到推力的方法，采用伪逆算法实现不同推进能力的螺旋桨的协调工作，有效的避免了相邻近的推进器之间推力干扰引起的推力衰减问题，解决了其他算法中的多奇异值问题，多次迭代计算问题，克服了收敛时间长的缺点，实现了船舶/石油平台动力定位推力分配计算的快速性，稳定性，高效性，易于在船舶、石油平台上实现。

本发明是一种方便实用的船舶及石油平台动力定位系统的推力分配方法专利附图】【附图说明】 [0007]图1为船舶/石油平台推进器分布与推进角度示意图； [0008]图2为船舶/石油平台推进器间距与偏差角关系图； [0009]图3为船舶/石油平台动力定位系统推力分配程序框图； [0010]图4为推力分配结果对比仿真图； [0011]图5为各螺旋桨推力分配结果仿真图具体实施方式】 [0012]实施例: [0013]本发明船舶及石油平台动力定位系统的推力分配方法，包括如下步骤: [0014]步骤一:推进器推进角的确定； [0015]步骤二:求取结构矩阵； [0016]步骤三:伪逆算法求推力分配； [0017]步骤四:推力限幅确定 [0018]上述步骤一推进器推进角的确定方法如下: [0019]I)计算全回转推进器的期望推进方向的方向角α [0020]取船舶所受合力的方向角作为全回转推进器期望推进方向的方向角α，以下称推进器推进角，则全回转推进器推进角α为:【权利要求】1.一种船舶及石油平台动力定位系统的推力分配方法，其特征在于包括如下步骤: 步骤一:推进器推进角的确定； 步骤二:求取结构矩阵； 步骤三:伪逆算法求推力分配； 步骤四:推力限幅确定。

2.根据权利要求1所述的船舶及石油平台动力定位系统的推力分配方法，其特征在于上述步骤一推进器推进角的确定方法如下: 1)计算全回转推进器的期望推进方向的方向角α 取船舶所受合力的方向角作为全回转推进器期望推进方向的方向角α，以下称推进器推进角，则全回转推进器推进角α为:定义τ = [FU，FV，NZ]T，FU，FV，NZ分别为船舶在纵荡、横荡和艏摇三个方向上的所受合力(矩)，α以船头正前方为零度，顺时针增加，一周为360° ；FU, Fv, Nz的正方向分别为船头正前方，船体正右侧，顺时针方向； 2)全回转推进器推进角α的优化 本发明中的推进器推进角度的计算方法和伪逆算法，能实现对连续变化的给定信号的快速跟踪，如果期望推进器推进角α和当前时刻α’之差满足下式:代替α ； 3)确定两个相邻的全回转推进器推进角 如果两个相邻的全回转推进器i和j不满足合理间距屯」>2 (D.+Dj)，引入推进器偏差角Δ a u，用于解决两个相近的全回转螺旋桨之间产生的推力干扰问题，Di^Dj分别为第1、第j个推进器的直径，则推进器推进角满足:其中α 1、α j分别为第1、第j个推进器的推进角；推进器偏差角Λ a 满足:3.根据权利要求1所述的船舶及石油平台动力定位系统的推力分配方法，其特征在于上述步骤二求取结构矩阵的方法是: 定义τ = Α(α )f，f = [fp…，fp, ，…，表示第i个推进器的推力；A为结构矩阵，由推进器的类型、位置分布以及推进器角度得到:其中，Lix, Liy分别代表第i个推进器在纵荡方向和横荡方向的力臂。

4.根据权利要求1所述的船舶及石油平台动力定位系统的推力分配方法，其特征在于上述步骤三伪逆算法求推力分配的方法是:根据优化的目标函数:mjnw = n|n.m ?以及约束条件:S.t.τ =Α(α).?.，求解推力f，其中W为权重矩阵，定义推力的拉格朗日函数为: L (f, λ ) = fTfff + λ T ( τ -A ( a ) f) (1.6) 推力的拉格朗日函数对变量f求导:由局部最优定理可知: f = (ff+ffT) -1A ( a )Τ λ (1.9) (W+WT)T>0时，u取最小值，其中W+WT为奇异矩阵，且Wn>0 (i = 1...η)，由边界条件可知: τ = A(a )f = A(a ) (ff+ffT) -1A ( a )Τ λ (2.0) 从而得到: λ = (A(a) (ff+ffT)^A(a)T)_1 τ (2.1) 代入上式，可得推力: f = (ff+ffT) ^1A ( a )T (A ( a ) (ff+ffT) ^1A ( a )T) -1 τ (2.2) 如果Α(α 存在，由式(2.2)，得 f = (ff+ffT)_1A(a )T(A(a )T)^ (ff+ffT)A(a )^1 τ = A ( a ) τ (2.3) 即，当Α(α)奇异矩阵时，(2.2)式可实现推力的逆分配，然而结构矩阵AU) —般为非奇异矩阵，因此我们利用W+WT实现推力的伪逆分配； 当W古I (I为单位矩阵)，(2.2)式为带权重的伪逆推力分配算法；当W = I，为等价或者无权重的伪逆推力分配算法。

5.根据权利要求1所述的船舶及石油平台动力定位系统的推力分配方法，其特征在于上述步骤四推力限幅的确定方法是:判断推进器推力f是否超过边界值，如果有推进器推力超出此推进器推力的边界值，则取此边界值作为此推进器的输出，实现推力分配。