船舶纵向运动多输入多输出不稳定系统的非线性反馈控制
Nonlinear Feedback Control of MIMO Unstable System for Ship Longitudinal Motion
通讯作者: 张显库,教授,博士生导师,电话(Tel.): 0411-84729572;E-mail:zhangxk@dlmu.edu.com.
责任编辑: 王一凡
收稿日期: 2022-03-22 修回日期: 2022-07-23 接受日期: 2022-07-27
基金资助: |
|
Received: 2022-03-22 Revised: 2022-07-23 Accepted: 2022-07-27
作者简介 About authors
曹添(1997-),硕士生,主要从事船舶运动控制方面的研究.
为了提高多输入多输出(MIMO)不稳定系统控制器的鲁棒性与节能效果,以大连海事大学教学实习船“育鲲”轮在波浪上纵摇和升沉运动的数学模型为研究对象,将数学模型的分子分母具有不稳定零极点过程使用根轨迹成形技术以及引用数学工具加权矩阵将其根轨迹拉到左半平面.之后运用闭环增益成形算法设计控制器的思路,以成形后的模型用一阶闭环增益成形算法设计鲁棒控制器,再引入双极性S函数驱动的非线性反馈代替原有的误差线性反馈.从仿真结果可以看出,在模型产生摄动时,非线性反馈控制改进了升沉和纵摇,控制能量大幅减少,对模型干扰具有鲁棒性能.该方案求解过程简单,并且整个系统控制器与被控对象物理意义明确,设计所得的控制器阶次较低,且补充完善了MIMO不稳定系统的控制器设计方法.
关键词:
In order to improve the robustness and energy saving of the controller for multiple input multiple output (MIMO) unstable systems, the mathematical model of pitching and the heave motion of the scientific ship “YUKUN” of Dalian Maritime University is taken as the research plant, and the root locus of the transfer function mathematical model with unstable zero-pole and zero-pole processes is pulled to the left half plane by using the root locus shaping technique and weighting matrix. The first-order closed-loop gain shaping algorithm is used to design the robust controller, and the nonlinear feedback driven by the bi-polar S function is introduced to replace the original error linear feedback. The simulation results show that the nonlinear feedback control improves heave and pitching when the model is perturbed, greatly reduces the control energy, and is robust to model disturbance. The solution process of the scheme is simple, and the physical significance of the controller and the controlled object of the whole system is clear. The order of the controller design is low, and the controller design method of the MIMO unstable system is supplemented and perfected.
Keywords:
本文引用格式
曹添, 张显库.
CAO Tian, ZHANG Xianku.
船舶运动控制是海洋工程的热点研究领域,是实现船舶自动导航的核心.船舶在海上航行时由于受到风浪流等环境因素、螺旋桨舵等船舶自身的各种推动和控制作用会产生复杂的运动.船舶的船体结构很复杂,与流体之间的作用也很复杂.并且,船舶在海上航行受到海浪等扰动的作用时,会产生航向与航迹的偏移以及船舶的摇摆.船舶的摇摆中,纵摇和横摇的影响尤为严重,而船舶纵摇控制与横摇控制是船舶运动控制的重点组成部分.
不稳定系统的控制问题在工业发展中占据重要位置,而不稳定系统控制器的设计至关重要.目前针对不稳定系统的控制研究,主要研究内容有不稳定系统的镇定、鲁棒控制器设计以及系统模型不稳定过程的控制实例仿真[7].主要思路是把不稳定的控制过程转换为稳定的控制过程.其中文献[8]中采用经典的内模控制法,能够将不稳定的过程转换为稳定的控制过程;文献[9]中也针对单输入单输出(SISO)不稳定系统的鲁棒控制器设计给出了系列研究成果.此外,船舶航向保持也是一个典型的SISO问题, 在文献[10]中提出了一种基于非线性正弦函数修饰的闭环增益成形算法的鲁棒积分反推控制方法,仿真结果表明控制器具有合理的操舵频率和良好的节能效果.文献[11]中基于二阶闭环增益成形算法设计出具有强鲁棒性的线性控制器,并使用双极性S函数对线性控制律进行修饰,进一步增强了系统抗高频干扰能力和节能效果.但是对于多输入多输出(MIMO)的不稳定系统,前述控制方法的控制效果并不好.文献[12]中以水翼艇的纵向运动数学为研究对象,对MIMO鲁棒控制器进行改善,提高了鲁棒性和节能效果,但其研究传递函数数学模型具有稳定的零极点,因此其控制器设计方法不适于具有不稳定极点的数学模型.而“育鲲”轮的纵摇模型属于MIMO不稳定系统,针对其鲁棒控制器的设计尤其是具有一对靠近虚轴不稳定极点的情况,还没有相关研究结果报道.
基于上述分析,提出一种具有良好减纵摇效果的简捷鲁棒控制器,并通过理论分析和多种条件下的仿真实验对其有效性进行验证.首先给出了“育鲲”轮纵摇严重的机理解释,并针对MIMO不稳定系统控制器设计难点提出一种新思路,最后使用双极性S函数非线性反馈技术代替线性反馈,进一步提升系统的鲁棒性与节能效果.
1 船舶模型
船舶航行时受到风浪等环境因素的影响,同时由于螺旋桨舵等船舶自身的各种推动和控制作用,会产生复杂的运动.船舶的船体结构与外形比较复杂,因此在水中运动时船舶和流体之间的作用也比较复杂,且船舶在海上航行时受到风浪等扰动的影响时会产生航向与航迹的偏移,其中船舶的横摇和纵摇更为严重.以牛顿动力学定理为基础,结合船舶运动控制的特点分析讨论船舶纵向运动方程的建立.图1所示为“育鲲”轮运动坐标系.
图1
船舶运动坐标系$G_0 x_b y_b z_b$,坐标系原点取在船的重心$G_0$上,坐标轴与船体中心的惯性主轴重合. $G_0 y_b$轴在船的中线面内且平行于基线,指向船首为正; $G_0 x_b$轴垂直于中线面,指向右舷为正; $G_0 z_b$轴在船的中线面内且垂直于$G_0 x_b$轴,向上为正.结合坐标系,船舶在任意时刻的运动可以分解为沿着3个坐标轴的直线运动与绕3个坐标轴的转动,因此形成往复运动分别为横摇$φ$、纵摇$θ$和艏摇$ψ$.
参考文献[13]对船舶在波浪中的运动特征的描述,得出船舶纵摇运动的数学模型如下:
式中:$h$为船舶的升沉位移; $θ$为纵摇角度; $k_hh$、$k_{hh_1}$ 、$k_{hθ_1}$ 、$k_{hθ_0}$、$k_{θθ_1}$、$k_{θθ_0 }$、$k_{θh_1}$、$k_{θh_0}$为流体动力系数; $F_h (t) $为深沉运动的波浪干扰力, $t$为时间; $M_θ (\dot{t})为纵摇运动的波浪干扰力矩;控制向量$U=[FM]^T$,分别为控制船舶升沉和纵摇外加的力和力矩.为了更直观地分析MIMO不稳定系统成形传递函数的设计,将波浪干扰力与波浪干扰力矩等效为一种简单的线性模型,即白噪声干扰.
为简化数学模型设计,做出如下基本假定:
(1) 假定船体是刚体,忽略它的弹性变性.
(2) 对于船舶减纵摇控制系统,忽略纵荡、横荡、首摇和横摇等船舶运动,只考虑纵摇与垂荡.
如果输出h和θ,则传递函数型数学模型G为
其中:s为拉普拉斯算子;i为虚数单位.
表1 “育鲲”轮参数
Tab.1
参数 | 满载取值 | 参数 | 满载取值 |
---|---|---|---|
两柱间长L/m | 105 | 航速v/kn | 16.7 |
船宽B/m | 18 | 排水体积Δ/m3 | 5735.5 |
重心距船舯距离xc | -0.51 | 船舯吃水dM/m | 5.2 |
舵叶面积AR/m2 | 11.46 | 方形系数Cb | 0.5595 |
2 研究方法
2.1 “育鲲”轮纵摇严重的机理分析
分析“育鲲”轮的纵摇数学模型的各项参数, G1中有一对虚根在右半平面,不稳定,且非常靠近虚轴.并且a11的零点也非常靠近虚轴,系统模型很容易发散.船舶的纵向运动模型是MIMO的多变量系统,尤其输入是由不规则波浪产生的力与力矩引起的纵摇,情况较为复杂.船舶纵摇模型的传递函数根轨迹如图2所示.
图2
图2
系统开环传递函数的根轨迹
Fig.2
Root trajectories of open-loop transfer functions of the system
由上图可以分析得出$G_{12} $、$G_{21}$、$G_{22}$的根都在右半平面,左半平面的极点也非常靠近虚轴,且$G_{12}$和$G_{21}$的根轨迹都是发散的,不稳定.由此可知,“育鲲”轮的纵向运动数学模型极易不稳定,这也是纵摇严重的机理.
2.2 不稳定系统的根轨迹成形原理
已知根轨迹理论是在复平面上由开环零极点的分布状况来研究闭环特征根的轨迹随开环增益的变化而变化的理论,不同开环增益决定了在根轨迹上有不同位置的闭环极点及相应的系统动态性能.控制系统的稳定性完全由其特征方程的根决定,而特征方程的根又与系统参数紧密相关,可见根轨迹成形原理就是基于变增益的概念而来.
根据上述原理,取被控对象G(s)根轨迹成形的权函数为
且设计权函数系数矩阵为λ,取a=50.
2.3 设计权函数参数并验证其稳定性
根据权函数成形原则,取
L11与 L22包含两个零点和一个极点,且再乘以权函数系数矩阵得到的广义被控对象P为
将根轨迹增益变小,乘以系数0.05,成形后得出根轨迹如图3所示.图中
图3
图3
成形后的广义被控对象根轨迹图
Fig.3
Root locus diagram of generalized controlled plant after forming
2.4 MIMO不稳定系统的鲁棒控制器设计
船舶减纵摇控制图如图4所示.图中:r为白噪声干扰;u为控制器输出信号;G'为不稳定系统成形后的稳态输出值.
图4
基于船舶在航行中受波浪的作用而产生纵摇和升沉运动,使用白噪声模拟波浪干扰输入,控制器K的原理为运用闭环增益成形算法设计控制器的思路,以简单的积分环节与不稳定系统成形后的稳态输出值G'构成的闭环传递函数,从而避免了高阶控制器鲁棒性差的问题.其中积分环节本身就可以消除持续恒定干扰和慢时变干扰对系统性能的影响,而且大大减小控制输出切换的频率和幅值,减少对成形后的广义被控对象P的影响,并且能有效减少波浪对船舶纵摇角和沉深位移的影响.
对于一个MIMO系统,在设计控制器K时,令其非对角线元素为0,则K阵的形式如下:
设A、B为成形后广义被控对象单位阶跃信号激励后的稳定输出值,进一步推导得出其传递函数为
最终化简得到:
式(10)为一阶闭环增益成形算法的情况[14], 属于一种简单的鲁棒控制算法.
根据上式推导得出$AC$=1, $BD$=1,所以可得出$K_{11}=\frac{C}{T_{11}s}$,$K_{22}=\frac{D}{T_{22}s}$K_11=C/(T_11 s),K_22=D/(T_22 s), $T_{11}$和$T_{22}$为控制器参数.整个控制系统等价于一个一阶闭环增益成形算法.闭环增益成形算法是利用$H_∞$鲁棒控制理论的结果使用具有工程意义的4个参数构造出闭环系统传递函数,进而反推出鲁棒控制器的算法[14].
2.5 双极性S函数非线性反馈改进
图5
图5
双极性S函数驱动的非线性反馈结构图
Fig.5
Nonlinear feedback structure driven by bipolar S-function
3 仿真实验和结果分析
为了验证所设计成形后广义被控对象的稳定性与控制器的有效性,在MATLAB中用Simulink进行仿真实验.
3.1 基于MIMO根轨迹成形模型稳定性实验验证
图6
图7
经过上图用Simulink搭建模块实验验证之后可知,使用根轨迹成形法以及乘以对称矩阵的设计思路是有效的.
3.2 基于船舶减纵摇鲁棒控制器性能验证
图8
图8
未加入非线性反馈函数比较控制效果图
Fig.8
Comparison control effect without adding nonlinear feedback function
图9
图9
加入非线性反馈函数比较控制效果图
Fig.9
Comparison control effect of adding nonlinear feedback function
为了定量比较控制效果,定义总升沉位移误差eTDDE与总纵摇角度误差eTPAE两个指标进行分析.定量分析结果如表2所示.
式中:$h(t) $为$t$时刻所对应的升沉位移; $θ(t) $为$t$时刻所对应的纵摇角度.
表2 闭环性能比较
Tab.2
控制方法 | θmax/(°) | hmax/m | eTDDE | eTPAE |
---|---|---|---|---|
被控之前 | 5.21 | 3.94 | 11.72 | 15.48 |
加入鲁棒控制器 | 2.43 | 0.33 | 1.44 | 9.73 |
非线性反馈改进 鲁棒控制器 | 1.09 | 0.16 | 0.78 | 4.35 |
为了进一步验证控制器的鲁棒性,给模型加入摄动,相当于在原模型上加入一个数值为0.15的纯滞后环节e-0.15s,可得摄动模型P'(s).图10给出了模型产生摄动时的仿真结果.
图10
图10
模型产生摄动时的仿真结果
Fig.10
Simulation results when the model generates perturbations
通过上述实验比较可知,对广义模型加入鲁棒控制器之后很明显地改进了船舶的升沉和纵摇效果,升沉位移最大值降为被控前最大值的8.0%,纵摇角最大值降为被控前最大值的46.6%,且综合对比eTDDE与eTPAE,分别改进了87.71%和37.14%.对控制器进行非线性反馈改进之后,升沉位移的最大值降为被控之前的4.0%,纵摇角最大值降为被控前的19.1%,对比eTDDE与eTPAE,分别改进93.4%与71.9%.
模型摄动时,鲁棒控制的升沉位移最大值为0.23 m,纵摇角的最大值为1.19°,且eTDDE与eTPAE分别为1.02与4.02.从表2中可以看出其控制效果基本与广义模型基本相同,说明控制器具有一定的鲁棒性.仿真结果表明,所设计的鲁棒控制器简捷有效,控制效果更优,而非线性反馈改进的鲁棒控制器具有明显的节能效果,补充完善了MIMO不稳定过程的控制器设计方法.
进一步从控制器的角度来分析,控制器K属于2×2对角矩阵形式的一阶惯性系统,在海上航行时高频的波浪干扰是引起船舶升沉和纵摇变化的主要原因, 而一阶鲁棒控制器的引入有效地抑制了高频干扰,使系统有很强的鲁棒性.
4 结语
基于MIMO不稳定系统提出了一种船舶减纵摇的新型节能控制方案.该方案首先以大连海事大学教学实习船“育鲲”轮的纵摇数学模型为研究对象,对“育鲲”轮纵摇严重的机理进行了分析.随后对系统模型使用根轨迹成形法,以及引用数学工具加权矩阵L(s)使不稳定系统的模型趋于稳定.控制器的设计思路是闭环增益成形算法的一种衍生,并且使用双极性S函数进行非线性改进.通过Simulink工具箱验证了在一般海况下控制器对高频波浪的抑制作用,且效果良好.之后对成形后广义被控对象模型加入非线性反馈改进的控制器进行仿真实验,并与原有控制器对广义被控对象的控制效果对比.实验结果表明,改进后控制器的模型最大升沉位移降为被控之前的4.0%,最大纵摇角降为被控之前的19.1%,且eTDDE与eTPAE分别改进了93.4%与71.9%.随后给出了模型产生摄动时的仿真结果,基本与广义模型相同,表明本文所设计控制器在一般海况下具有较好的减纵摇效果,并且非线性反馈改进的控制器具有明显的节能效果.最后对本文控制器具有良好的控制效果的原因给出了理论分析.本文所设计研究方案简单有效且节能,弥补了MIMO不稳定系统控制器设计的空白.在未来的研究中,将结合本文的控制算法设置特定的纵倾角控制装置,并将该装置运用到船舶自身的动力控制系统当中,产生抵抗纵摇的阻尼力矩以及恢复力矩以达到减缓纵向运动的目的.
参考文献
大型油轮艏摇混沌现象的仿真与滑模控制
[J]. ,
Modeling and sliding mode control for chaotic yawing phenomenon of large oil tanker
[J]. ,
Robust integral backstepping and terminal synergetic control of course keeping for ships
[J]. ,DOI:10.1016/j.oceaneng.2020.108532 URL [本文引用: 1]
Global robust adaptive path-tracking control of underactuated ships under stochastic disturbances
[J]. ,DOI:10.1016/j.oceaneng.2015.10.038 URL [本文引用: 1]
反步鲁棒镇定及其在舵减横摇中的应用
[J]. ,
Back-steeping robust stabilization and its application in rudder roll stabilization
[J]. ,
船舶纵摇-升沉运动系统的解耦
[J]. ,
The decoupling of ship pitching and heaving motion
[J]. ,
船舶在波浪上纵向运动与控制研究综述
[J]. ,
Review of longitudinal motion and control of ships on waves
[J]. ,
Control of unstable delayed recycling systems
[J]. ,DOI:10.1002/apj.v14.2 URL [本文引用: 1]
结构不稳定系统的智能控制方法及其应用
[J]. ,
Intelligient control method of structurally unstable system and its application
[J]. .
Robust control of pressure for LNG carrier cargo handling system via mirror-mapping approach
[J]. ,
Improvement of integrator backstepping control for ships with concise robust control and nonlinear decoration
[J]. ,DOI:10.1016/j.oceaneng.2019.106349 URL [本文引用: 1]
Energy saving of course keeping for ships using CGSA and nonlinear decoration
[J]. ,DOI:10.1109/Access.6287639 URL [本文引用: 2]
基于非线性反馈的水翼艇纵向运动鲁棒控制
[J]. ,
Robust control of longitudinal motion for hydrofoils based on nonlinear feedback
[J]. .
Vertical plane motion control of an S-SWATH vehicle with flapping foil stabilisers sailing in waves
[J]. ,DOI:10.1016/j.oceaneng.2016.05.004 URL [本文引用: 2]
Linear reduction of backstepping algorithm based on nonlinear decoration for ship course-keeping control system
[J]. ,DOI:10.1016/j.oceaneng.2017.10.017 URL [本文引用: 1]
/
〈 | 〉 |