林靜正,方勇純,盧 彪,郝運嵩,曹海昕
(南開大學(xué)人工智能學(xué)院機器人與信息自動化研究所,天津 300350)
近些年,伴隨著人類社會的高速繁榮,加大資源開采,擴張活動領(lǐng)域是必然的發(fā)展趨勢.作為覆蓋地球71%表面積的重要戰(zhàn)略資源,海洋逐漸受到關(guān)注,海上作業(yè)活動也日益頻繁,比如貨物的吊運輸送、設(shè)備的投放回收、船只間的補給等.在這些作業(yè)任務(wù)中,船用起重機扮演著不可或缺的角色.因此,對其展開研究具有非常重要的實際意義和廣闊的應(yīng)用前景.
與陸上起重機類似,船用起重機是一種典型的欠驅(qū)動系統(tǒng).目前,對陸上起重機系統(tǒng)的研究已取得較大的進展:國內(nèi)外學(xué)者基于輸入整形[1]、滑模控制[2]、自適應(yīng)控制[3]、預(yù)測控制[4]等理論提出了多種控制器設(shè)計思路.然而,對船用起重機控制的研究還未成熟,仍面臨著許多難題.首先,船用起重機系統(tǒng)內(nèi)部存在摩擦力,未建模動態(tài),參數(shù)不確定等干擾,使得基于模型設(shè)計的控制器難以保證在實際應(yīng)用中達到預(yù)期的精度;更為重要的是,船用起重機需要在惡劣的海洋環(huán)境中進行作業(yè),不可避免地會受到海浪、大風(fēng)這類未知外部干擾的影響,使負載的定位精度下降,擺動加劇,給控制器的設(shè)計帶來了極大的挑戰(zhàn).
為了解決上述問題,國內(nèi)外學(xué)者陸續(xù)投入到對船用起重機的研究中來.在文獻[5]中,Fang等人利用歐拉-拉格朗日法對船用起重機進行了建模,并設(shè)計了非線性控制器使?fàn)顟B(tài)變量較好地跟蹤參考軌跡.在此基礎(chǔ)上,Lu等[6]通過狀態(tài)變量轉(zhuǎn)換,將船體運動引起的動力學(xué)方程與起重機自身動力學(xué)方程結(jié)合,簡化了控制器的設(shè)計.在文獻[7]中,Hong等討論了在海面較平穩(wěn)時海浪對船用起重機狀態(tài)的影響.一些學(xué)者基于滑??刂评碚揫8],模糊控制理論[9],神經(jīng)網(wǎng)絡(luò)算法[10]設(shè)計控制器來完成控制任務(wù).Küchler等[11]通過預(yù)測海浪來進行補償,同樣對海浪特性進行研究的還有文獻[12],然而,他們只關(guān)注了海浪的周期性質(zhì),并沒有充分考慮海浪的非周期部分,在真實環(huán)境下效果可能無法到達預(yù)期.此外,許多控制器都假設(shè)所建模型是完全準(zhǔn)確的,在控制器中包含大量的系統(tǒng)參數(shù).而在很多應(yīng)用場景下無法獲得精確的參數(shù)信息,因此會大大降低控制效果[13].
根據(jù)上述分析,本文提出了一種基于迭代學(xué)習(xí)和神經(jīng)網(wǎng)絡(luò)的控制方法,通過補償未知干擾,使起重臂和繩長到達期望位置,同時抑制負載擺動.本文首先將船用起重機動力學(xué)模型拆分為有驅(qū)和無驅(qū)部分以便后續(xù)控制器設(shè)計;接著通過周期估計器估計未知周期,并經(jīng)過迭代學(xué)習(xí)補償未知周期干擾;之后利用雙層神經(jīng)網(wǎng)絡(luò)對非周期干擾進行估計補償;在對未知干擾進行補償后,基于反饋線性化設(shè)計控制輸入.通過Lyapunov分析方法證明了期望平衡點是全局有界的.最后在所搭建的模擬船用起重機實驗平臺上進行了大量實驗,充分驗證了所設(shè)計控制器的有效性和魯棒性.本文的主要貢獻為
1) 提出了一種基于迭代學(xué)習(xí)和神經(jīng)網(wǎng)絡(luò)的控制方法.在有持續(xù)未知干擾的情況下,可以從理論上保證系統(tǒng)狀態(tài)的全局一致有界性.
2) 對船用起重機系統(tǒng)受到的未知干擾進行了充分的考慮,將其拆分成周期與非周期兩部分,利用迭代學(xué)習(xí)和神經(jīng)網(wǎng)絡(luò)分別進行補償.使系統(tǒng)對干擾的估計更加準(zhǔn)確完備,增強了控制器的魯棒性.
3) 基于Stewart平臺搭建了船用起重機實驗系統(tǒng),可以真實地模擬船用起重機工作環(huán)境.通過設(shè)計大量實驗,充分驗證了在惡劣環(huán)境下所提控制器的有效性.
本文其他部分安排如下:第2節(jié)對海浪干擾和船用起重機系統(tǒng)進行了建模分析,并描述了控制目標(biāo);第3節(jié)分別對周期和非周期干擾進行估計補償,并利用反饋線性化設(shè)計了一種基于迭代學(xué)習(xí)和神經(jīng)網(wǎng)絡(luò)的控制器;第4節(jié)通過Lyapunov方法進行了穩(wěn)定性分析;第5節(jié)展示了在模擬船吊平臺上的實驗結(jié)果;最后,在第6節(jié)對本文工作做了總結(jié)與展望.
本節(jié)分析了船用起重機受到的干擾,將未知干擾分為周期與非周期干擾.基于此給出了船用起重機在未知干擾下的動力學(xué)模型,并描述了具體的控制目標(biāo).
船用起重機系統(tǒng)受到的干擾主要由海浪引起.通常將這類海浪干擾dp視為由N個具有不同幅值、周期以及相位的正弦函數(shù)疊加的形式[11-12]:
其中:Ai,fi,φi分別代表第i個正弦函數(shù)的幅值、頻率和相位.定義T為的最小公倍數(shù),則式(1)中的干擾具有周期性:
本文將這部分干擾視為周期性的,且結(jié)構(gòu)和參數(shù)未知的一種干擾.需要說明的是,出于安全考慮,在海況極其惡劣的情況下,船只停止作業(yè).因此本文考慮的船用起重機工作時所受海浪干擾有界,作如下合理假設(shè).
假設(shè)1未知周期干擾具有已知上界Md,其周期T未知,但同樣具有已知上界TM[15-16].
當(dāng)海況平穩(wěn)時,海浪可看作是若干正弦波疊加的模型,但當(dāng)海面有風(fēng)浪時,干擾會體現(xiàn)出非周期性質(zhì)[17].并且,系統(tǒng)內(nèi)部存在各種摩擦力、參數(shù)不確定等未知干擾.因此,在對船用起重機所受干擾做了詳細考慮分析后,本文將未知干擾分為周期未知干擾與非周期未知干擾兩部分,分別進行估計補償.
圖1是船用起重機的示意圖,由船體、起重臂、吊繩、負載組成.Ts={Os,xs,ys,zs}表示動坐標(biāo)系,Te={Oe,xe,ye,ze}表示大地坐標(biāo)系.由海浪引起的船體俯仰運動以及升沉運動分別為α(t),h(t).基于歐拉-拉格朗日方程可以得到如下的動力學(xué)模型[14]:
圖1 船用起重機示意圖Fig.1 Schematic model of the offshore boom crane
其中:
為了便于隨后的分析,把式(3)拓展為以下形式:
本節(jié)中,首先設(shè)計了周期估計器對未知周期進行預(yù)測,并通過迭代學(xué)習(xí)控制補償未知周期干擾;接著設(shè)計了雙層神經(jīng)網(wǎng)絡(luò)以及權(quán)重更新律,來補償未知非周期干擾;最后基于反饋線性化設(shè)計控制器.
本小節(jié)在上述對未知干擾估計結(jié)果的基礎(chǔ)上進行控制器的設(shè)計.對式(18)中的r求導(dǎo)并代入式(10)中:
根據(jù)式(16)(19)(25)(30),設(shè)計神經(jīng)網(wǎng)絡(luò)與迭代學(xué)習(xí)控制相結(jié)合的反饋控制器:
本節(jié)將對閉環(huán)系統(tǒng)期望平衡點的穩(wěn)定性進行分析.
定理1對于受到海浪等未知干擾的欠驅(qū)動船用起重機系統(tǒng)(10),通過所設(shè)計的迭代學(xué)習(xí)更新律(16)(19)以及神經(jīng)網(wǎng)絡(luò)權(quán)重更新律(30),能夠?qū)ξ粗蓴_進行有效估計和補償.在此基礎(chǔ)上,提出的控制器(32)能夠保證在有持續(xù)未知干擾的情況下,使起重臂以及繩長跟蹤誤差漸近收斂,負載擺角全局有界:
為了證明定理1,設(shè)計如下李雅普諾夫候選函數(shù):
廣東省9家“三甲”醫(yī)院藥師對互聯(lián)網(wǎng)藥事服務(wù)認(rèn)知情況的調(diào)查分析 ……………………………………… 汪小惠等(2):159
其中:第1項積分項是為了解決迭代學(xué)習(xí)的學(xué)習(xí)誤差,第2項為誤差信號相關(guān)項,后兩項是為了解決神經(jīng)網(wǎng)絡(luò)權(quán)重的估計誤差所設(shè)計.對式(37)求導(dǎo),并代入學(xué)習(xí)律(19)-(21)(30)以及閉環(huán)系統(tǒng)方程(34),化簡可得
至此,根據(jù)式(52)(54),可以得出結(jié)論:在有持續(xù)未知干擾的情況下,起重臂q1(t)和繩長q2(t)漸近收斂,負載擺角q3(t)全局有界.對定理1的證明結(jié)束.
本節(jié)通過在模擬船吊實驗平臺上進行實驗,對所提出的方法進行了驗證.所搭建的實驗平臺如圖2所示.一個縮小比例的懸臂起重機安裝在一個Stewart六自由度平臺上,該平臺能夠提供所需要的俯仰、升沉運動來模擬船用起重機的工作環(huán)境.懸臂和吊繩分別由兩個帶有編碼器的伺服電機驅(qū)動,提供力/力矩的同時實時讀回幅角φ(t)與繩長L(t).另外,在懸臂頂端安裝有一個編碼器用來測量擺角θ(t).為了保證控制的實時性,控制周期設(shè)定為10 ms.
圖2 船用起重機實驗平臺Fig.2 Experimental testbed
實驗平臺的各項參數(shù)如下:
所設(shè)計的控制器的各項控制參數(shù)如下:
神經(jīng)網(wǎng)絡(luò)節(jié)點數(shù)為N=1,周期初值以及神經(jīng)網(wǎng)絡(luò)各權(quán)值初值設(shè)為
為了證明所提方法的優(yōu)勢,本節(jié)分為兩個子小節(jié).在第1小節(jié)的實驗中,將本文設(shè)計的方法與其他方法在相同條件下進行對比;在第2小節(jié)中,對所提方法加入各類干擾來驗證魯棒性.
本節(jié)將所設(shè)計的控制器與線性二次型控制器(linear quadratic regulator,LQR)和文獻[20]提出的自適應(yīng)增益滑??刂破?以下用SMC表示)進行對比.將式(3)進行適當(dāng)線性化后,設(shè)計Q矩陣為
R為矩陣,R=diag{1,1},利用MATLAB工具箱求解里卡蒂微分方程,得到控制器輸入
SMC控制器控制參數(shù)設(shè)計為k11=k21=20,k12=k22=5,γ1=γ2=3,λ1=λ2=0.1.
為了方便控制器之間的比較,定義如下性能指標(biāo):1)q3max:最大負載擺角的絕對值;2)q3var:負載擺角的方差;3)tq1:懸臂到達穩(wěn)定的時間;4)tq2:繩長到達穩(wěn)定的時間;5)xmean:x與xd之差的平方的平均數(shù);6)ymean:y與yd之差的平方的平均數(shù).
實驗結(jié)果如圖3-4和表1所示.從圖3中可以看到,在有持續(xù)未知擾動的情況下,本文方法具有更好的控制效果.具體而言,本文方法控制起重臂和繩長到達期望位置所需時間更短;而對比方法到達期望位置所需時間較長,并且在期望位置附近有較明顯的振蕩.同時,本文所提方法能夠更好地抑制負載擺角,而對比方法由于不能很好地補償船體引起的干擾,表現(xiàn)出負載在期望位置附近來回搖擺,各項性能指標(biāo)均劣于所提方法(具體指標(biāo)見表1).因此證明了本文方法對未知干擾進行估計和補償?shù)挠行?更適用于實際船用起重機作業(yè)場景.
圖3 本文設(shè)計方法,LQR方法,SMC方法[20]的實驗結(jié)果Fig.3 Experimental results of the proposed method,LQR method,and the SMC method[20]
表1 各類方法的性能指標(biāo)Table 1 Performance indices of different methods
圖4的前兩張圖表示伺服電機的輸出力/力矩.可以看到,本文方法的輸出變化更平緩,峰值更小,而對比方法抖振明顯,增加了安全隱患.圖4剩余部分為對周期的估計,對周期干擾以及非周期干擾的補償結(jié)果.在實驗中給定俯仰和升沉運動的周期分別為5 s,10 s.結(jié)果證明了本文方法能夠快速準(zhǔn)確地估計干擾周期.
圖4 本文設(shè)計算法,LQR算法,SMC算法[20]的實驗結(jié)果Fig.4 Experimental results of the proposed method,LQR method,and the SMC method[20]
為了驗證本文控制器的魯棒性,設(shè)計了兩組實驗.
實驗1在時間t=20 s時將俯仰運動與升沉運動的干擾改變?yōu)?/p>
結(jié)果如圖5-6所示.從圖5可以看到,當(dāng)干擾的周期以及幅值在控制過程中發(fā)生改變時,本文方法仍有很好的控制效果,負載的擺動始終控制在3°以內(nèi),且能夠?qū)崿F(xiàn)負載的定位.從圖6中可以看到,在t=20 s時干擾α(t)和h(t)的幅值和周期都發(fā)生了變化,周期估計器可以很好地跟蹤并且估計出干擾的周期,進而對干擾進行補償,圖中也很好地說明了干擾的周期與周期估計器得到的周期是一致的.證明了該方法對干擾的變化有較好的魯棒性,在實際應(yīng)用中具有重要意義.
圖5 干擾周期改變情況下的實驗結(jié)果Fig.5 Experimental results under the condition of changing disturbance period
圖6 干擾周期改變情況下控制器對周期干擾估計補償結(jié)果Fig.6 Controller’s estimation and compensation for period disturbance under the condition of changing disturbance period
實驗2將控制算法中包含的實驗平臺各項參數(shù)修改為m=10 kg,mp=1 kg,J=4 kg·m2;而實際參數(shù)為m=5 kg,mp=0.5 kg,J=0.4 kg·m2.
結(jié)果如圖7-8所示.圖7給出了所設(shè)計控制器在參數(shù)不準(zhǔn)確的情況下的控制效果.可以看到,即使存在系統(tǒng)參數(shù)不準(zhǔn)確問題,控制器依然能夠在短時間內(nèi)控制起重臂q1,繩長q2到達到期望位置,同時抑制負載的擺動并實現(xiàn)負載定位.圖8左側(cè)一列給出的是參數(shù)正確情況下控制器對非周期干擾的補償,而右側(cè)一列為參數(shù)不正確情況下的補償圖.通過左右兩列圖的對比可以發(fā)現(xiàn),當(dāng)系統(tǒng)參數(shù)設(shè)定不準(zhǔn)確時,所帶來的干擾可以通過神經(jīng)網(wǎng)絡(luò)來進行補償.說明了對于船用起重機系統(tǒng),考慮未知干擾的非周期部分是必要的,也證明了神經(jīng)網(wǎng)絡(luò)對非周期干擾的補償能力.在實際作業(yè)中,由于設(shè)備的尺寸及安裝的位置等原因,常常無法精確地獲得船用起重機所有的系統(tǒng)參數(shù),而所提方法不依賴于模型參數(shù)的準(zhǔn)確性,具有良好的魯棒性,應(yīng)用前景廣闊.
圖7 參數(shù)不確定情況下的實驗結(jié)果Fig.7 Experimental results under parameters uncertainty
圖8 參數(shù)不確定情況下控制器對非周期干擾估計補償對比Fig.8 Comparison of controller’s estimation and compensation for non-periodic disturbance under parameters uncertainty
本文針對未知干擾下船用起重機的控制問題提出了一種基于迭代學(xué)習(xí)和神經(jīng)網(wǎng)絡(luò)的方法.該方法不僅能夠?qū)ξ粗芷诟蓴_進行周期估計,基于迭代學(xué)習(xí)做出補償,還能通過雙層神經(jīng)網(wǎng)絡(luò)對非周期性干擾進行估計補償,提升了系統(tǒng)的抗干擾能力.利用Lyapunov方法進行了嚴(yán)格的分析,從理論上證明了期望平衡點的全局有界性.最后,在基于Stewart的船用起重機實驗平臺上進行了大量實驗,充分證明了所提方法在有未知干擾、摩擦、參數(shù)不確定等情況下的有效性和魯棒性.未來會將該控制問題提升到三維層面,在考慮俯仰、升沉、扭轉(zhuǎn)運動的情況下完善控制器設(shè)計.
致謝
感謝廣東省機器人與智能系統(tǒng)重點實驗室,中國科學(xué)院深圳先進技術(shù)研究院,對本工作的支持!