1.一種微創(chuàng)外科手術機器人遠心定位執(zhí)行機構設計方法,所述執(zhí)行機構包括第一連桿
(1)、第二連桿(2)和第三連桿(3),第一連桿(1)的一端與基座連接,第一連桿(1)的另一端
與第二連桿(2)的一端連接,第三連桿(3)的一端與第二連桿(2)的另一端連接,第一連桿
(1)的橫截面是由兩個正方形組成的回字形,第一連桿(1)的外徑a為81mm,第一連桿(2)的
內徑b為71mm,第三連桿(3)的橫截面是由兩個正方形組成回字形,第三連桿(3)的內徑a1為
73mm,第三連桿(3)的外徑b1為63mm,第二連桿(2)的橫截面為工字形,所述工字形的高度h
為81mm,所述工字形頂邊的長度bb為70mm,所述工字形中間的寬度t2為6mm,所述工字形頂
邊的寬度t1為6mm;其特征在于:所述一種微創(chuàng)外科手術機器人遠心定位執(zhí)行機構設計方法
是通過如下步驟實現(xiàn)的:
步驟一、遠心定位執(zhí)行機構的正運動學公式為:
公式(1)中P表示遠心定位執(zhí)行機構執(zhí)行端的位置坐標,P
x表示遠心定位執(zhí)行機構執(zhí)行
端在X軸上的坐標,P
y表示遠心定位執(zhí)行機構執(zhí)行端在Y軸上的坐標,P
z表示遠心定位執(zhí)行機
構執(zhí)行端在Z軸上的坐標,d表示遠心定位執(zhí)行機構執(zhí)行端的平移量,c
1=cosθ
1,θ
1表示旋轉
關節(jié)1運動角度,s
1=sinθ
1,c
2=cosθ
2,s
2=sinθ
2,θ
2表示表示旋轉關節(jié)2運動角度,cα
0=cos
α
0,sα
0=sinα
0,α
0=(γ
0+0.5π),γ
0表示旋轉關節(jié)1與旋轉關節(jié)2軸線夾角,cα
1=cosα
1,sα
1=sinα
1,α
1表示表示旋轉關節(jié)1與旋轉關節(jié)2軸線夾角,cα
2=cosα
2,sα
2=sinα
2,α
2表示表示
旋轉關節(jié)2與平移關節(jié)3軸線夾角;
步驟二、遠心定位執(zhí)行機構的逆運動學公式為:
k=(p
ysα
0-p
zcα
0+dcα
1cα
2)/(dsα
1sα
2)......(3)
θ
2=arccos(k)orθ
2=2π-arccos(k)......(4)
θ
1=acot(k
1/k
2)orθ
1=pi+acot(k
1/k
2)......(5)
k
1=p
xsα
2s
2-(p
ycα
0+p
zsα
0)(cα
1sα
2c
2+sα
1cα
2)......(6)
k
2=p
x(cα
1sα
2c
2+sα
1cα
2)+scα
1sα
2c
2+sα
2s
2(p
ycα
0+p
zsα)......(7)
公式(2)、(3)、(4)、(5)、(6)、(7)中,d表示遠心定位執(zhí)行機構執(zhí)行端的平移量,θ
1表示旋
轉關節(jié)1運動角度,θ
2表示表示旋轉關節(jié)2運動角度;
步驟三、遠心定位執(zhí)行機構的雅克比矩陣為:
公式(8)中J表示遠心定位機構前兩個旋轉運動的雅克比矩陣;;
步驟四、優(yōu)化遠心定位執(zhí)行機構的靈活性全局性能指標和全局剛度性能指標:
靈活性性能指標為:
公式(9)、(10)中m表示雅克比矩陣的行秩,n表示雅克比矩陣的列秩,ω表示正矩陣,tr
表示表示矩陣的跡,;
靈活性全局性能指標為:
ζ=∫
wk
Jdw/∫
wdw......(11)
公式(11)中ζ表示靈活性全局性能指標,k
J表示靈活性能指標,w表示工作空間;;
串聯(lián)機械臂第i部分基于基坐標的等效柔度矩陣為:
公式(12)中,
![]()
表示串聯(lián)機械臂第i部分基于基座標的等效柔度矩陣,
![]()
表示機械臂
第i部分相對于基座標的旋轉變換矩陣,
iP
T表示串聯(lián)機械臂第i部分基于基座標的位置向
量,P=(P
x,P
y,P
z),(P×)=(0 -P
z P
y;P
z 0 -P
x;-P
y P
x 0),
![]()
表示串聯(lián)機械臂第i部分相對
自身參考坐標系的等效柔度矩陣;
整體的剛度矩陣為:
公式(13)中K表示剛度矩陣,B表示基座標,i表示第i個坐標系,C表示柔度矩陣,K
t表示3
×3平移剛度矩陣,K
tr表示3×3平移旋轉耦合剛度矩陣,K
rt表示3×3旋轉平移剛度矩陣,K
r表示3×3旋轉剛度矩陣;
平移剛度和旋轉剛度條件數指標分別為:
k
t=||K
t||||(K
t)
-1||......(14)
k
r=||K
r||||(K
r)
-1||......(15)
公式(14)和(15)中
![]()
i表示t,r;
全局平移剛度和旋轉剛度指標分別為:
公式(16)和(17)中S
t表示全局平移剛度指標,S
r表示全局旋轉剛度指標;
全局剛度綜合性能指標為:
S=(S
r+S
t)/2......(18);
步驟五、計算兩種工作模式下的推導碰撞幾率:
設l
r(l
r表示切口到人的肩部的長度)的最小長度為660mm,根據公式計算ll:
ll=l
rcosγ
1......(19)
γ
1∈(0 1/3π)......(20)
公式(19)和(20)中γ
1表示旋轉關節(jié)1與水平面的夾角,由公式(19)和公式(20)可推出
當γ
1越大,ll就越小;
第一種工作模式的極限情況下,η可由以下公式求出:
公式(21)中η表示第一種工作模式下第一連桿與豎直面夾角,w
2表示人頭部寬度,l表示
第一連桿長度在豎直平面的投影長度,h
1表示基座離人頭部的高度,h
2表示人頭部高度,患
者的最大肩寬為489mm,如果l≥250mm且η≤θ
1≤π時,患者和機器人碰撞可以規(guī)避;l
r1(l
r1表
示第二旋轉關節(jié)與切口的距離)由如下公式求出:
l
r1=l/cosα
1+llcosα
1......(22)
公式(22)中α表示第一旋轉關節(jié)與第二旋轉關節(jié)的夾角;
根據手術工作空間的要求,α
2最小為0.25π,α
2最大為0.5π,α
2的平均值為0.375π,l
r1中
間值為395mm;當θ
2等于0.5π時,θ
1分別是β
1、β
2和α
2時的工作狀態(tài),當第三個平移軸在虛線左
側時,不會出現(xiàn)機器人與機器人之間的碰撞,當在右側會發(fā)生碰撞;但是當π-α
2≤θ
1≤π時,
θ
2≠90°時,就不會發(fā)生碰撞;因此為了避免碰撞,θ
2需要旋轉一定的角度,所要旋轉的角度
由下式求的:
B'D=ABsin(α
2)......(23)
CD=-ABcos(α
2)tan(θ
1)......(24)
∠B'DC=CD/B'D=-cot(α
2)tan(θ
1)......(25)
公式(23)、(24)、(25)、(26)中α
2表示第二旋轉關節(jié)與平移關節(jié)的夾角,θ
1表示第一旋轉
角度運動的角度;
由上面分析可以得到碰撞機率等于零的情況如下式所示:
第二種工作模式中η
1(η
1表示表示第二種工作模式下第一連桿與豎直面夾角)通過以下
公式得出:
公式(28)中h
3表示基座離人頭部的距離,h
4表示人的頭部高度,w
1表示病床的寬度,w
3表
示人胸腔的寬度;
第二種工作模式的極限情況下,不發(fā)生碰撞的條件為:
公式(28)中η
1表示第二種工作模式下第一連桿與豎直面夾角,
碰撞機率由下式得到:
CP=1-(1-N
1/M
1)*(1-N
2/M
2)......(29)
公式(29)中N
1表示第一種工作模式下碰撞機率為零次數,N
2表示第二種工作模式下碰
撞機率為零次數,M
1表示第一種工作模式下總的次數,M
2表示表示第二種工作模式下總的次
數。