楊詩怡, 張峰峰, 范立成, 匡紹龍, 孫立寧

(蘇州大學(xué) 機電工程學(xué)院機器人與微系統(tǒng)研究中心, 江蘇 蘇州 215000)

?

放療床多目標(biāo)協(xié)調(diào)機構(gòu)參數(shù)優(yōu)化研究

楊詩怡, 張峰峰, 范立成, 匡紹龍, 孫立寧

(蘇州大學(xué) 機電工程學(xué)院機器人與微系統(tǒng)研究中心, 江蘇 蘇州 215000)

隨著精確放療技術(shù)在腫瘤治療中的廣泛應(yīng)用,放療床的研發(fā)也受到人們越來越多的關(guān)注.為了達(dá)到使放療床的工作空間和結(jié)構(gòu)尺寸能夠滿足實際使用要求的目的,提出針對6-HTRT并聯(lián)構(gòu)型放療床的一種多目標(biāo)協(xié)調(diào)的機構(gòu)參數(shù)優(yōu)化方法.研究過程中首先對放療床進(jìn)行了運動學(xué)分析,結(jié)合運動學(xué)逆解約束條件,對雙端虎克鉸形式的并聯(lián)放療床進(jìn)行工作空間的精確分析,并得到工作空間隨各結(jié)構(gòu)參數(shù)的變化規(guī)律,為放療床的參數(shù)優(yōu)化提供依據(jù).其次,針對提高虎克鉸強度及剛度問題,提出了一種虎克鉸初始安裝角度優(yōu)化方法.之后考慮工作空間、放療床受力、雅可比矩陣條件數(shù),同時構(gòu)建關(guān)于條件數(shù)全域均值、條件數(shù)波動情況的目標(biāo)函數(shù),對放療床進(jìn)行參數(shù)的具體優(yōu)化.最終得到了放療床各結(jié)構(gòu)參數(shù),經(jīng)驗證使得放療床的雅可比矩陣條件數(shù)在工作空間內(nèi)較小且分布均勻.由上述可知,所研究的放療床多目標(biāo)協(xié)調(diào)機構(gòu)參數(shù)優(yōu)化方法在一定范圍內(nèi)可以減小雅克比矩陣條件數(shù),使經(jīng)過參數(shù)優(yōu)化后的放療床得到合理的工作空間和結(jié)構(gòu)尺寸，滿足實際使用的要求.

放療床; 工作空間; 參數(shù)優(yōu)化; 雅可比矩陣

世界上每年因腫瘤死亡的人數(shù)約500萬.我國每年新發(fā)現(xiàn)的腫瘤患者約100萬人以上[1].目前,手術(shù)、放射和藥物治療是3種主要治療方法.隨著放射治療技術(shù)的快速進(jìn)步,腫瘤放療向著精確放療方向快速發(fā)展,從而取得了更好的治療效果[2].

常規(guī)放療床如西門子550TxT、CyberKnife放射治療床[3-4]均為串聯(lián)結(jié)構(gòu),只能實現(xiàn)3個自由度的運動,且剛度較差,定位精度不高,承載能力較差,空間自由度受到限制.而并聯(lián)機構(gòu)由于剛度大、精度高、動態(tài)性能好、承載能力強等優(yōu)點,在機器人領(lǐng)域得到了越來越廣泛的重視[5],本文提出的六自由度放療床可以有效滿足放療設(shè)備使用要求,達(dá)到精確放療的效果.

放療床用于對患者進(jìn)行定位,因此對于其平穩(wěn)性和安全性有嚴(yán)格要求.六自由度并聯(lián)機構(gòu)結(jié)構(gòu)參數(shù)復(fù)雜,工作空間小,因此,通過分析各參數(shù)對放療床受力的影響,對放療床進(jìn)行結(jié)構(gòu)設(shè)計和參數(shù)優(yōu)化就成為極其重要的環(huán)節(jié).針對放療床的使用要求,本文提出一種多目標(biāo)協(xié)調(diào)的機構(gòu)參數(shù)優(yōu)化方法,考慮工作空間、放療床受力、雅克比矩陣條件數(shù),以及構(gòu)建關(guān)于條件數(shù)全域均值、條件數(shù)波動情況的目標(biāo)函數(shù),進(jìn)行參數(shù)的具體優(yōu)化.

1　運動學(xué)模型

根據(jù)精確放射治療的要求,選用如圖1所示的6-HTRT結(jié)構(gòu)形式.該構(gòu)型沿著導(dǎo)軌方向的運動即縱向運動的范圍僅受導(dǎo)軌長度限制,可以輕易得到放大,且結(jié)構(gòu)整體寬度較小,可以與床板寬度保持一致,結(jié)構(gòu)的整體體積也較小.因此,該并聯(lián)構(gòu)型能夠很好地滿足放療床的使用要求.

圖1　6-HTRT放射治療床Fig.1　Radiotherapy bed of 6-HTRT

1.1坐標(biāo)系建立

圖2　放療床坐標(biāo)系Fig.2　Coordinate system of radiotherapy bed

設(shè)動平臺中心到鉸鏈中心的矢徑與該組鉸鏈軸線間的夾角為θp,將θp定義為動平臺鉸鏈分布角度.

用Φpi表示動平臺坐標(biāo)系{P}中x正半軸到矢徑PPi沿逆時針方向需要旋轉(zhuǎn)的角度，則有下述公式:

(1)

定平臺虎克鉸中心點Bi沿著基座坐標(biāo)系的y軸方向作平移運動,將其移動量Si作為放療床的輸入量.動平臺虎克鉸中心點Pi在動平臺坐標(biāo)系中的矢量描述為

(2)

式中：r為動平臺半徑,mm.

定平臺鉸鏈中心點Bi在{B}中的矢量描述為

(3)

xBi與yBi的表達(dá)式如下:

(4)

式中：λ為相鄰2條導(dǎo)軌中心線之間的距離,mm；

ki與動平臺中心點的空間位置和姿態(tài)有關(guān).

設(shè)動平臺坐標(biāo)系繞基座坐標(biāo)系3個主軸的轉(zhuǎn)角分別為θx,θy,θz,得到的合成旋轉(zhuǎn)變換矩陣來表示動平臺坐標(biāo)系相對于基座坐標(biāo)系的姿態(tài),具體表達(dá)式如式(5):

(5)

式中:Sθx=sinθx,Sθy=sinθy,Sθz=sinθz,

Cθx=cosθx,Cθy=cosθy,Cθz=cosθz.

(6)

1.2運動學(xué)逆解

在對放療床進(jìn)行位置控制時,當(dāng)給定床板在空間內(nèi)某一位姿,需要利用其位姿逆解來得到各個滑塊的輸入值,進(jìn)而得到電機輸入的參考值.對于本文研究的放療床,當(dāng)給定動平臺幾何中心的位置和姿態(tài)后,通過結(jié)構(gòu)的幾何關(guān)系求出各個滑塊的位移,進(jìn)而得到電機的轉(zhuǎn)角，即為其位姿逆解.設(shè)點Pi在基座坐標(biāo)系中的位置矢量為

(7)

(8)

設(shè)向量P為動平臺中心點在基座坐標(biāo)系中的矢量描述,由此可得動平臺虎克鉸中心點在基座坐標(biāo)系中的矢量描述為

BPi=BQP·pPi+P.

(9)

根據(jù)上述表達(dá)式可得到放療床的位姿逆解yBi,對于表達(dá)式中的正負(fù)號的選取,結(jié)合了動平臺幾何中心的位置和各滑塊所在的位置.

1.3影響工作空間的因素

放療床的位置工作空間可描述為,在不考慮動平臺姿態(tài)轉(zhuǎn)角情況下,即動平臺平動時動平臺中心點能夠到達(dá)的最大空間范圍[9-10].

1)影響放療床工作空間的主要結(jié)構(gòu)參數(shù):r,為動平臺虎克鉸中心點Pi所在圓的半徑;l(桿長參數(shù)),為定平臺虎克鉸中心點Bi與動平臺虎克鉸中心點Pi間的距離;λ,為相鄰導(dǎo)軌間的距離;θp,為動平臺虎克鉸鉸鏈分布角度.機構(gòu)的工作空間可由r,l,λ,θp完全確定,因此,對以上4個參數(shù)進(jìn)行初選、調(diào)整和優(yōu)化來滿足所需位置和姿態(tài)工作空間以及放療床性能的要求.

2)逆解約束條件包括:①各虎克鉸的轉(zhuǎn)角應(yīng)在允許的轉(zhuǎn)動范圍之內(nèi),即動平臺中心運動到工作空間內(nèi)的各個位姿,各虎克鉸之間均不發(fā)生干涉;②連桿之間不發(fā)生干涉,即動平臺在運動過程中相鄰兩桿之間的最短距離應(yīng)大于連桿的直徑[11];③定平臺滑塊的移動應(yīng)在允許的行程范圍內(nèi);④放療床凸性限制,即各滑塊沿導(dǎo)軌方向應(yīng)在其相應(yīng)動平臺鉸鏈點的外側(cè).

2　放療床參數(shù)優(yōu)化

由于受到放射治療使用過程和系統(tǒng)整體布局的限制,要求放療床床板寬度及整體外形具有較小的尺寸,但并聯(lián)機構(gòu)本身工作空間較小,因此,需要通過優(yōu)化參數(shù)來滿足工作空間.基于以上限制,本文提出一種多目標(biāo)協(xié)調(diào)機構(gòu)參數(shù)優(yōu)化方法,考慮工作空間、放療床受力情況、雅可比矩陣條件數(shù)[12-13],以及構(gòu)建關(guān)于條件數(shù)全域均值、條件數(shù)波動情況的目標(biāo)函數(shù),進(jìn)行參數(shù)優(yōu)化.

2.1虎克鉸初始安裝角度優(yōu)化

在動平臺所需的工作空間確定之后,可以通過優(yōu)化虎克鉸的初始安裝角度的方法來減小虎克鉸繞自身動坐標(biāo)系x′軸和y′所需要的轉(zhuǎn)角極限值,這樣可以有效減小虎克鉸的尺寸,進(jìn)一步提高虎克鉸的強度和剛度,在相同的受力情況下減小虎克鉸的變形,提高放療床的整體剛度[14-15].

先假設(shè)定平臺各虎克鉸繞基座坐標(biāo)系x,z軸的初始安裝角度均為零度,式(10)所示為動平臺處于任意位姿時定平臺每對虎克鉸繞其自身坐標(biāo)系的擺角α和β:

(10)

因此,遍歷放療床整個使用工作空間,求得各個位姿的定平臺虎克鉸的2個擺角,并找出每個虎克鉸的擺角在整個遍歷空間內(nèi)所取得的最大值和最小值.在此,取最大與最小值的平均值作為虎克鉸的初始安裝角度,公式如下:

(11)

(12)

得到定平臺虎克鉸的初始安裝角度后,進(jìn)行相應(yīng)的圓整,數(shù)據(jù)見表1.

表1定平臺虎克鉸初始安裝角度

Table 1Initial mounting angle of the Hooke joint of fixed platform

虎克鉸序號繞x軸初始安裝角度/(°)繞z軸初始安裝角度/(°)定平臺虎克鉸1-412定平臺虎克鉸2412定平臺虎克鉸335-7.5定平臺虎克鉸4355.5定平臺虎克鉸5-355.5定平臺虎克鉸6-35-7.5

同理,按照上述方法得到動平臺各對虎克鉸的初始安裝角度后,進(jìn)行相應(yīng)的圓整,數(shù)據(jù)見表2.

表2動平臺虎克鉸初始安裝角度

Table 2Initial mounting angle of the moving platform Hooke joint

虎克鉸序號繞x軸初始安裝角度/(°)繞z軸初始安裝角度/(°)動平臺虎克鉸1-41-2動平臺虎克鉸241-2動平臺虎克鉸3347.5動平臺虎克鉸434-5.5動平臺虎克鉸5-34-5.5動平臺虎克鉸6-347.5

沒考慮虎克鉸初始安裝角度時,在工作空間內(nèi),以動平臺各對虎克鉸為例,其繞自身坐標(biāo)系軸線的轉(zhuǎn)角范圍如表3所示.

表3動平臺虎克鉸繞自身軸線的轉(zhuǎn)角范圍(不考慮虎克鉸初始安裝角度情況下)

Table 3Angle range of the moving platform Hooke joint around it’s own axis (without considering Hooke joint initial mounting angle)

虎克鉸序號繞x'軸的轉(zhuǎn)角范圍/(°)繞y'軸的轉(zhuǎn)角范圍/(°)動平臺虎克鉸1-72.2～-8.9-24.6～20.8動平臺虎克鉸28.9～72.2-24.6～20.8動平臺虎克鉸3-5.1～72.2-17.7～32.6動平臺虎克鉸4-3.9～72.8-29.7～18.7動平臺虎克鉸5-72.8～3.9-29.7～18.7動平臺虎克鉸6-72.2～5.1-17.7～32.6

表4動平臺虎克鉸繞自身軸線的轉(zhuǎn)角范圍(考慮虎克鉸初始安裝角度)

Table 4Angle range of the moving platform Hooke joint around it’s own axis (considering Hooke joint initial mounting angle)

虎克鉸序號繞x'軸的轉(zhuǎn)角范圍/(°)繞y'軸的轉(zhuǎn)角范圍/(°)動平臺虎克鉸1-31.2～32.1-22.6～22.8動平臺虎克鉸2-32.1～31.2-22.6～22.8動平臺虎克鉸3-39～38.2-25.5～25.1動平臺虎克鉸4-37.9～38.8-24.2～24.2動平臺虎克鉸5-38.8～37.9-24.2～24.2動平臺虎克鉸6-38.2～39.1-25.2～25.1

通過上述分析得到,虎克鉸繞z軸的初始安裝角度較小,考慮繞z軸的初始安裝角度后虎克鉸的轉(zhuǎn)角極限減小程度較小,且會給零件的加工和裝配帶來困難,因此,本文只考慮虎克鉸繞x軸的初始安裝角度,在保證滿足工作空間的條件下,減小了虎克鉸繞自身坐標(biāo)系x軸的轉(zhuǎn)角極限值,因此,可以有效減小虎克鉸的結(jié)構(gòu)尺寸.優(yōu)化后虎克鉸尺寸參數(shù)如表5所示.

表5優(yōu)化后虎克鉸參數(shù)

Table 5　Optimized Hooke joint parameters　mm

2.2放療床幾何參數(shù)初步優(yōu)化

2.2.1動平臺半徑r

首先給定各幾何參數(shù)一個初值,使這幾個參數(shù)構(gòu)成的放療床模型滿足其所需工作空間的要求.然后,固定除動平臺半徑r以外的其他幾何參數(shù),進(jìn)一步通過工作空間分析可得,增大r可以增大放療床的工作空間.但由于其實際的使用要求,床板寬度不能超過800 mm,進(jìn)而r的取值受到一定的限制.

利用上述指標(biāo)對放療床的r進(jìn)行優(yōu)化.如圖3(a)所示,增大動平臺半徑r可以減小放療床雅可比矩陣的條件數(shù),且動平臺的高度越大,機構(gòu)的雅可比矩陣條件數(shù)越小.

桿件所受的最大力在一定范圍內(nèi)隨著動平臺半徑的增大而減小,如圖3(b)所示.因此,在床板寬允許的范圍內(nèi),增大放療床的動平臺半徑,可以起到增大工作空間、減小機構(gòu)雅可比矩陣條件數(shù)以及減小桿件所受最大力的作用.進(jìn)而,可以初步確定r較好的取值范圍為300～400 mm.

圖3　動平臺半徑對放療床性能影響Fig.3　Influence of moving platform radius on radiotherapy bed performance

2.2.2導(dǎo)軌間距λ

這種模式較適于大型肉牛場。建設(shè)沼氣系統(tǒng)，以糞污混合物、污水為原料生產(chǎn)沼氣，沼渣、沼液用于農(nóng)田或深加工成有機肥和液體肥，沼氣可發(fā)電，也可作為燃?xì)馐褂?。大?guī)模沼氣可考慮發(fā)電并網(wǎng)或生產(chǎn)生物天然氣，中等規(guī)?？煽紤]電、氣自用，或者供附近百姓使用。該模式優(yōu)點是生產(chǎn)綠色能源；缺點是周圍要有消納沼渣、沼液的農(nóng)田，如再建設(shè)沼渣、沼液長期貯存設(shè)施，投資較高。

固定除導(dǎo)軌間距λ以外的其他幾何參數(shù),進(jìn)一步通過工作空間分析可得,當(dāng)導(dǎo)軌間距λ略小于動平臺半徑時,放療床的工作空間較大.如圖4(a)所示,導(dǎo)軌間距λ越小,動平臺的高度越高,放療床的雅克比矩陣條件數(shù)越小.

桿件所受的最大力在一定范圍內(nèi)隨導(dǎo)軌間距的變化規(guī)律如圖4(b)所示,導(dǎo)軌間距在200～300 mm范圍內(nèi)變化時,桿件所受的最大力較小.因此,可以確定導(dǎo)軌間距較好的取值范圍為200～300 mm.

圖4　導(dǎo)軌間距對放療床性能影響Fig.4　Influence of guide spacing on radiotherapy bed performance

2.2.3動平臺鉸鏈分布角度θp

固定除θp以外的其他幾何參數(shù),進(jìn)一步通過工作空間分析可得,θp對放療床工作空間的影響很小,但主要對工作空間內(nèi)雅可比矩陣條件數(shù)大小及分布產(chǎn)生影響.如圖5(a),30°～90°與50°～70°的鉸鏈分布角度時雅可比矩陣條件數(shù)較小,并且隨著動平臺中心高度的增加條件數(shù)逐漸減小.在MATLAB中得到放療床30°～90°的分布角度比50°～70°分布桿件所受的最大力要小,因此初步選30°～90°作為放療床的鉸鏈分布角度.

圖5　動平臺鉸鏈分布角度和桿長對放療床性能的影響Fig.5　Influence of moving platform hinge distribution angle and rod length on radiotherapy bed performance

2.2.4桿長l

固定除桿長l以外的其他幾何參數(shù),進(jìn)一步通過工作空間分析可得,增大桿長l可以較大程度地增大放療床的工作空間,但結(jié)構(gòu)的整體尺寸會變得很大.根據(jù)放療床的使用要求,床板前段下方要留出一定的空間放置影像板,用于X射線的成像,以便實時跟蹤病灶部位,因此導(dǎo)軌長度受到制約,此情況下,在一定范圍內(nèi)減小連桿長度反而可以增大放療床的工作空間.結(jié)合工作空間、雅可比矩陣條件數(shù)以及連桿受力對l進(jìn)行優(yōu)化顯得尤為重要,如圖5(b),連桿長度在600～700 mm范圍內(nèi)時,連桿所受的最大力較小.因此,可以確定連桿長度較好的取值范圍為600～700 mm.

2.3放療床幾何參數(shù)最終確定

以上通過對工作空間、雅可比矩陣條件數(shù)和連桿受力情況進(jìn)行的分析,得到放療床各幾何參數(shù)較佳的取值范圍.在此,通過建立目標(biāo)函數(shù)對各參數(shù)進(jìn)行最終的確定[16].將雅可比矩陣條件數(shù)的全域均值作為放療床性能的評價指標(biāo),其反映了放療床在整個工作空間內(nèi)性能的平均水平[17-18],具體表達(dá)式如式(13)所示，

(13)

(14)

圖6　機構(gòu)A與機構(gòu)B性能比較Fig.6　Performance comparison between institution A and institution B

(15)

通過MATLAB編程,對上述得到的各參數(shù)較佳的取值范圍進(jìn)行遍歷,求解各參數(shù)組合下ηs的值.此時,參數(shù)優(yōu)化變?yōu)榍蠼饽繕?biāo)函數(shù)的最小值問題,即

(16)

通過上述優(yōu)化得到各參數(shù)的取值,對得到的值進(jìn)行圓整,最終得到放療床的參數(shù)，見表6.

表6　放療床最終參數(shù)

2.4放療床工作空間截面內(nèi)雅可比矩陣條件數(shù)分布

根據(jù)最終確定的幾何參數(shù),在MATLAB中得到放療床雅可比矩陣條件數(shù)在工作空間截面內(nèi)的分布,如圖7.由此可以看出,雅可比矩陣條件數(shù)在工作空間內(nèi)較小且總體分布比較均勻,只有在工作空間截面內(nèi)的邊界部位，雅可比矩陣的條件數(shù)稍大些.

3　結(jié)　論

為了實現(xiàn)放療床的精確治療,本文提出了一種多目標(biāo)協(xié)調(diào)的機構(gòu)參數(shù)優(yōu)化方法,對參數(shù)進(jìn)行具體優(yōu)化,最終得到了放療床各參數(shù),使得放療床在工作空間內(nèi)雅可比矩陣條件數(shù)較小且分布均勻.

本文首先建立了放療床的運動學(xué)模型,分析了各對虎克鉸的運動范圍,求得了放療床在確定位姿下各虎克鉸的確切擺角.對雙端虎克鉸形式的并聯(lián)放療床進(jìn)行了工作空間的精確分析,得到了機構(gòu)各參數(shù)對工作空間的影響,為放療床的參數(shù)優(yōu)化提供依據(jù).之后提出了一種優(yōu)化虎克鉸初始安裝角度的方法,在滿足工作空間的條件下,可以有效減小虎克鉸繞各軸的轉(zhuǎn)角極限,進(jìn)而減小虎克鉸的尺寸.根據(jù)放療床使用特性的要求,提出了一種多目標(biāo)協(xié)調(diào)的機構(gòu)參數(shù)優(yōu)化方法,考慮工作空間、放療床受力、雅可比矩陣條件數(shù),以及構(gòu)建關(guān)于條件數(shù)全域均值、條件數(shù)波動情況的目標(biāo)函數(shù),進(jìn)行參數(shù)的具體優(yōu)化,最終得到了放療床各參數(shù),使得放療床在工作空間內(nèi)雅可比矩陣條件數(shù)較小且分布均勻.

圖7　雅可比矩陣條件數(shù)分布Fig.7　Jacobian matrix condition number distribution

[1] 楊瑞杰,王俊杰.放射治療質(zhì)量保證與患者安全：2013北美放射腫瘤年會進(jìn)展[J].中國醫(yī)療設(shè)備,2014(29):163-164．YANG Rui-jie,WANG Jun-jie.Radiotherapy quality assurance and patient safety：the progress of 2013 ASTRO annual meeting[J].China Medical Device,2014(29):163-164．

[2] 甄鑫.自適應(yīng)放射治療中的圖像變形配準(zhǔn)關(guān)鍵技術(shù)研究[D].廣州:南方醫(yī)科大學(xué)生物醫(yī)學(xué)工程學(xué)院,2013:1-10．

ZHEN Xin.Research of deformable image registration for adaptive radiation therapy[D].Guangzhou:Southern Medical University,School of Biomedical Engineering,2013:1-10．

[3] SPEZI E,FERRI A.Dosimetric characteristics of the Simens IGRT carbon fiber tabletop[J].Med Dos,2007,32(4):295-298．

[4] SINCLAIR J,CHANG S D,GIBBS I C,et al.Multisession Cyberknife radiosurgery for intramedullary spinal cord arteriovenous malformations[J].Neurosurgery,2006,58(6):1081-1089．

[5] 彭忠琦.并聯(lián)機構(gòu)的發(fā)展及應(yīng)用[J].光機電信息,2011,28(12):45-50．

PENG Zhong-qi.Development and application of parallel mechanism[J].Ome Information,2011,28(12):45-50．

[6] 黃鵬,汪勁松,王立平,等.3-PRS并聯(lián)機構(gòu)誤差運動學(xué)分析及辨識[J].清華大學(xué)學(xué)報(自然科學(xué)版),2010,50(11):1811-1814．HUANG Peng,WANG Jing-song,WANG Li-ping,et al.Kinematical error analysis and identification of a 3-PRS parallel mechanism[J].Journal of Tsinghua University (Science & Technology),2010,50(11):1811-1814．

[7] 謝志江,史浩明.6-HUS并聯(lián)機構(gòu)位置逆解與運動學(xué)優(yōu)化設(shè)計[J].機械設(shè)計,2011,12(12):26-30．XIE Zhi-jiang,SHI Hao-ming.Inverse positions solution and kinematics optimal design based on a 6-HUS parallel mechanism[J].Journal of Machine Design,2011,12(12):26-30．

[8] 皮陽軍,王驥,胡玉梅.并聯(lián)六自由度機構(gòu)運動學(xué)與動力學(xué)標(biāo)定對比[J].哈爾濱工程大學(xué)學(xué)報,2014,35(11):1422-1426．

PI Yang-jun,WANG Ji,HU Yu-mei.Comparison of dynamic and kinematic calibrations for the 6-DOF parallel mechanism[J].Journal of Harbin Engineering University,2014,35(11):1422-1426．

[9] 于凌濤,孫立寧,杜志江.并聯(lián)機器人胡克鉸工作空間的研究與應(yīng)用[J].機械工程學(xué)報,2006,42(8):120-124．YU Ling-tao,SUN Li-ning,DU Zhi-jiang.Study and application of workspace on Hooke joint in parallel robot[J].Journal of Mechanical Engineering,2006,42(8):120-124．

[10] 劉玉斌,趙杰,楊永剛,等.一種新型6-PRRS并聯(lián)機器人工作空間分析[J].機械與電子,2007(2):53-56．LIU Yu-bin,ZHAO Jie,YANG Yong-gang,et al.Workspace analysis of a novel 6-PRRS parallel robot[J].Machinery & Electronics,2007(2):53-56．

[11] DASGUPTA B,MRUTHYUNJAYS T S.The stewart platform manipulator:a review[J].Mech Mach Theory,2000,35(1):15-40．

[12] 孫立寧,于暉,祝宇虹,等.機構(gòu)影響系數(shù)和并聯(lián)機器人雅克比矩陣的研究[J].哈爾濱工業(yè)大學(xué)學(xué)報,2002,34(6):810-814．

SUN Li-ning,YU Hui,ZHU Yu-hong,et al.Influence coefficient of mechanism and Jacobia matrix of parallel robot[J].Journal of Harbin Institute of Technology,2002,34(6):810-814．

[13] 蘭陟,李振亮,李亞,等.基于旋量理論的5-DOF上肢康復(fù)機器人雅克比矩陣求解[J].機械設(shè)計,2011,28(5):51-53,74．

LAN Zhi,LI Zhen-liang,LI Ya,et al.Calculation of Jacobin matrix of a 5-DOF upper limb re-habilitation robot based on screw theory[J].Journal of Machine Design,2011,28(5):51-53,74．

[14] 蓋永軍,王靜泉.虎克鉸在6-UPU運動平臺中的應(yīng)用[J].電子測試,2014(z1):52-54．

GAI Yong-jun,WANG Jing-quan.Application of Hooke Joints in 6-UPU motion platform[J].Electronic Test,2014(z1):52-54．

[15] 趙強,閻紹澤.雙端虎克鉸型六自由度并聯(lián)機構(gòu)的動力學(xué)模型[J].清華大學(xué)學(xué)報(自然科學(xué)版),2005,45(5):610-613．

ZHAO Qiang,YAN Shao-ze.Dynamic model of a 6-DOF parallel mechanism with Hooke’s joints at both chain ends[J].Journal of Tsinghua University(Science and Technology),2005,45(5):610-613．

[16] 孫小勇,鄭彬,鮑捷,等.高速6-PSS并聯(lián)機器人參數(shù)優(yōu)化設(shè)計[J].農(nóng)業(yè)機械學(xué)報,2015,46(5):372-378．

SUN Xiao-yong,ZHEN Bin,BAO Jie,et al.Parameter optimization design of high-speed 6-PSS parallel robot[J].Transactions of the Chinese Society for Agricultural Machinery,2015,46(5):372-378．

[17] 黃府,劉會議.基于雅可比矩陣條件數(shù)的并聯(lián)機構(gòu)參數(shù)優(yōu)化[J].機械,2012,39(12):41-45．

HUANG Fu,LIU Hui-yi.Parameters optimization of parallel mechanism based on the condition of Jacobian matrix[J].Machinery,2012,39(12):41-45．

[18] 魏永庚,袁明,石勇,等.基于雅克比矩陣的并聯(lián)機構(gòu)結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計[J].黑龍江大學(xué)工程學(xué)報,2015(1):87-91．

WEI Yong-kang,YUAN Ming,SHI Yong,et al.Structure parameters optimization design of parallel mechanism based on Jacobian matrix[J].Journal of Engineering of Heilongjiang University,2015(1):89-91．

Multi target coordinated mechanism parameter optimization method for radiotherapy bed

YANG Shi-yi, ZHANG Feng-feng, FAN Li-cheng, KUANG Shao-long, SUN Li-ning

(Robotics and Microsystems Center, College of Mechanical and Electrical Engineering, Soochow University, Suzhou 215000, China)

With the wide application of precise radiotherapy techniques in the treatment of tumor,the research of radiotherapy bed is receiving more and more attention.To achieve the aim that the workspace and the dimensions of the structure of the radiotherapy bed can meet the requirement in actual use,the multi target coordinated mechanism parameter optimization method for radiotherapy bed was introduced.In the research,the kinematics analysis of the radiotherapy bed was firstly carried out,the precise analysis of the workspace was taken on the double universal joints parallel radiotherapy bed,combining with the other inverse solution constraints.Then the influence of workspace by changing the rule of geometric parameters was obtained,which could provide the basis for the parameter optimization of the radiotherapy bed.For the problem of improving the strength and rigidity of the Hooke joint,a method of optimizing the initial mounting angle of the Hooke joint was proposed.Considering the workspace,the forces on it and the Jacobian matrix condition number,the objective function of the universe average of condition number and the fluctuation of condition number was constructed at the same time,the specific optimization of radiotherapy bed parameters was carried out.Finally,the parameters of the radiotherapy bed were obtained with the smaller condition number with uniform distribution of the Jacobian matrix in the whole workspace.According to the above,the multi target coordinated parameters optimization method for radiotherapy bed can reduce the Jacobian matrix condition number within limits.With the method,the radiotherapy bed can get the reasonable workspace and structure after the parameters are optimized,and can meet the requirement of actual use.

radiotherapy bed; workspace; parameter optimization; Jacobian matrix

2015-10-27.

國家高技術(shù)研究發(fā)展計劃(863計劃)資助項目(2012AA041602).

楊詩怡(1993—),女,江蘇蘇州人,碩士生,從事醫(yī)療機器人技術(shù)研究,E-mail:525406324@qq.com.

通信聯(lián)系人：張峰峰(1979—),男,副教授,博士,從事醫(yī)療機器人技術(shù)研究,E-mail:zhangfengfeng@suda.edu.cn.http://orcid.org//0000-0001-5250-9421

10.3785/j.issn. 1006-754X.2016.03.010

TP 242

A

1006-754X(2016)03-0256-08

本刊網(wǎng)址·在線期刊：http://www.journals.zju.edu.cn/gcsjxb