基于Hermite插值的管片拼裝機(jī)軌跡規(guī)劃

呂志斌，秦東晨，朱強(qiáng)，陳江義

(鄭州大學(xué)機(jī)械與動力工程學(xué)院，河南鄭州 450001)

0 前言

近幾年，隨著國內(nèi)區(qū)域經(jīng)濟(jì)一體化的不斷發(fā)展，軌道交通建設(shè)越來越熱火朝天。同時，各大城市不斷加大對地下空間的開發(fā)利用，使對周邊環(huán)境影響小、施工速度快、工程質(zhì)量優(yōu)良、施工安全環(huán)保、適用范圍廣的盾構(gòu)施工法越來越受到歡迎[1]，為盾構(gòu)機(jī)打開了一個空前繁榮的市場。近年來,盾構(gòu)機(jī)的國產(chǎn)化取得了長足的進(jìn)展,國內(nèi)多家單位都能完成盾構(gòu)機(jī)的獨(dú)立設(shè)計及制造[2]。管片拼裝機(jī)是盾構(gòu)機(jī)的重要組成部分，在管片拼裝過程中，盾構(gòu)停止掘進(jìn)，主機(jī)處于靜止?fàn)顟B(tài)[3]。在盾構(gòu)向前掘進(jìn)一環(huán)后，按預(yù)定要求將混凝土管片拼裝成環(huán)，形成襯砌，以此來支護(hù)剛開挖的隧道表面[4]。其拼裝質(zhì)量、精度及效率直接影響隧道施工安全和進(jìn)程[5]。

目前，KOSUGE等[6]進(jìn)行了面向任務(wù)的盾構(gòu)管片拼裝力控系統(tǒng)的研究；程永龍等[7]提出了微型盾構(gòu)半自動拼裝技術(shù)；劉飛香[8]研究了智能化抓取以及智能化拼裝方法；王明斗等[9]對矩形盾構(gòu)的串聯(lián)式管片拼裝機(jī)進(jìn)行了軌跡優(yōu)化。另外隧道股份上海隧道工程有限公司研發(fā)了中國首臺推拼同步盾構(gòu)機(jī)[10]；李杰等人[11]對管片拼裝機(jī)提升系統(tǒng)同步性能進(jìn)行了研究。綜合來看，這些學(xué)者在管片拼裝機(jī)的液壓仿真以及微調(diào)機(jī)構(gòu)等方面做出了大量工作，但是對其進(jìn)行軌跡優(yōu)化的研究較少。王明斗等采用的是B樣條曲線，它在兩點(diǎn)間最高只能構(gòu)造二階函數(shù)。文中的研究對象是成都11號線盾構(gòu)機(jī)的管片拼裝機(jī)，為雙提升臂拼裝機(jī)。結(jié)合Hermite插值算法對其進(jìn)行時間最優(yōu)的軌跡優(yōu)化，提出一種新的拼裝軌跡規(guī)劃方法，可以在條件較少的情況下構(gòu)造高階函數(shù)，使軌跡平滑，從而在降低啟動沖擊的同時，提升管片拼裝機(jī)的拼裝效率。

1 雙提升臂拼裝機(jī)

雙提升臂拼裝機(jī)的粗調(diào)機(jī)構(gòu)是由平移系統(tǒng)、回轉(zhuǎn)系統(tǒng)和提升系統(tǒng)三部分組成。回轉(zhuǎn)系統(tǒng)搭載在平移系統(tǒng)之上，隨平移系統(tǒng)一起運(yùn)動；提升系統(tǒng)搭載在回轉(zhuǎn)系統(tǒng)之上，從而實(shí)現(xiàn)平移和回轉(zhuǎn)。其中平移系統(tǒng)和提升系統(tǒng)由對稱的雙液壓缸來驅(qū)動；回轉(zhuǎn)系統(tǒng)由對稱的雙液壓馬達(dá)驅(qū)動，在保證驅(qū)動力的同時，確保其精確度。管片拼裝機(jī)的三維模型如圖1所示。管片拼裝時，平移系統(tǒng)控制水平位移，回轉(zhuǎn)系統(tǒng)控制回轉(zhuǎn)角度，提升系統(tǒng)控制管片高度。通過粗調(diào)機(jī)構(gòu)到達(dá)指定位置附近，再由微調(diào)機(jī)構(gòu)調(diào)整管片位姿偏差，完成管片拼裝。

圖1 管片拼裝機(jī)三維模型Fig.1 Three dimensional model of segment assembler

由于對接機(jī)器視覺，其輸出參數(shù)為直角坐標(biāo)系下的三維向量(X,Y,Z)，而管片拼裝機(jī)的輸入?yún)?shù)為(l,θ,d)。其中：l為平移系統(tǒng)的位移，θ為回轉(zhuǎn)系統(tǒng)的回轉(zhuǎn)角度，d為提升系統(tǒng)的位移。因此需要求管片拼裝機(jī)運(yùn)動學(xué)方程的逆解。根據(jù)管片拼裝機(jī)的結(jié)果易得其運(yùn)動學(xué)方程為

(1)

式中：r1為機(jī)械結(jié)構(gòu)在徑向上的偏移距離；r2為機(jī)械結(jié)構(gòu)在軸向上的偏移距離。

因此其運(yùn)動學(xué)方程的逆解為

(2)

2 軌跡規(guī)劃方法

由于管片質(zhì)量達(dá)到了4 t，因此其空載和負(fù)載的受力有很大不同，空載時可以提高管片拼裝機(jī)運(yùn)行速度，負(fù)載時降低其運(yùn)行速度，故將管片拼裝的動作分為抓取、遞送、安裝以及回原點(diǎn)4個部分，從而充分發(fā)揮各階段的特點(diǎn)，優(yōu)化各階段的軌跡，并且消除累計誤差對管片拼裝機(jī)自動拼裝的影響。

2.1 軌跡規(guī)劃方案

在始末速度優(yōu)化中常用的方法為擺線輪廓速度法[12]，但是通過計算發(fā)現(xiàn)采用Hermite插值法將進(jìn)一步降低沖擊力。以啟動時間1 s、最大速度1 m/s為例，對比兩者的速度和加速度曲線如圖2所示。

圖2 兩者速度(a)和加速度(b)對比

通過對比擺線輪廓速度法，發(fā)現(xiàn)在相同時間內(nèi)提升到相同的速度三次Hermite插值法的加速度較小，減小了4.5%左右。也就是說在負(fù)載的情況下，啟動沖擊力減小了2 kN左右。

因此本文作者選用Hermite插值法對管片拼裝機(jī)進(jìn)行軌跡規(guī)劃，其流程如圖3所示。

圖3 軌跡規(guī)劃流程Fig.3 Trajectory planning flow chart

為了提高管片拼裝效率，本文作者在保證加速度安全的前提下，采用三次Hermite插值算法將其速度平穩(wěn)而快速地提升。由于每個動作都是由管片拼裝機(jī)的平移、回轉(zhuǎn)以及提升系統(tǒng)共同作業(yè)完成的，因此可以將其分解為3條速度-時間曲線。首先基于最大加速度的限制比較出各自的最快完成時間，選取用時最長的作為時間基準(zhǔn)，其他兩個動作基于這個時間基準(zhǔn)和已知的位移量。為了保證拼裝機(jī)高效率和平穩(wěn)的運(yùn)行，其運(yùn)動軌跡曲線一般需要具有C2連續(xù)性[13]，因此采用五次Hermite插值算法規(guī)劃出路徑，從而使其在達(dá)到最優(yōu)速度的前提下進(jìn)一步降低其加速度。

2.2 加速度計算

將已知條件列表1如下。

表1 已知參數(shù)Tab.1 Known parameters

對于平移系統(tǒng)，只需考慮推力。

由F=ma可得，負(fù)載工況下

(3)

空載工況下

(4)

但是由于管片拼裝機(jī)的額定功率為55 kW，經(jīng)過計算在抓取階段實(shí)際功率超出額定功率。經(jīng)調(diào)整，空載工況下的加速度也改為4.17 m/s2。

對于回轉(zhuǎn)系統(tǒng)，為了保證所求的角加速度在整個運(yùn)動過程中都能達(dá)到，必須考慮其極限工況。即當(dāng)回轉(zhuǎn)角度為±90°時，重力扭矩最大，總扭矩最小。由T=Jα可得，負(fù)載工況下

(5)

空載工況下

(6)

其中：T1為管片重力產(chǎn)生的扭矩；T2為提升系統(tǒng)自重產(chǎn)生的扭矩；J1為管片繞回轉(zhuǎn)中心的轉(zhuǎn)動慣量；J2為提升系統(tǒng)繞回轉(zhuǎn)中心的轉(zhuǎn)動慣量；J3為回轉(zhuǎn)系統(tǒng)繞回轉(zhuǎn)中心的轉(zhuǎn)動慣量。

經(jīng)過計算在回原點(diǎn)階段實(shí)際功率超出額定功率。經(jīng)調(diào)整，空載工況下的加速度改為2.8 rad/s2。

對于提升系統(tǒng)，只需考慮提升力與重力。

由F=ma可得，負(fù)載工況下

12.84 m/s2

(7)

空載工況下

(8)

由于提升系統(tǒng)用時不長，通常不作為基準(zhǔn)，不會滿負(fù)荷運(yùn)轉(zhuǎn)，因此加速度不做調(diào)整。

2.3 軌跡規(guī)劃方法

為了降低起停沖擊，提高管片拼裝機(jī)的運(yùn)行平穩(wěn)性，將起始位置和末尾位置的速度以及加速度設(shè)為零。

在平移過程中，假設(shè)其運(yùn)行時間為t1，最高速度為v1，先用三次Hermite插值法表達(dá)出速度與時間的表達(dá)式

(9)

同理，回轉(zhuǎn)和提升的表達(dá)式分別為

(10)

(11)

然后根據(jù)算出的最大加速度和輸入的位移量，分別對多項(xiàng)式進(jìn)行求導(dǎo)和積分，解出時間t和速度v。

通過對比三者的時間，選取時間最長的為基準(zhǔn)動作，可得出該動作速度與時間的曲線。對多項(xiàng)式進(jìn)行求導(dǎo)和積分，可得到加速度與時間的曲線以及位移與時間的曲線。

對于其他兩個動作，根據(jù)已知的位移量和時間基準(zhǔn)t，采用五次Hermite插值算法分別對位移進(jìn)行分析，可計算出其余兩個動作位移與時間的曲線，求導(dǎo)可得到速度與時間的曲線以及加速度與時間的曲線。

3 計算結(jié)果

通過管片拼裝機(jī)的設(shè)計圖紙，確定了各個動作的平移量、回轉(zhuǎn)角度和提升量。其中液壓缸伸出為正，縮回為負(fù)；順時針回轉(zhuǎn)為正，逆時針回轉(zhuǎn)為負(fù)。

以K管片為例[14]，其中遞送時需注意，當(dāng)提升系統(tǒng)快達(dá)到管片厚度時，平移系統(tǒng)和回轉(zhuǎn)系統(tǒng)均要到達(dá)指定位置，以防碰到已裝好的管片。粗調(diào)參數(shù)列于表2中。

表2 K管片粗調(diào)參數(shù)Tab.2 K segment rough adjustment parameters

運(yùn)用MATLAB算法得出管片拼裝機(jī)的Hermite插值函數(shù)，將其導(dǎo)入ADAMS進(jìn)行運(yùn)動學(xué)仿真，繪制出K管片整體動作管片的位移、速度以及加速度的曲線如圖4所示。

圖4 運(yùn)動學(xué)曲線

從圖4中可以看出：在K管片拼裝過程中，管片拼裝機(jī)運(yùn)行平穩(wěn)。拼裝分速度最大出現(xiàn)在回原點(diǎn)階段的Y軸方向，約為2.7 m/s；合速度最大出現(xiàn)在回原點(diǎn)階段，約為2.7 m/s。加速度曲線連續(xù)且無突變，能夠使運(yùn)動平穩(wěn)，從而顯著降低了系統(tǒng)沖擊。其最大分加速度出現(xiàn)在回原點(diǎn)階段X軸方向，約為5 m/s2。最大合加速度出現(xiàn)在回原點(diǎn)階段，約為5 m/s2。

通過已知數(shù)據(jù)計算出各動作的大致運(yùn)行時間，列于表3中。

表3 運(yùn)動時間 單位：s

由表3可以看到，目前常用的擺線輪廓法作業(yè)各個系統(tǒng)在運(yùn)行期間的總用時大致為40.639 s；即使采用3個系統(tǒng)并行運(yùn)行，其總用時大致也要26.097 s；不限速度后的總用時大致為10.377 s。而采用文中的方法其總用時大致為9.519 s。相比于現(xiàn)在遙控拼裝常用的擺線輪廓法作業(yè)，總用時降低了30.82 s，約75.8%；就算采用并行運(yùn)行，總用時也降低了16.278 s，約62.4%。相比于不限速度的擺線輪廓法，總用時降低了0.858 s，約8.3%。因此文中采用的五次Hermite插值法拼裝效率最高。

由于大多學(xué)者沒有對此進(jìn)行深入研究，因此沒有過多的數(shù)據(jù)對比。且由于管片拼裝機(jī)的微調(diào)整和后續(xù)環(huán)節(jié)主觀性太大，無法用具體數(shù)據(jù)進(jìn)行精確分析，文中做了相應(yīng)的簡化，僅做了管片拼裝機(jī)粗調(diào)整的軌跡規(guī)劃，但也大大提升了管片拼裝機(jī)的拼裝效率。

4 結(jié)論

(1)建立了管片拼裝機(jī)的運(yùn)動學(xué)模型，并基于Hermite插值算法提出了一種應(yīng)用于管片拼裝機(jī)的軌跡規(guī)劃方法，其加速度曲線連續(xù)，運(yùn)行平穩(wěn)，能大大提高管片拼裝機(jī)的拼裝效率并有效降低沖擊。

(2)編輯了基于Hermite算法的MATLAB程序，輸入平移量、旋轉(zhuǎn)量、提升量以及動作序號可得到平移、回轉(zhuǎn)以及提升系統(tǒng)的運(yùn)動函數(shù)。將其導(dǎo)入ADAMS仿真軟件，得出管片拼裝時的運(yùn)動學(xué)曲線。不僅通過仿真得出了直觀的運(yùn)動學(xué)曲線；還可以為以后接入機(jī)器視覺輸入、實(shí)現(xiàn)自動化施工[15]和智能控制打下基礎(chǔ)。

(3)采用Hermite插值算法規(guī)劃了一條平穩(wěn)的拼裝路徑，相比于擺線輪廓法減小了4.5%的啟動沖擊，更是將拼裝效率提升了約75.8%。就算將來引入智能化拼裝，采用擺線輪廓法并行運(yùn)行，拼裝效率也提升了約62.4%；即使打破速度限制，拼裝效率也提升了約8.3%。