許煒東,張 軍,仲銘杰,楊一帆

(1.江蘇大學(xué) 電氣信息工程學(xué)院,江蘇 鎮(zhèn)江 212013; 2.江蘇航空職業(yè)技術(shù)學(xué)院,江蘇 鎮(zhèn)江 212013)

1 引言

高速多體船(包括兩體船和三體船)是當前新型船舶的重要發(fā)展趨勢之一，與單體船相比具有更小的興波阻力、橫向穩(wěn)定性好、裝載量大、抗打擊能力強，在軍事、航運領(lǐng)域有著廣泛的應(yīng)用前景[1]。然而，高速多體船細長的側(cè)體結(jié)構(gòu)使得縱向恢復(fù)力矩變小，造成船舶在高速航行時會產(chǎn)生較大的升沉和縱搖幅度，出現(xiàn)“埋首”運動[2]，嚴重破壞船體結(jié)構(gòu)。劇烈的升沉和縱搖運動產(chǎn)生較大的垂向加速度，造成乘客暈船、儀器失靈。因此，如何保證高速多體船的垂向運動穩(wěn)定，減少升沉和縱搖運動幅度，成為該領(lǐng)域研究的重點。目前，在多體船船體上安裝主動式T型翼和壓浪板兩類附體可實現(xiàn)閉環(huán)減搖控制，T型翼可以提高船體的附加質(zhì)量和阻尼，壓浪板用于改變船體縱傾[3]。

安裝T型翼和壓浪板的高速多體船是多輸入輸出系統(tǒng)，升沉和縱搖運動存在強耦合，并且航速較高導(dǎo)致水動力參數(shù)變化劇烈，模型具有快時變性、強不確定性等特點[4-5]。文獻[4]中設(shè)計解耦的比例微分(PD)控制律實現(xiàn)雙體船的垂向穩(wěn)定控制，但是需要精確辨識水動力學(xué)參數(shù)。文獻[5]中提出擴張觀測器和PD結(jié)合的解耦控制策略，提高了魯棒性。文獻[6]中設(shè)計定量反饋控制律，但是控制性能比較保守。文獻[7]中提出有限時間控制器和有限時間觀測器復(fù)合的控制策略，增強了多體船抑制海浪擾動的能力。上述研究主要從鎮(zhèn)定角度來設(shè)計解耦或者魯棒控制器，沒有考慮垂向控制性能的優(yōu)化，而垂向控制目的是實現(xiàn)升沉和縱搖變化幅度和加速度越小越好。

針對升沉和縱搖耦合的多體船垂向運動模型，預(yù)測控制不僅能有效處理多變量耦合，并且實現(xiàn)控制性能優(yōu)化。文獻[8]中提出多體船的解析預(yù)測控制，避免在線復(fù)雜的數(shù)值優(yōu)化。文獻[9]中對隨機海浪擾動力和力矩進行建模，采用卡爾曼濾波器估計狀態(tài)，提出顯式預(yù)測控制。然而，上述預(yù)測控制只考慮升沉和縱搖運動幅度的均方差最小，沒有考慮兩者的暫態(tài)響應(yīng)過程。實際上在海浪擾動激勵下升沉和縱搖的輸出不僅幅度大，波動變化率也大，造成垂向加速度很大，嚴重影響高速多體船的耐波性。針對上述問題，提出基于增量預(yù)測控制的垂向運動穩(wěn)定方法，其創(chuàng)新點為：

1) 分析隨機海浪擾動的頻域和時域特性，基于積分擾動模型在線估計海浪擾動力和力矩；為減少海浪平均擾動力和擾動力矩的影響，將高速多體船舶的垂向運動模型轉(zhuǎn)化為增量模型；

2) 為減少隨機海浪擾動下的升沉和縱搖幅度和波動變化率，將升沉和縱搖變化率和幅

度同時引入到優(yōu)化目標函數(shù)中，推導(dǎo)出解析的預(yù)測控制律。

具體章節(jié)安排如下：首先，建立升沉和縱搖描述的高速多體船垂向運動模型，分析隨機海浪的時域特性。其次，將高速多體船的垂向運動模型轉(zhuǎn)化為增量模型，基于滾動優(yōu)化策略推導(dǎo)預(yù)測控制律。最后，通過仿真驗證所提算法的有效性。

2 高速多體船的垂向運動模型

高速多體船的船體是相對于中間截面對稱的，T型翼和壓浪板分別在船艏底部和船艉。在航行中通過調(diào)節(jié)T型翼和壓浪板攻角，產(chǎn)生恢復(fù)力和力矩，抵消海浪引起的擾動力和力矩，抑制升沉和縱搖劇烈運動。當在康明斯方程中考慮流體記憶效應(yīng)時[10]，基于升沉和縱搖兩自由度描述的高速多體船垂向運動模型為

(1)

(2)

式中：ρ為海水密度；S1、S2分別為T型翼和壓浪板的有效面積；CL1和CL2分別為T型翼和壓浪板的升力系數(shù)；V表示流體相對于水翼的速度；α1、α2和lT-foil、lflap分別表示T型翼和壓浪板的攻角和力臂。

為后續(xù)設(shè)計方便，將式(1)與式(2)寫成下面的狀態(tài)空間形式

(3)

從式 (3) 可以看出，高速多體船模型存在升沉和縱搖運動的強耦合，以及存在強的海浪擾動。高速多體船垂向運動控制目的是在海浪擾動下，自動調(diào)節(jié)減搖附體的攻角使得升沉和縱搖運動的幅度和加速度越小越好。

3 海浪擾動力和力矩估計

高速多體船的升沉和縱搖幅度變化劇烈，主要原因是受到海浪擾動力和力矩作用，本節(jié)研究四級海浪產(chǎn)生擾動的時域和頻域特性，為后續(xù)的減搖控制設(shè)計提供基礎(chǔ)。本研究中采用ITTC單參數(shù)譜作為海浪仿真功率譜函數(shù)[1,11]，即

(4)

式中:h1/3為海浪有義波高；ω為海浪主導(dǎo)頻率。根據(jù)功率譜反演和隨機信號合成理論，得到作用于兩體船的擾動力和擾動力矩，如圖1、圖2所示。

圖1 多體船的海浪擾動力Fig.1 Disturbance wave force of the multihull

圖2 多提船的海浪擾動力矩Fig.2 Disturbance wave moment of the multihull

可以看出多體船的擾動力和力矩是一種平穩(wěn)隨機過程。通過統(tǒng)計數(shù)據(jù)知道海浪擾動力和擾動力矩均值分別為0.429 7×103N，2.028 2×103N·m，均不為零，并且海浪變化率Δd(k)≠0，這會影響下一節(jié)的預(yù)測控制模型準確性，需要在線實時估計d(k)，在文獻[10]中假設(shè)Δd(k)=0是有局限性的。

海浪擾動力和力矩的估計常采用白噪化方法，將白噪聲作為系統(tǒng)新的輸入，對系統(tǒng)狀態(tài)進行擴張，采用卡爾曼濾波器估計狀態(tài)[12]。但是，該方法需要辨識成型濾波器的參數(shù)，并且卡爾曼濾波維數(shù)較大，計算復(fù)雜。這里將海浪擾動力和力矩作為擴張狀態(tài)，基于積分擾動模型理論近似估計海浪擾動力和力矩[13]。將高速多體船模型方程式(3) 進行離散化可得

(5)

d(k+1)=d(k)+ζk(k)

(6)

其中：ζk(k)為零均值白噪聲；噪聲方差為Qr。將海浪擾動作為狀態(tài)的一部分，擴張后的狀態(tài)為

擴展的隨機線性模型如下

xe(k+1)=Aexe(k)+Beu(k)+Bewζ(k)

(7)

圖3 估計的海浪擾動力Fig.3 Estimated disturbance wave force

圖4 估計的海浪擾動力矩Fig.4 Estimated disturbance wave moment

4 減少變化率的增量預(yù)測控制

高速多體船的垂向運動模型方程式 (3) 是多變量，并且垂向運動控制目的是實現(xiàn)升沉和縱搖輸出和加速度越小越好，本節(jié)提出具有優(yōu)化功能的垂向穩(wěn)定預(yù)測控制[13-14]。由上一節(jié)分析知道，海浪擾動力和擾動力矩均值都不為零，需要有積分作用來提高抑制海浪擾動能力。因此，本節(jié)采用增量預(yù)測模型設(shè)計預(yù)測控制律。

Δx(k+1)=AdΔx(k)+BdΔu(k)+BwdΔd(k)

(8)

由于當前k時刻的擾動d是已知的，但未來時刻k+1，k+2…卻是未知并且不可測的，所以假設(shè)在k時刻之后的可測擾動近似不變，即Δd(k+i)=0,i=1,2,…,Np。

對于系統(tǒng)的輸出y(k)，有

y(k+1)=y(k)+CdAdΔx(k)+

CdBdΔu(k)+CdBwdΔd(k)

(9)

將y(k)擴展為一個新的狀態(tài)變量，擴展后的狀態(tài)空間模型為

(10)

對k時刻，定義Np步預(yù)測輸出和Nc步控制輸入

YNp(k+1)=[y(k+1|k),y(k+2|k),…,y(k+Np|k))]T

ΔU(k)=[Δu(k),Δu(k+1),…,Δu(k+Nc-1)]T

則基于第k步系統(tǒng)狀態(tài)，未來k+Np步的預(yù)測輸出為

YNp(k+1)=Fxu(k)+ΘΔU(k)+GΔd(k)

(11)

式中，

高速多體船的垂向運動穩(wěn)定控制不僅局限于減少縱搖和升沉運動幅度，還需要減少縱搖和升沉的變化率，其原因是過大的縱搖和升沉變化率會引起多體船垂向加速度過大，進而使得暈船率過大。這里將預(yù)測時域內(nèi)的縱搖和升沉變化率也引入優(yōu)化目標函數(shù)，提出了下面的預(yù)測控制目標函數(shù)

J=(Yref(k+1)-YNp(k+1))TQ(Yref(k+1)-

YNp(k+1))+ΔU(k)TRΔU(k)+

(YNp(k+1)-YNp(k))TK(YNp(k+1)-YNp(k))

(12)

ΔU(k)=(ΘTQΘ+R+ΦTKΦ)-1[ΘTQ(Yref(k+1)-

Fxu(k)-GΔd(k))-ΦTK(Γxu(k)+ΨΔd(k))]

(13)

式中，

為了提高預(yù)測控制的魯棒性，只取ΔU(k)中的第一個元素作為實際的控制增量，在下一采樣時刻的控制增量根據(jù)下一采樣時刻的系統(tǒng)狀態(tài)重新計算。令矩陣H=[I2×20 … 0]，則當前時刻的MPC控制量

Δu(k)=H(ΘTQΘ+R+ΦTKΦ)-1[ΘTQ(Yref(k+1)-

Fxu(k)-GΔd(k))-ΦTK(Γxu(k)+ΨΔd(k))]

(14)

為了驗證上述閉環(huán)增量MPC的穩(wěn)定性，將增量MPC控制器方程式 (14) 代入系統(tǒng)方程式 (10) 可得

xu(k+1)=BuH(ΘTQΘ+R+ΦTKΦ)-1·

ΘTQYref(k+1)+(Au-BuH(ΘTQΘ+R+ΦTKΦ)-1·

(ΘTQF+ΦTKΓ))xu(k)+

(Bwu-BuH(ΘTQΘ+R+ΦTKΦ)-1·

(ΘTQG+ΦTKΨ))Δd(k)

(15)

閉環(huán)系統(tǒng)的穩(wěn)定性通過矩陣Au-BuH(ΘTQΘ+R+ΦTKΦ)-1(ΘTQF+ΦTKΓ)的特征值來判斷。當加權(quán)矩陣為Q=diag(6.2×103，4.8×104)，R=diag(3.4×10-3，2.4×10-5)，K=diag(1.2×106，10×107)，計算可得6個特征值分別是0.963 7+0.138 4i，0.963 7-0.138 4i，0.976 7+0.067 5i，0.976 7-0.067 5i，0.995 2+0.000 0i，0.999 9+0.000 0i均位于單位圓內(nèi)，因此離散系統(tǒng)是閉環(huán)穩(wěn)定。

5 仿真分析

在仿真中高速多體船以兩體船為例，取采樣周期T=0.05 s，動力學(xué)模型參數(shù)為[12]：

為驗證增量預(yù)測控制減搖的有效性，采用四級海浪進行仿真，仿真條件同文獻[9]，進行下面仿真：

1) 分別對多體船進行無控制器和增量預(yù)測控制的仿真，升沉位移和縱搖角圖5和圖6所示。無控時升沉位移均值為0.240 6 m，方差為0.089 1 m2，縱搖角均值為 0.736 3 deg，方差為0.014 6 deg2。在加入增量預(yù)測控制后，升沉位移均值為0.097 1 m，方差為0.015 4 m2，縱搖角均值為0.382 0 deg，方差為0.004 2 deg2。通過對比可以看出，在增量預(yù)測控制作用下多體船升沉位移減少了59%，縱搖角減少了47%，減搖控制效果明顯。

圖5 升沉運動Fig.5 Comparison of heave motion

圖6 縱搖運動Fig.6 Comparison of pitch angle motion

2) 為進一步驗證所提增量預(yù)測控制性能，與傳統(tǒng)預(yù)測控制減搖進行比較[10]，結(jié)果如圖7和圖8所示，可以看出所提的預(yù)測控制使得升沉位移幅度進一步減少16%，縱搖運動幅度減少約11%，波動變化率明顯變緩，增強了多體船抑制海浪擾動的能力。

圖7 升沉運動Fig.7 Comparison of heave motion

圖8 縱搖運動Fig.8 Comparison of pitch angle motion

6 結(jié)論

提出基于增量預(yù)測控制的高速多體船垂向運動穩(wěn)定方法，抑制過大的升沉和縱搖運動幅度。仿真結(jié)果表明：

1) 采用積分擾動模型在線估計海浪擾動力和力矩，將估計值引入到預(yù)測模型中，增強多體船抑制海浪擾動的能力。

2) 基于增量模型設(shè)計預(yù)測控制減搖方法，使得升沉位移減少59%，縱搖角減少47%。