楊文輝 黃道平
(華南理工大學(xué)自動化科學(xué)與工程學(xué)院 廣州 510640)
圖1 船舶轉(zhuǎn)向操作
Job Van Amerongen在對船舶進(jìn)行研究的時候,提出了船舶轉(zhuǎn)向的理想操作[1],如圖1所示。船舶轉(zhuǎn)向操作分成三個階段,0~t1:轉(zhuǎn)向啟動;t1~t2:轉(zhuǎn)向保持;t2~t3:轉(zhuǎn)向結(jié)束。不同使用者對船舶轉(zhuǎn)向的共同要求:1)船舶的轉(zhuǎn)向速度可以調(diào)節(jié);2)船舶達(dá)到預(yù)定航向的時候,不會出現(xiàn)超調(diào)現(xiàn)象。
典型的船舶轉(zhuǎn)向控制方法包括模型參考自適應(yīng)控制,前饋加LQ控制器、模型跟蹤器,以及基于Lyapunov穩(wěn)定性理論的設(shè)計,在文獻(xiàn)[2]提出了一種由低通濾波器、前饋控制器、反饋控制器組成的轉(zhuǎn)向控制系統(tǒng),并且按轉(zhuǎn)向角的大小,分成三個不同的控制區(qū)域,分別進(jìn)行控制。文獻(xiàn)[3]擴(kuò)展了文獻(xiàn)[2]中對轉(zhuǎn)向角分區(qū)控制的方法,將轉(zhuǎn)向角度分成9個區(qū)域。這些控制器設(shè)計都有自己的優(yōu)點(diǎn)。本文提出了一種航向成形算法,給定航向 ψd,經(jīng)成形算法后變成理想轉(zhuǎn)向參考信號ψr,并作為控制系統(tǒng)的參考航向與經(jīng)低通濾波器之后的測量航向值之差Δ ψ作為控制器的輸入信號??刂葡到y(tǒng)的設(shè)計,如圖2所示。其中低通濾波器的作用是濾除高頻波浪干擾的信號的影響。
圖2 船舶轉(zhuǎn)向控制系統(tǒng)
式(2)廣泛應(yīng)用于船舶自動舵的控制器設(shè)計中。求解該模型需要知道的8個船舶參數(shù):航速V(m/s)、船舶兩柱間長L(m)、船寬B(m)、滿載吃水T(m)、船舶方形系數(shù) Cb、排水量(m3)、船舶重心距中心距離 xC(m)、舵葉面積 Aδ(m2),具體計算可參考文獻(xiàn)[4]。
在研究船舶的操縱運(yùn)動時候一般只是考慮平面運(yùn)動。日本學(xué)者野本謙作通過大量的實(shí)船試驗(yàn),對船舶的航向控制運(yùn)動方程進(jìn)行了簡化,得到低頻范圍內(nèi)的舵角δ,到航向 ψ的一階模型[4~5]:
通過分析,可以發(fā)現(xiàn),在船舶轉(zhuǎn)向中,存在兩種情況,情況1:如圖3所示,可以將轉(zhuǎn)向過程分成三個階段,分別是一個加速、勻速、減速的過程。情況1多出現(xiàn)在轉(zhuǎn)向角度較大的情況下。情況2:如圖4所示,轉(zhuǎn)向過程沒有勻速階段,這種情況多出現(xiàn)在轉(zhuǎn)向角度小,且轉(zhuǎn)向速度大時,例如在一些直航道上面的緊急調(diào)整航向。但并不否定轉(zhuǎn)向角度小就不能出現(xiàn)情況1所述的三個轉(zhuǎn)向過程,此時只要選擇一個較小的轉(zhuǎn)向速度仍然會出現(xiàn)情況1所述的過程。下面分別介紹兩種情況下航向成形算法。
為了簡化問題,假設(shè)情況1中出現(xiàn)的三個階段分別是:0~t1:勻加速;t1~t2:勻速;t2~t3:勻減速,且加速和減速階段的加速度a(°/s2)大小的相同,這樣加速和減速過程就是對稱的。當(dāng)轉(zhuǎn)向角度ψd,轉(zhuǎn)向速度rd,由操作人員設(shè)定之后,就可以通過求解方程(3)~(5),得到時間 t1、t2、t3的值。
關(guān)于加速度a的取值,從圖5中可以發(fā)現(xiàn)船舶轉(zhuǎn)向的加速過程是一個加速度逐漸減小的變加速過程,如果直接根據(jù)這樣的時間-速度關(guān)系來設(shè)計,則加速度a的值就會很小,考慮到船舶的大慣性特點(diǎn),這里建議將a取在0.1附近,以加速船舶的轉(zhuǎn)向啟動。
關(guān)于轉(zhuǎn)向速度r,在船舶轉(zhuǎn)向中,其轉(zhuǎn)向速率受到舵角的影響,就長時間來看,舵角與穩(wěn)態(tài)的轉(zhuǎn)向速度存在下面的關(guān)系:其中:ψ表示航向角,r表示船舶轉(zhuǎn)向速率,K為野本模型中的增益系數(shù),δ舵角。不同舵角下船舶轉(zhuǎn)向速度,如圖5所示。同時因?yàn)榇嬖诙娼窍拗?船舶的最大轉(zhuǎn)向速度將由船舶允許的最大舵角決定,通常船舶最大舵角介于 25°~35°之間,當(dāng)船舶的野本模型參數(shù)得到之后,就可以確定船舶的最大轉(zhuǎn)向速率:
而且在實(shí)際船舶轉(zhuǎn)向操作中發(fā)現(xiàn),當(dāng)轉(zhuǎn)向角度大時,操作者希望使用較大的轉(zhuǎn)向速度,以縮短轉(zhuǎn)向的時間,此時舵機(jī)可能滿舵,以達(dá)到較大的轉(zhuǎn)向速度。當(dāng)轉(zhuǎn)向角度較小的時候,操作者希望使用較小的轉(zhuǎn)向速度,以精確控制船舶達(dá)到設(shè)定航向。
從舵機(jī)使用的角度考慮,為了減少舵機(jī)的磨損,希望舵機(jī)能夠平穩(wěn)的操作,避免出現(xiàn)劇烈的操舵行為。從這個角度看,也是希望轉(zhuǎn)向角度較小時,使用小的轉(zhuǎn)向速度,以減少舵機(jī)的劇烈操作。基于這樣的考慮,在非緊急情況下,可以下面的公式設(shè)定轉(zhuǎn)向速度:其中:α速度系數(shù),通過表1得到。當(dāng)出現(xiàn)緊急情況的時候,可以直接設(shè)定轉(zhuǎn)向速度為rmax。
表1 轉(zhuǎn)向速度系數(shù)表
情況2中的船舶轉(zhuǎn)向過程只存在加速和減速的過程,假設(shè):0~t1:勻加速;t1~t2:勻減速,且加速度a(°/s2)大小相同,當(dāng)轉(zhuǎn)向角度 ψd給定之后,就可以由式(10)、(11)計算得到時間t1、t2的值。
加速度a的取值,同情況1。此時的航向成形算法:
基于閉環(huán)增益成形算法的魯棒PID控制器[6],對于二階對象:
其中:T1是根據(jù)系統(tǒng)的有效截止頻率1/T1選擇的時間常數(shù)。
船舶難免會受到高頻波浪干擾的作用,而這些干擾,會引起高頻的操舵行為,所以在設(shè)計控制系統(tǒng)的時候,應(yīng)該對測量的航向進(jìn)行濾波,以排除高頻波浪干擾的作用,減少舵機(jī)的磨損。通常波浪的頻率介于0.05~0.2Hz之間[7],為此可以考慮設(shè)計一個低通濾波器。
以5446TEU集裝箱船COSCO Shanghai號[8]為仿真模型,船舶額定船速 U=24.5knots=12.6028m/s,其它參數(shù)如表2所示。
表2 5446TEU集裝箱船Shanghai號
計算得到:K=0.0793,T=41.35(s),取Tf=10,T1=50,對船舶進(jìn)行大角度轉(zhuǎn)向和小角度轉(zhuǎn)向仿真研究,結(jié)果如圖6、7所示。從仿真結(jié)果可以看出,經(jīng)過航向成形算法改進(jìn)之后,船舶的轉(zhuǎn)向平穩(wěn),沒有出現(xiàn)超調(diào),舵機(jī)的使用更加合理,轉(zhuǎn)向角度大舵角大,轉(zhuǎn)向角度小舵角小。
本文基于理想船舶轉(zhuǎn)向操作,設(shè)計的航向成形算法,其原理簡單,易于實(shí)現(xiàn),船舶轉(zhuǎn)向過程更加的平穩(wěn),沒有出現(xiàn)超調(diào)的現(xiàn)象,而且轉(zhuǎn)向速度可以根據(jù)轉(zhuǎn)向角度的大小或者用戶要求進(jìn)行選擇,符合理想船舶轉(zhuǎn)向要求。
[1]Van Amerongen,J.,Van Nauta Lemke,H.R.Criteria for optimum steering of ships[C]//Proceedings of Symposium on ship steering automatic control,Genoa,Italy,1980
[2]賈欣樂,蔣丹東,張顯庫.船舶轉(zhuǎn)向控制器研究[J].大連海事大學(xué)學(xué)報,1998,24(1):23~28
[3]張顯庫,趙翔宇.船舶轉(zhuǎn)向的魯棒控制及其優(yōu)化設(shè)計[J].哈爾濱工程大學(xué)學(xué)報,2006,27(3):319~322
[4]張顯庫,賈欣樂.船舶運(yùn)動控制[M].北京:國防工業(yè)出版社,2006
[5]李殿璞.船舶運(yùn)動與建模[M].哈爾濱:哈爾濱工程大學(xué)出版社,1999
[6]張顯庫,賈欣樂.求PID參數(shù)新方法[J].工程與電子技術(shù),200,22(8):4~5
[7]Fossen,T.I.Guidance and Control of Ocean Vehicles[M].Chichester:John Wiley&Sons,1999
[8]徐黎黎.船舶運(yùn)動智能控制及其虛擬仿真的研究[D].大連:大連海事大學(xué)碩士學(xué)位論文,2004