蘇召斌，林 森

（中交天津港航勘察設(shè)計研究院有限公司，天津 300457）

絞刀運動姿態(tài)的歐拉角表示

蘇召斌，林 森

（中交天津港航勘察設(shè)計研究院有限公司，天津 300457）

傳統(tǒng)的絞刀運動往往被簡化為橫移加旋轉(zhuǎn)的二維運動方式，這種運動簡化方法難以準(zhǔn)確地定義絞刀的位置和姿態(tài)。以剛體運動學(xué)為出發(fā)點，嚴(yán)格推導(dǎo)了絞刀運動的歐拉表示方法，以全新的視角描述了絞刀的運動。對繞任意軸轉(zhuǎn)動的歐拉旋轉(zhuǎn)進行了公式推導(dǎo)，能夠簡單而準(zhǔn)確地描述絞刀的位置和姿態(tài)變化。使用絞刀橫移偏轉(zhuǎn)角度、下放角度和旋轉(zhuǎn)角度定義絞刀的位置和姿態(tài)，并給出了絞刀和刀齒位置姿態(tài)的表達(dá)形式。

絞刀；歐拉角；坐標(biāo)變換

0 引言

在港口建設(shè)、航道疏浚及吹填工程中，絞吸式挖泥船被廣泛應(yīng)用[1]，絞刀是絞吸式挖泥船的核心裝備。絞吸挖泥船在估算挖泥產(chǎn)量時，將被挖泥土看成切層厚度、步進距離和擺寬的乘積，切層厚度指絞刀切削塊的豎向尺寸，步進距離指切削塊的縱向尺寸，擺寬就是指切削區(qū)間的橫向尺寸。事實上，絞刀的運動是三種轉(zhuǎn)動的疊加，分別為繞著定位鋼樁的轉(zhuǎn)動、繞著橋架耳軸的轉(zhuǎn)動和絞刀軸線的轉(zhuǎn)動[2]。絞刀距離定位鋼樁約為100 m，距離絞刀耳軸約40 m，而絞刀直徑僅僅約為3 m，經(jīng)過簡化以后，普遍將絞刀運動分解為沿著水平方向的平動和沿著絞刀軸線的轉(zhuǎn)動[3]。在絞吸挖泥船施工的一個擺寬區(qū)間，絞刀繞著定位鋼樁的擺動幅度將近，目前對絞刀的運動簡化與實際工作狀態(tài)并不相符。

目前，描述空間結(jié)構(gòu)姿態(tài)表示的有PRY角、歐拉角和四元數(shù)等，其中歐拉角是用來描述轉(zhuǎn)動剛體的3個獨立角參量，描述絞刀的運動更適合[4-6]。

1 坐標(biāo)變換

1.1 絞刀運動分析

絞刀整體運動一般被看做橫移和轉(zhuǎn)動的復(fù)合運動，其中絞刀距離鋼樁軸線的距離遠(yuǎn)大于絞刀的直徑，繞鋼樁的旋轉(zhuǎn)運動被簡化為橫移運動。如圖1所示，當(dāng)絞吸挖泥船向前推移一個步進距離，絞刀進入下一個往復(fù)運動過程。在每一個往復(fù)運動過程中，鋼樁均是在固定位置。如圖2所示，除了繞鋼樁轉(zhuǎn)動外，絞刀的運動還包括：繞橋架耳軸的旋轉(zhuǎn)運動和繞絞刀軸線的旋轉(zhuǎn)運動。因此，絞刀的運動是上述三種轉(zhuǎn)動的復(fù)合運動。

圖1 絞刀繞鋼樁運動示意圖Fig.1 Diagrammatic sketch of the cutter rotated around the steel pile

圖2 絞刀繞橋架耳軸轉(zhuǎn)動示意圖Fig.2 Diagrammatic sketch of the cutter rotated around the bridge trunnion

1.2 三維坐標(biāo)變換

剛體的運動可以等效為一次平動和一次轉(zhuǎn)動。剛體上的任意一點P在固定坐標(biāo)系A(chǔ)中的齊次坐標(biāo)為AP=（xA，yA，zA，1），在坐標(biāo)B系中的齊次坐標(biāo)為BP=（xB，yB，zB，1），則存在關(guān)系式：

式中：繞x軸旋轉(zhuǎn)的矩陣為Rx（琢）；繞y軸旋轉(zhuǎn)的矩陣為Ry（茁）；繞z軸旋轉(zhuǎn)的矩陣為Rz（酌）。對于右手坐標(biāo)系逆時針為正，即繞x軸旋轉(zhuǎn)時，由y軸轉(zhuǎn)向z軸為正；繞y軸旋轉(zhuǎn)時，由z軸轉(zhuǎn)向x軸為正；繞z軸旋轉(zhuǎn)時，由x軸轉(zhuǎn)向y軸為正。

1.3 定點旋轉(zhuǎn)的歐拉變換

歐拉角是用來唯一地確定定點轉(zhuǎn)動物體位置的3個一組獨立參量，由章動角茲，進動角鬃和自轉(zhuǎn)角漬組成。按照旋轉(zhuǎn)所繞軸次序的不同，共有12種不同的歐拉角，6種對稱型歐拉角為：XYX，XZX，YXY，YZY，ZXZ和 ZYZ。6種非對稱型歐拉角為：XYZ、XZY、YXZ、YZX、ZXY和ZYX。

如圖3所示的ZXZ型歐拉角，即坐標(biāo)軸先繞z軸旋轉(zhuǎn)鬃角，再繞新的x軸（即N軸）旋轉(zhuǎn)茲角，最后繞新的z軸（即z憶軸）旋轉(zhuǎn)漬角，剛體旋轉(zhuǎn)矩陣表示為：

同理，其它5種對稱型歐拉角對應(yīng)的旋轉(zhuǎn)矩陣可以按照旋轉(zhuǎn)的先后順序求解。

圖3 定點轉(zhuǎn)動的歐拉角示意圖Fig.3 Diagrammatic sketch of Euler angle rotated round the fixed point

當(dāng)歐拉角類型為ZXY時，即坐標(biāo)軸先繞y軸旋轉(zhuǎn)鬃角，再繞新的x軸旋轉(zhuǎn)茲角，最后繞新的z軸旋轉(zhuǎn)漬角，剛體旋轉(zhuǎn)矩陣為：

同理，其它5種非對稱型歐拉角對應(yīng)的旋轉(zhuǎn)矩陣可以按照旋轉(zhuǎn)的先后順序求解。

1.4 繞任意軸轉(zhuǎn)動的歐拉變換

記固定坐標(biāo)系為o-xyz，空間物體繞任意軸轉(zhuǎn)動的可以按照如下步驟進行：

1）將軸線平移到坐標(biāo)原點寅2）將軸線繞固定坐標(biāo)系的x軸旋轉(zhuǎn)琢角至xoz平面寅3）將軸線繞固定坐標(biāo)系的y軸旋轉(zhuǎn)茁角至z軸寅4）繞固定坐標(biāo)系的z軸旋轉(zhuǎn)茲角度寅5）將軸線繞固定坐標(biāo)系的x軸旋轉(zhuǎn)-茁角至xoz平面寅6）將軸線繞固定坐標(biāo)系的x軸旋轉(zhuǎn)-琢角至xoz平面寅7）將軸線平移到原位置。

則空間物體繞任意軸的歐拉變換可以看做一系列的旋轉(zhuǎn)變換，按照式（12）進行變換：

式中：T（a1，b1，c1，a2，b2，c2，茲）為旋轉(zhuǎn)變換矩陣，旋轉(zhuǎn)軸線用起點 p1（a1，b1，c1）和 p2（a2，b2，c2）終點來表示；茲為旋轉(zhuǎn)角度。軸線的方向可以表示為（a，b，c）：

同理，步驟3）為軸線繞y軸順時針旋轉(zhuǎn)-茁角度的變換：

其變換矩陣為：

步驟4）的變換矩陣參照式（6），步驟5）的變換矩陣為式（18）逆矩陣，步驟6）的變換矩陣為式（16）逆矩陣，步驟7）的變換矩陣為式（14）逆矩陣，則把上述7個矩陣順次連乘就是繞任意軸的歐拉變換矩陣：

2 絞刀的運動

2.1 坐標(biāo)系建立

以絞刀抬出水面并且絞刀軸指向船艏方向為基準(zhǔn)狀態(tài)，建立固定參考系M。其中以絞刀軸為z軸并以指向輪轂方向為正，以大圈下平面為xoy平面。同時，建立隨動坐標(biāo)系N，則絞刀在基準(zhǔn)狀態(tài)時，坐標(biāo)系N和坐標(biāo)系M重合，相對于船體而言，x軸朝右，y軸朝下，z軸朝前，如圖4所示。

圖4 絞吸挖泥船F(xiàn)ig.4 Cutter suction dredger

2.2 絞刀姿態(tài)確定

絞吸挖泥船施工時，絞刀首先被下放到一定深度處，然后繞著鋼樁水平擺動，同時絞刀也繞自身軸線轉(zhuǎn)動。為了表述方便，記絞吸挖泥船橫向擺動的偏轉(zhuǎn)角為鬃，一般在正負(fù)60毅之間；絞刀的下放角茲，一般小于50毅；絞刀的轉(zhuǎn)動角度為漬，可取任意角度。

絞刀的運動表示方式和歐拉角次序有關(guān)。絞刀在施工時，絞刀刀齒分別繞著鋼樁軸線、耳軸軸線和絞刀軸線運動。其中鋼樁軸線的優(yōu)先級最高，耳軸軸線次之，絞刀軸線最低，因此歐拉角的次序ZXY。只要給定具體的鬃、茲和漬的大小，就可以按照式（11）確定絞刀的姿態(tài)。

2.3 絞刀位置確定

記鋼樁軸線到耳軸軸線的距離為A，耳軸軸線到參考系M的距離為B，則位置矩陣為

綜上所述，絞刀的運動位置和姿態(tài)可以描述為：

式中：s茲表示 sin茲；c茲表示 cos茲；s鬃 表示 sin鬃；c鬃 表示 cos鬃；s漬 表示 sin漬；c漬 表示 cos漬。

同理，記任意刀齒初始時刻的柱坐標(biāo)為（R，準(zhǔn)，z），則刀齒的位置和姿態(tài)可以描述為：

式中：XA是轉(zhuǎn)動后刀齒在固定參考系M下的X坐標(biāo)；YA是轉(zhuǎn)動后刀齒在固定參考系M下的Y坐標(biāo)；ZA是轉(zhuǎn)動后刀齒在固定參考系M下的Z坐標(biāo)。

2.4 實例驗證

假設(shè) A=70，B=30，歐拉角為（10，-30，90），即絞刀向右偏轉(zhuǎn)10毅，下放絞刀20毅，轉(zhuǎn)動90毅，記絞刀上的刀齒柱坐標(biāo)為（1.5，20，0），笛卡爾坐標(biāo)為（1.410，0.513，0），絞刀大圈中心的坐標(biāo)為（0，0，0），將各參數(shù)帶入式（22）中，則偏轉(zhuǎn)后絞刀大圈中心的笛卡爾坐標(biāo)為（16.667，15.000，-5.477），該刀齒新的笛卡爾坐標(biāo)為（18.076，14.487，-5.477）。

3 結(jié)語

1）以剛體運動學(xué)為出發(fā)點，討論了在直角坐標(biāo)系中的歐拉坐標(biāo)變換方程。對于絞刀運動描述而言，該方法比傳統(tǒng)的橫移加旋轉(zhuǎn)的輪擺式運動方式更加準(zhǔn)確。

2）推導(dǎo)出繞任意軸的歐拉角旋轉(zhuǎn)方程，利用絞刀鋼樁-耳軸-絞刀軸的先后繞轉(zhuǎn)次序，推導(dǎo)出絞刀姿態(tài)和位置的關(guān)系表達(dá)形式。

3）根據(jù)旋轉(zhuǎn)次序和旋轉(zhuǎn)方向的不同，利用歐拉坐標(biāo)變換方法描述了絞刀刀齒的運動軌跡，給出了任意刀齒在任意歐拉角時的位置和姿態(tài)表達(dá)形式。

在給定絞刀的下放角度、轉(zhuǎn)動偏角和絞刀自轉(zhuǎn)角度時，能夠唯一確定絞刀的位置和姿態(tài)，同時也能確定每一個刀齒的位置。在后續(xù)的研究中，對軌跡方程進行時間求導(dǎo)，能夠得出各刀齒的速度，進而求出各刀齒的運動加速度，能夠為絞刀的優(yōu)化設(shè)計提供技術(shù)支持。

[1] 馮永軍，陳新權(quán)，丁金鴻，等.絞吸挖泥船橋架發(fā)展概述[J].中國港灣建設(shè)，2015（8）：1-5.FENG Yong-jun,CHEN Xin-quan,DING Jin-hong,et al.Development overview of bridge ladder on cutter suction dredgers[J].China Harbour Engineering,2015(8):1-5.

[2] 李金貴，李進軍，楊建華，等.污染底泥精確疏浚技術(shù)[J].中國港灣建設(shè)，2004（6）：11-14，20.LI Jin-gui,LI Jin-jun,YANG Jian-hua,et al.Technology for accurate dredging of contaminated subsoil[J].China Harbour Engi原neering,2004(6):11-14,20.

[3]王翰濤，魯力群，王秀景.絞吸式挖泥船絞刀系統(tǒng)技術(shù)現(xiàn)狀及發(fā)展[J].機電設(shè)備，2014（6）：10-16.WANG Han-tao,LU Li-qun,WANG Xiu-jing.Technology status and development of cutter system of cutter suction dredger[J].Mechanical and Electrical Equipment,2014(6):10-16.

[4]張帆，曹喜濱，鄒經(jīng)湘.一種新的全角度四元數(shù)與歐拉角的轉(zhuǎn)換算法[J].南京理工大學(xué)學(xué)報，2002,26（4）：376-380.ZHANG Fan,CAO Xi-bin,ZOU Jing-xiang.A new large-scale transformation algorithm of quaternion to Euler angle[J].Journal of Nanjing University of Science and Technology,2002,26(4):376-380.

[5] 林積昶，韓國柱，吳松，等.基于虛擬維修系統(tǒng)的配件模型實時移動方法[J].計算機工程，2010,36（10）：254-255.LIN Ji-chang,HAN Guo-zhu,WU Song,et al.Real-time move method for accessory model based on virtual maintenance system[J].Computer Engineering,2010,36(10):254-255.

[6] 徐旭松，楊將新，曹衍龍，等.基于齊次坐標(biāo)變換的制造誤差建模研究[J].浙江大學(xué)學(xué)報：工學(xué)版，2008，42（6）：1 021-1 026.XU Xu-song,YANG Jiang-xin,CAO Yan-long,et al.Manufactur原ing error modeling based on homogeneous transformation[J].Journal of Zhejiang University:Engineering Science,2008,42(6):1 021-1 026.

Expression of the moving orientation of cutter by Euler angle

SU Zhao-bin,LIN Sen
（CCCC Tianjin Port&Waterway Prospection&Design Research Institute Co.,Ltd.,Tianjin 300457,China）

The traditional sports of cutter is often simplified as two-dimensional motion with transverse and rotation,this simplified method can not accurately define the movement and posture of cutter.As the starting point of the rigid body kinematics,we rigorously derived the Euler expression of cutter movement,and described the movement of the cutter from a new perspective.We deduced the Euler rotation to any axis,and can describe the position and attitude of the cutter simply and effectively.Uses the sliding angle,downward angle and revolving angle can define cutter position and attitude,and gives the expression of the position and attitude of cutter and tooth.

cutter；Euler angle；coordinate transformation

U615.351

A

2095-7874（2017）12-0062-05

10.7640/zggwjs201712014

2017-05-17

2017-07-13

蘇召斌（1986— ），男，山東鄄城人，碩士，工程師，固體力學(xué)專業(yè)，主要從事疏浚挖泥機具研究和巖土切削數(shù)值模擬工作。E-mail：suprc@163.com