基于LQR限寬水域中KVLCC2操縱運動控制研究

邵 闖，馬 寧,2，顧解忡,2

(1. 上海交通大學 船舶海洋與建筑工程學院，上海 200240；2. 上海交通大學 海洋工程國家重點實驗室，上海 200240)

基于LQR限寬水域中KVLCC2操縱運動控制研究

邵 闖1，馬 寧1,2，顧解忡1,2

(1. 上海交通大學 船舶海洋與建筑工程學院，上海 200240；2. 上海交通大學 海洋工程國家重點實驗室，上海 200240)

在航道寬度受限制的水域中，船舶會受到岸壁效應的影響，橫向力與首搖力矩將發生變化，這會對船舶的航行安全產生不利的影響。鑒于此問題，本文應用現代控制理論最優控制LQR方法，對在限制水域中航行的超大型油輪KVLCC2的操縱運動進行控制研究。為便于LQR控制器的設計，采用線性狀態空間形式的操縱運動方程，基于數值模擬獲取的相應線性水動力系數，計算出使目標函數值最小的增益矩陣K，從而得到滿足最優控制規律的時域舵角變化，實現對不同寬度水域中船舶運動的最優控制，并與極點配置控制法作比較，驗證LQR控制器的優越性。結果表明，當船岸距離d/L≥=1.2時，船舶基本不受岸壁效應的影響，控制幅度極??；當岸壁距離d/L=0.25時，擺舵角度將超過6°，同時船舶前進速度也將下降，下降幅度將超過前進速度的10%，岸壁效應明顯。

LQR；限寬水域；最優控制

0 引 言

隨著船舶的大型化，海峽通道、運河航道等船舶航行的水域寬度相對變窄，岸壁效應隨之產生，船舶的水動力性能相較于在無限寬水域中航行時會發生顯著變化，這也造成限制水域中船舶操縱運動的控制愈加困難。限寬水域船舶操縱運動中，船舶橫蕩與首搖運動之間具有較強的耦合性，這需要尋求多輸入多輸出（MIMO）方法來控制船舶運動。船舶控制的主要方法有PID控制、魯棒控制和模糊控制[1 – 3]。但在MIMO控制領域，上述方法控制效果一般且算法復雜。Mucha P 與Thomas[4 – 5]則利用LQR控制方法實現對限寬水域中船舶操縱運動的較好控制，展現了LQR控制方法在MIMO控制領域較好的兼容性，為復雜狀況下船舶操縱運動的控制提出了一個新的思路。

LQR（Linear Quadratic Regulater）即線性二次型調節器，其控制對象是以狀態空間形式給出的線性系統，而目標函數為對象狀態和控制輸入的二次型函數。LQR的最優解具有標準的解析式，可以形成簡單的線性狀態反饋控制規律，容易構成最優反饋控制，較好地兼顧了系統的魯棒穩定性和快速性，在工程中便于實現。

本文以超大型油輪KVLCC2為研究對象，基于LQR控制方法對限寬水域中KVLCC2操縱運動控制問題進行研究。本文采用線性狀態空間形式的操縱運動方程，基于數值模擬獲取的相應線性水動力系數，計算出使目標函數值最小的增益矩陣K，從而得到滿足最優控制規律的時域舵角變化，實現對不同寬度水域中船舶運動的最優控制，并與極點配置控制方法作比較，結果驗證了LQR控制器的有效性和優越性。

1 數學模型

1.1 線性船舶操縱運動方程

設船舶排水量為?，z軸慣性矩為Iz，X，Y分別為船舶在縱向與橫向所受外力，N為船舶所受豎向力矩。u為船舶前進速度，v為船舶橫移速度，r為船舶轉首角速度，xg為船舶重心縱向坐標，ψ為船舶航向角，δR為轉舵角度。

基于LQR控制器的控制對象是以狀態空間形式給出的線性系統，本文采用首揺-橫蕩二自由度船舶線性操縱運動方程。非限制水域下線性船舶操縱運動方程為：

1.2 岸吸力與力矩

在限寬水域中，由于岸壁效應的影響，船舶將受到岸吸力與力矩的作用。這些岸吸力與力矩主要取決于岸壁距離與前進速度的大小。Thomas[5]岸吸力與力矩的計算公式為：

其中L，B和T分別為船長、船寬與吃水。

在式（1）的右端分別迭加岸吸力YS與岸吸力矩NS，可得修正后的船舶線性操縱運動方程

2 控制方法

2.1 LQR控制

線性二次型最優設計旨在尋求最優控制u（t），使二次型目標函數J取得最小值。

考慮系統的狀態方程為：

式中：x為系統的狀態向量；u為控制向量；D為系統的干擾矩陣；A為由系統運動方程推導而得的控制對象矩陣；B為控制矩陣，用來表示控制向量u對狀態向量x的線性作用。

而系統的二次型性能指標為：

式中：Q為對狀態向量的加權矩陣，是半正定對稱矩陣；R為對控制向量的加權矩陣，是正定對稱矩陣。這 2 個量決定了系統誤差與控制能量消耗之間的相對重要性。

為使J最小，由最小值原理得到最優控制為：

其中P為代數矩陣Riccati方程的解。

綜上，最優反饋系數矩陣K為：

2.2 極點配置法

極點配置[7 – 8]控制法是通過比例環節的反饋把定常線性系統的極點移置到預定位置的一種控制方法。極點配置的實質是用比例反饋去改變原系統的運動模式，以滿足設計規定的性能要求。

3 水動力系數

為了獲取船舶的水動力導數，本文基于計算流體力學（CFD）技術并開發了混合動網格技術，對限寬水域中的平面運動機構（PMM）試驗進行數值模擬，從而獲取不同岸壁距離下KVLCC2的水動力導數。劉晗等[9 – 10]將上述數值計算結果與循環水槽PMM試驗結果進行對比，兩者相似程度較高，表明本文所用數值計算水動力系數具有較好的精度。

本文選用 4 種岸壁距離下KVLCC2的水動力導數，岸壁距離分別為：d/L=0.25，d/L=0.35，d/L=0.50，d/L=1.20。圖2給出了岸吸力與矩隨岸壁距離變化的試驗結果?？梢钥吹?，隨著岸壁距離的增加，岸吸力與力矩快速下降。

4 狀態空間方程

4.1 線性狀態空間方程

其中：

4.2 積分反饋

當船舶航行于限制水域中，LQR控制器將產生一個穩定的舵角來抵消岸吸力與岸吸力矩。但當船舶達到穩定狀態時，船舶已經與原始航跡產生了偏移。為了消除這一穩態漂移偏差，將在狀態空間方程中加入一個狀態變量yi，

則狀態空間方程為：

為了展現積分反饋對于控制效果的影響，本文對有穩態誤差項方程和無穩態誤差項方程分別進行了LQR控制仿真模擬，結果如圖3所示?？梢钥闯?，無積分反饋項的方程只能保證航向控制，而有積分反饋項的方程可以同時保證航向與航跡的控制。

5 控制模擬分析

本文分別對不同岸壁距離下船舶控制進行了模擬，岸壁距離分別選取了 4 種情況：d/L=0.25，d/L=0.35，d/L=0.50，d/L=1.20。

5.1 LQR與極點配置的比較

為了驗證LQR控制較極點配置控制的優越性，本文分別利用這2種控制方法對岸壁距離d/L=0.25工況下的船舶運動控制進行了仿真。對于極點配置控制，主導極點的選取對控制效果的影響至關重要，它將決定系統的阻尼特性與穩定性，本文最終選定的極點為–1.1，–0.85，–0.35，–0.08+0.096*j，–0.08–0.096*j。

表 1 狀態變量與控制變量最大期望值Tab. 1 Maximum desired state and control variables

由圖4仿真結果可看出，對于船舶橫移的控制，LQR控制器相較于極點配置輸出明顯較小的最大橫移偏差，且回復至原有航跡更趨平緩；對于輸入舵角，均符合物理現實，但極點配置控制的舵角顯示了較大的振蕩，而LQR控制器的舵角輸入更趨平緩，操舵的執行性更好。綜上，相較極點配置控制，LQR控制器具有較好的魯棒穩定性，顯現了LQR控制器的的有效性和優越性。

5.2 LQR控制結果分析

圖5展現了不同岸壁距離下LQR控制仿真結果。結果表明，當岸壁距離d/L＞1.2時，船舶趨于穩定時，基本沒有發生橫移，且舵角基本上保持為 0，船舶基本不受岸壁效應的影響，控制幅度極?。浑S著岸壁距離的逐漸減小，船舶最后的穩定舵角與降速也越來越大，當岸壁距離d/L=0.25時，舵角將超過6°，速度下降幅度超過前進速度的10%，岸壁效應明顯。

6 結 語

本文基于線性狀態空間形式的船舶操縱運動方程，根據最優控制理論，設計了LQR控制器并對不同寬度水域中船舶運動進行了控制模擬，并與極點配置控制法作比較，結果驗證了LQR控制器的有效性和優越性。研究結果表明，隨著船舶逐漸靠近岸壁，船舶所受橫向力與首搖力矩將變大，使用本文設計的LQR控制器通過擺舵與降速可有效矯正與控制船舶運動軌跡，達到最優控制的目的。

[ 1 ]張顯庫. 船舶運動控制[M]. 北京: 國防工業出版社, 2006.

[ 2 ]王先洲. 船舶及潛艇操縱中的魯棒控制研究[D]. 武漢: 華中科技大學, 2006.

[ 3 ]李鐵山, 楊鹽生, 洪碧光, 等. 船舶航跡控制魯棒自適應模糊設計[J]. 控制理論與應用, 2007, 24(3): 445–448.

[ 4 ]MUCHA P, MOCTAR O E. Ship-bank interaction of a large tanker and related control problems[C]// ASME 2013, International Conference on Ocean, Offshore and Arctic Engineering. 2013: V005T06A072.

[ 5 ]THOMAS B S, SCLAVOUNOS P D. Optimal-control theory applied to ship maneuvering in restricted waters[J]. Journal of Engineering Mathematics, 2007, 58(1–4): 301–315.

[ 6 ]BRYSON A E, HO Y C, SIOURIS G M. Applied optimal control: optimization, estimation, and control[M]. distributed by Halsted Press, 1975.

[ 7 ]CHILALI M, GAHINET P, APKARIAN P. Robust pole placement in LMI regions[J]. IEEE Transactions on Automatic Control, 2000, 2(12): 2257–2270.

[ 8 ]KAUTSKY J, NICHOLS N K, DOOREN P V. Robust pole assignment in linear state feedback[J]. International Journal of Control, 1985, 41(5): 1129–1155.

[ 9 ]LIU H, MA N, GU X. Maneuvering prediction of a VLCC model based on CFD simulation for PMM tests by using a circulating water channel[C]// ASME 2015 34th International Conference on Ocean, Offshore and Arctic Engineering. American Society of Mechanical Engineers, 2015.

[10]LIU H, MA N, SHAO Chuang. Numerical simulation of planar motion mechanism test and hydrodynamic derivatives of a ship in laterally restricted water [J]. Journal of Shanghai Jiao Tong University, 2016, 50(01): 115–122.

Control study of KVLCC2 maneuvering in restricted water based on LQR

SHAO Chuang1, MA Ning1,2, GU Xie-chong1,2
(1. School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiaotong University, Shanghai 200240, China; 2. State Key Laboratory of Ocean Engineering, Shanghai Jiaotong University, Shanghai 200240, China)

To investigate the hydrodynamic characters of a large ship in laterally restricted water, In restricted waters, ships will be affected by bank effect, sway force and yaw moment will change, which will have adverse effects on the control of ship navigation. In this paper, the LQR optimal control theory is applied to a large oil tanker KVLCC2 to achieve control design of ship maneuvering in restricted waters. In order to permit the design of a linear quadratic (LQ) controller, the maneuvering equations of motion are cast in a linear state-space framework, and CFD numerical simulation methods are used to compute the corresponding hydrodynamic coefficients. Based on the above, the paper calculates gain matrix K to minimize cost function and obtains the time-domain rudder angle satisfying optimal control law, which achieves optimal control of ship maneuvering in restricted waters with different widths. Simulation results demonstrate the effectiveness and superiority of the LQR controller compared with traditional pole placement method under the same control parameters. The results show that when the largest non-dimensional separation distance is d/L=1.2 the ship is merely disturbed and control action is very small. And the case of d/L=0.25 sees steady rudder angles beyond 6° and also forward speed reduction is 10% of the forward speed.

LQR；restricted water；optimal control

U661

A

1672 – 7649(2017)07 – 0073 – 05

10.3404/j.issn.1672 – 7649.2017.07.015

2016 – 09 – 21；

2016 – 11 – 04

邵闖(1989 – )，男，碩士研究生，主要從事限寬水域中船舶操縱及控制研究。