Xem mẫu
- 84
Robot c«ng nghiÖp
ch−¬ng VII
§éng lùc häc Robot
(Dynamic of Robot)
7.1. NhiÖm vô vµ ph−¬ng ph¸p ph©n tÝch ®éng lùc häc robot
Nghiªn cøu ®éng lùc häc robot lµ c«ng viÖc cÇn thiÕt khi ph©n tÝch còng nh− tæng
hîp qu¸ tr×nh ®iÒu khiÓn chuyÓn ®éng. ViÖc nghiªn cøu ®éng lùc häc robot th−êng gi¶i
quyÕt hai nhiÖm vô sau ®©y :
1/ X¸c ®Þnh momen vµ lùc ®éng xuÊt hiÖn trong qu¸ tr×nh chuyÓn ®éng. Khi ®ã qui
luËt biÕn ®æi cña biÕn khíp qi(t) coi nh− ®· biÕt.
ViÖc tÝnh to¸n lùc trong c¬ cÊu tay m¸y lµ rÊt cÇn thiÕt ®Ó chän c«ng suÊt ®éng c¬,
kiÓm tra ®é bÒn, ®é cøng v÷ng, ®¶m b¶o ®é tin cËy cña robot.
2/ X¸c ®Þnh c¸c sai sè ®éng tøc lµ sai lÖch so víi qui luËt chuyÓn ®éng theo ch−¬ng
tr×nh. Lóc nÇy cÇn kh¶o s¸t Ph−¬ng tr×nh chuyÓn ®éng cña robot cã tÝnh ®Õn ®Æc tÝnh ®éng
lùc cña ®éng c¬ vµ c¸c kh©u.
Cã nhiÒu ph−¬ng ph¸p nghiªn cøu ®éng lùc häc robot, nh−ng th−êng gÆp h¬n c¶ lµ
ph−¬ng ph¸p c¬ häc Lagrange, cô thÓ lµ dïng ph−¬ng tr×nh Lagrange - Euler. §èi víi c¸c
kh©u khíp cña robot, víi c¸c nguån ®éng lùc vµ kªnh ®iÒu khiÓn riªng biÖt, kh«ng thÓ bá
qua c¸c hiÖu øng träng tr−êng (gravity effect), qu¸n tÝnh (initial), t−¬ng hæ (Coriolis), ly
t©m (centripetal)... mµ nh÷ng khÝa c¹nh nÇy ch−a ®−îc xÐt ®Çy ®ñ trong c¬ häc cæ ®iÓn; C¬
häc Lagrange nghiªn cøu c¸c vÊn ®Ò nªu trªn nh− mét hÖ thèng khÐp kÝn nªn ®©y lµ nguyªn
lý c¬ häc thÝch hîp ®èi víi c¸c bµi to¸n ®éng lùc häc robot.
7.2. C¬ häc Lagrange víi c¸c vÊn ®Ò ®éng lùc cña robot.
Hµm Lagrange cña mét hÖ thèng n¨ng l−îng ®−îc ®Þnh nghÜa :
L=K-P (7.1)
Trong ®ã : K lµ tæng ®éng n¨ng cña hÖ thèng
P lµ tæng thÕ n¨ng
K vµ P ®Òu lµ nh÷ng ®¹i l−îng v« h−íng nªn cã thÓ chän bÊt cø hÖ to¹ ®é thÝch hîp
nµo ®Ó bµi to¸n ®−îc ®¬n gi¶n. §èi víi mét robot cã n kh©u, ta cã :
n
n
P = ∑ Pi
K = ∑ Ki vµ
i =1
i =1
ë ®©y, Ki vµ Pi lµ ®éng n¨ng vµ thÕ n¨ng cña kh©u thø i xÐt trong hÖ to¹ ®é chän.Ta
biÕt mçi ®¹i l−îng Ki vµ Pi lµ mét hµm sè phô thuéc nhiÒu biÕn sè:
& &
Ki = K(qi, q i ) vµ Pi = P(qi, q i )
Víi qi lµ to¹ ®é suy réng cña khíp thø i. NÕu khíp thø i lµ khíp quay th× qi lµ gãc
quay θi, nÕu lµ khíp tÞnh tiÕn th× qi lµ ®é dµi tÞnh tiÕn di.
Ta ®Þnh nghÜa : Lùc t¸c dông lªn kh©u thø i (i=1, 2,..., n) víi quan niÖm lµ lùc tæng
qu¸t (Generalized forces), nã cã thÓ lµ mét lùc hoÆc mét momen (phô thuéc vµo biÕn khíp
qi lµ tÞnh tiÕn hoÆc quay), ®−îc x¸c ®Þnh bëi:
d ∂L ∂L
Fi = − (7.2)
dt ∂q i ∂q i
&
TS. Ph¹m §¨ng Ph−íc
- 85
Robot c«ng nghiÖp
Ph−¬ng tr×nh nÇy ®−îc gäi lµ ph−¬ng tr×nh Lagrange-Euler, hay th−êng ®−îc gäi t¾t
lµ ph−¬ng tr×nh Lagrange.
7.3. VÝ dô ¸p dông :
XÐt mét robot cã hai kh©u nh− h×nh vÏ, C¸c kh©u cã chiÒu dµi lµ d1 vµ d2 víi c¸c
khèi l−îng t−¬ng øng m1 vµ m2 qui ®æi vÒ ®Çu mót cña kh©u. Robot ®−îc ®Æt th¼ng ®øng
chÞu gia tèc träng tr−êng g. C¸c khíp chuyÓn ®éng quay víi c¸c biÕn khíp θ1 vµ θ2. TÝnh
lùc tæng qu¸t.
Qua vÝ dô nÇy, chØ víi mét mèi liªn kÕt hai
y kh©u, c¸c vÊn ®Ò ®Æt ra ®Òu ®· cã mÆt
g = 9,81m/s2 trong qu¸ tr×nh nghiªn cøu ®éng lùc häc,
vµ do ®ã, vÝ dô nªu trªn cã thÓ më réng ®Ó
¸p dông trong nh÷ng tr−êng hîp phøc t¹p
O0 x
x1 x2 h¬n. §èi víi kh©u 1 :
1 1 2&2
K 1 = m1 v 1 = m1 d 1 θ 1
2
(7.3)
z 2 2
θ1 P1 = -m1gd1cosθ1 (7.4)
m1
y1 §èi víi kh©u 2 :
VÒ to¹ ®é :
x2 = d1sinθ1 + d2sin(θ1 + θ2)
θ2
y2 m2
y2 = -d1cosθ1 - d2cos(θ1 + θ2)
ChiÒu cao thÕ n¨ng :
h = d1cosθ1 + d2cos(θ1 + θ2)
v2 = x2 + y2
2 2 2
& &
VÒ mÆt vËn tèc :
d & & &
x 2 = x 2 = d 1 cos(θ1 )θ 1 + d 2 cos(θ1 + θ 2 )(θ 1 + θ 2 )
&
Víi
dt
d & & &
y 2 = y 2 = d 1 sin(θ1 )θ1 + d 2 sin(θ1 + θ 2 )(θ1 + θ 2 )
&
dt
[ ]
2&2 &2 && &2 &2 & &
v 2 = d 1 θ1 + d 2 (θ 1 + 2θ1θ 2 + θ 2 ) + 2d 1d 2 cos(θ 2 )(θ 1 + θ 1θ 2 )
2 2
§éng n¨ng vµ thÕ n¨ng sÏ lµ :
[ ]
1 1 2&2 &2 && &2 &2 & &
K 2 = m 2 v 2 = m 2 d 1 θ1 + d 2 (θ 1 + 2θ1θ 2 + θ 2 ) + 2d 1d 2 cos(θ 2 )(θ 1 + θ 1θ 2 )
2
(7.5)
2
2 2
P2 = − m 2 g[d 1 cos(θ1 ) + d 2 cos(θ 1 + θ 2 )] (7.6)
7.4. Hµm Lagrange vµ lùc tæng qu¸t :
¸p dông hµm Lagrange cho vÝ dô trªn, ta cã :
L = (K1 + K2) - (P1 + P2)
1 1
2&2 &2 && & &2 & &
L = ( m1 + m 2 )d 1 θ1 + m 2 d 2 (θ1 + 2θ1θ 2 + θ 2 ) + m 2 d 1d 2 cos θ 2 (θ 1 + θ 1θ 2 ) +
2 2
2 2
+ ( m1 + m 2 )gd 1 cos θ1 + m 2 gd 2 cos(θ1 + θ 2 ) (7.7)
Khi tÝnh lùc tæng qu¸t, c¸c biÕn cña hÖ : q1 = θ1 vµ q2 = θ2.
§èi víi kh©u 1 :
∂L ∂L 2& & & & &
= = ( m1 + m 2 )d 1 θ1 + m 2 d 2 (θ 1 + θ 2 ) + 2 m 2 d 1d 2 cos θ 2 θ1 + m 2 d 1d 2 cos θ 2 θ 2
&
∂q 1 ∂θ1
2
&
TS. Ph¹m §¨ng Ph−íc
- 86
Robot c«ng nghiÖp
d ∂L
= ( m1 + m 2 )d 1 θ1 + m 2 d 2 (&&1 + && 2 ) − 2 m 2 d 1d 2 sin θ 2 θ 2 θ 1 + 2 m 2 d 1d 2 cos θ 2 θ 1 −
2 && && &&
θθ
2
&
dt ∂θ1
& &&
− m d d sin θ θ 2 + m d d cos θ θ
212 22 212 22
∂L ∂L
= = − ( m1 + m 2 )gd 1 sin θ1 − m 2 gd 2 sin(θ1 + θ 2 )
∂q 1 ∂θ1
VËy :
d ∂L ∂L &&
F1 = − = [( m1 + m 2 )d 1 + m 2 d 2 + 2 m 2 d 1d 2 cos θ 2 ]θ 1 +
2
&
dt ∂θ 1 ∂θ 1
2
&& && &2
+[ m 2 d 2 + m 2 d 1d 2 cos θ 2 ]θ 2 − 2 m 2 d 1d 2 sin θ 2 θ 2 θ 1 − m 2 d 1d 2 sin θ 2 θ 2 +
2
(7.8)
+ ( m1 + m 2 )gd 1 sin θ 1 + m 2 gd 2 sin(θ 1 + θ 2 )
Muèn cho kh©u 1 quay ®−îc mét gãc θ1 th× ®éng c¬ ph¶i t¹o ra mét lùc tæng qu¸t ≥
F1. Lùc tæng qu¸t nÇy cã ®Æc tÝnh phi tuyÕn, lµ hîp t¸c dông cña nhiÒu yÕu tè (non linear
and cuppling).
T−¬ng tù, ®Ó tÝnh lùc tæng qu¸t cña kh©u thø hai , ta cã :
∂L & & &
= m 2 d 2 θ1 + m 2 d 2 θ 2 + m 2 d 1d 2 cos θ 2 θ1
& 2 2
∂θ 2
d ∂L
= m 2 d 2 && 1 + m 2 d 2 θ 2 + m 2 d 1d 2 cos θ 2 θ1 − m 2 d 1d 2 sin θ 2 θ 1θ 2
&& && &&
2θ
& 2
dt ∂θ 2
∂L
= − m2 d1d 2 sin(θ 2 )θ&1θ&2 − − m2 d1d 2 sin(θ 2 )θ&12 − m2 gd 2 sin(θ1 + θ 2 )
vµ
∂θ 2
VËy :
d ∂L ∂L && 2& &
& − ∂θ = [m2 d 2 + m2 d1d 2 cosθ 2 ]θ1 + m2 d 2 θ 2
F2 = 2
dt ∂θ 2 (7.9)
2
− m d d sin(θ )θ& 2 + m gd sin(θ + θ )
2 1 2 2 1 2 2 1 2
§Ó ph©n tÝch ý nghÜa c¸c thµnh phÇn trong biÓu thøc tÝnh lùc tæng qu¸t, ta viÕt l¹i
c¸c biÓu thøc F1, F2 nh− sau :
F1 = D11&& 1 + D12 θ 2 + D111θ 1 + D122 θ 2 + D112 θ 1θ 2 + D121θ 1θ 2 + D1
&& &2 & && &&
θ 2
F2 = D12 &&1 + D 22 θ 2 + D 211θ 1 + D 222 θ 2 + D 212 θ1θ 2 + D 221θ 1θ 2 + D 2
&& &2 & && &&
θ 2
HiÖu øng HiÖu øng HiÖu øng HiÖu øng
qu¸n tÝnh ly t©m t−¬ng hæ träng tr−êng
Effective inertias Centripetal effect Coriolis effect Gravity
(Trong ®ã : D111 = 0; D222 = 0; D112 = D121 = D212 = D221 =-m2d1d2sinθ2 ...)
Trong c¸c biÓu thøc trªn, c¸c hÖ sè d¹ng Dii hoÆc D ij thÓ hiÖn hiÖu øng qu¸n tÝnh t¹i
khíp i hoÆc j g©y ra bëi gia tèc t¹i khíp i hoÆc j. C¸c sè h¹ng cã d¹ng D θ& 2 lµ lùc ly t©m ijj j
t¸c ®éng lªn khíp i g©y ra bëi vËn tèc t¹i khíp j. Sè h¹ng d¹ng D ijkθ& jθ&k + D ikjθ&kθ& j lµ lùc
Cariolis t¸c ®éng lªn khíp thø i g©y ra do vËn tèc t¹i khíp j vµ k. Sè h¹ng cã d¹ng Di lµ lùc
träng tr−êng t¸c ®éng lªn khíp i.
TS. Ph¹m §¨ng Ph−íc
- 87
Robot c«ng nghiÖp
7.5. Ph−¬ng tr×nh ®éng lùc häc robot :
XÐt kh©u thø i cña mét robot cã n kh©u. TÝnh lùc tæng qu¸t Fi cña kh©u thø i víi
khèi l−îng vi ph©n cña nã lµ dm. Lùc tæng qu¸t Fi ®ãng vai trß rÊt quan träng khi x©y dùng
s¬ ®å khèi ®Ó thiÕt lËp hµm ®iÒu khiÓn cho robot cã n bËc tù do.
7. 5. 1. VËn tèc cña mét ®iÓm trªn robot :
Mét ®iÓm trªn kh©u thø i ®−îc m« t¶ trong hÖ to¹ ®é c¬ b¶n lµ :
r = Ti. ir (7.10)
Trong ®ã : r lµ to¹ ®é cña ®iÓm xÐt ®èi víi kh©u thø i, ir kh«ng thay ®æi theo thêi
i
gian. Ti lµ ma trËn chuyÓn ®æi tõ kh©u thø i vÒ hÖ to¹ ®é gèc : Ti = A1A2...Ai. Nh− vËy r lµ
mét hµm cña thêi gian t.
Tèc ®é cña vi khèi l−îng dm ®−îc tÝnh bëi c«ng thøc :
⎛ i ∂T ⎞
dr d i
= Ti r = ⎜ ∑ i q j ⎟ i r
&= ⎜ j=1∂q & ⎟
r (7.11)
⎝ ⎠
dt dt j
Khi tÝnh b×nh ph−¬ng cña vËn tèc nÇy ta cã :
& . & = ∑ r 2 ( x o , y o , z o ) = Tr ( & & T )
&&&
rr rr (7.12)
y
Kh©u i
i
r
dm
Ti
r
x
O0
z
H×nh 7.1. Kh¶o s¸t tèc ®é cña vi khèi l−îng dm.
Víi rT lµ chuyÓn vÞ vect¬ vµ Tr lµ viÕt t¾t cña Trace (vÕt cña ma trËn) :
⎡ a 11 a 12 ... a 1n ⎤
⎢a a 22 ... a 2 n ⎥ n
Trace ⎢ ⎥ = ∑a
21
⎢ ... ... ... ... ⎥ i =1 ii
⎢ ⎥
⎣a n1 a n 2 a 11 a nn ⎦
Hay :
⎡x 2 ⎤
⎡ x⎤
⎢ ⎥
⎢ y ⎥ . [x y z ] = y2
⎢ ⎥
⎢⎥
⎢ z2 ⎥
⎢z⎥
⎣⎦ ⎣ ⎦
Do vËy
d d
& 2 = Tr ( &. & T ) = Tr ( Ti .i r. Ti T .i r T )
r rr
dt dt
TS. Ph¹m §¨ng Ph−íc
- 88
Robot c«ng nghiÖp
⎡ i ∂T ⎤
i ∂T
T
= Tr ⎢ ∑ i q j i r. ∑ i q k i r T ⎥
& &
⎣ j=1∂q j k =1 ∂q k
⎢ ⎥
⎦
⎡ i i ∂Ti i i T ∂Ti ⎤
T
= Tr ⎢∑∑ &&
q j qk ⎥
rr . (7.13)
⎢ j =1 k =1 ∂q j ∂qk ⎥
⎣ ⎦
7. 5. 2. TÝnh ®éng n¨ng cña vi khèi l−îng dm.
Ký hiÖu Ki lµ ®éng n¨ng cña kh©u thø i. dKi lµ ®éng n¨ng cña vi khèi l−îng dm ®Æt
t¹i vÞ trÝ ir trªn kh©u thø i.
1 ⎡ i i ∂Ti i i T ∂Ti ⎤
T
dK i = Tr ⎢ ∑ ∑ q j q k ⎥ dm
&&
rr .
2 ⎣ j=1k =1 ∂q j ∂q k
⎢ ⎥
⎦
1 ⎡ i i ∂Ti i ⎤
∂T T
= Tr ⎢ ∑ ∑ ( r. dm.i r T ). i q jq k ⎥
&& (7.14)
2 ⎢ j=1k =1 ∂q j ∂q k ⎥
⎣ ⎦
Vµ do ®ã ®éng n¨ng cña kh©u thø i sÏ lµ :
1 ⎡ i i ∂Ti ⎤
∂Ti T
K i = ∫ dK = Tr ⎢∑ ∑ ( ∫ r. r dm ). &&
i iT
q j qk ⎥ (7.15)
2 ⎢ j =1 k =1 ∂q j Khau i ∂q k ⎥
⎣ ⎦
Khau i
∫
Ji = i
r.i r T dm gäi lµ ma trËn gi¶ qu¸n tÝnh (Pseudo inertia matrix).
§Æt
Khau i
ý nghÜa "gi¶ qu¸n tÝnh" ®−îc sö dông v× khi thiÕt lËp ®Çy ®ñ c¸c phÇn tö cña ma trËn Ji ta
cã thÓ liªn hÖ víi c¸c kh¸i niÖm "m«men qu¸n tÝnh ®éc cùc" vµ tr×nh bµy c¸c phÇn tö cña Ji
gièng nh− c¸c phÇn tö cña m«men qu¸n tÝnh ®éc cùc. Ta xÐt mèi quan hÖ nÇy nh− sau :
Theo ®Þnh nghÜa ta cã :
⎡ i x 2 dm
∫ x ydm ∫ x zdm ∫ xdm ⎤
∫
ii ii i
⎢ii ⎥
⎢ ∫ x ydm ∫ y dm ∫ y zdm ∫ ydm ⎥
i2 ii i
J i = ∫ i r.i r T dm = Ji = ⎢ i i ⎥ (7.16)
⎢∫ x zdm ∫ i y i zdm ∫ i z 2 dm ∫ i zdm ⎥
Khau i
⎢ i xdm
∫ ydm ∫ zdm ∫ dm ⎥
⎣∫
i i
⎦
B©y giê ta nh¾c l¹i m«men qu¸n tÝnh ®éc cùc cña
mét vËt thÓ bÊt kú nh− h×nh vÏ.
yx
ω Theo ®Þnh nghÜa ta cã :
I xx = ∫ ( y 2 + z 2 )dm
z
I yy = ∫ x 2 + z 2 )dm
I zz = ∫ ( x 2 + y 2 )dm
H×nh 7.2 : M«men qu¸n tÝnh ®éc cùc
1 1 1
x2 = − ( y2 + z2 ) + (x2 + z2 ) + (x2 + y2 )
Vµ v× :
2 2 2
∫ x dm =( −I xx + I yy + I zz ) / 2 ; .v.v…
2
VËy :
Ngoµi ra ta cßn cã :
I xy = ∫ xydm ; I yz = ∫ yzdm ; I xz = ∫ xzdm
mx = ∫ xdm ; my = ∫ ydm ; mz = ∫ zdm
TS. Ph¹m §¨ng Ph−íc
- 89
Robot c«ng nghiÖp
§èi chiÕu víi ma trËn gi¶ qu¸n tÝnh Ji, ta cã thÓ tr×nh bµy Ji nh− sau :
⎡ − I xx + I yy + I zz ⎤
⎢ mx ⎥
I yx I zx
2
⎢ ⎥
I xx − I yy + I zz
⎢ my⎥
I xy I zy
ji = ⎢ ⎥ (7.17)
2
I xx + I yy − I zz
⎢ ⎥
⎢ mz ⎥
I yz I yz
2
⎢ m⎥
⎣ ⎦
mx my mz
Nh− vËy ý nghÜa biÓu tr−ng cña Ji ®· râ.
1 ⎡ i i ∂Ti ∂Ti ⎤
T
K i = Tr ⎢∑∑ &&
q j qk ⎥
VËy ta cã : Ji (7.18)
2 ⎢ j =1 k =1 ∂q j ∂q k ⎥
⎣ ⎦
Cuèi cïng, §éng n¨ng cña mét robot cã n kh©u ®−îc tÝnh :
n
K = ∑ Ki (7.19)
i =1
7. 5. 3. TÝnh thÕ n¨ng cña robot :
ThÕ n¨ng cña kh©u i cã khèi l−îng mi, träng t©m ®−îc x¸c ®Þnh bëi vect¬ ri (vect¬
biÓu diÔn träng t©m cña kh©u i trong hÖ to¹ ®é c¬ b¶n) lµ :
Pi = -mi. g. ri = -mi. g. Ti iri (7.20)
Trong ®ã, vect¬ gia tèc träng tr−êng g ®−îc biÓu diÔn d−íi d¹ng mét ma trËn cét :
⎡gx ⎤ ⎡ 0 ⎤
⎢g ⎥ ⎢ 0 ⎥
g = ⎢ y⎥ = ⎢ ⎥
⎢ g z ⎥ ⎢ − 9,8⎥
⎢⎥⎢ ⎥
⎣0⎦ ⎣ 0 ⎦
ThÕ n¨ng cña toµn c¬ cÊu robot n kh©u ®éng sÏ lµ :
n
P = − ∑ mi gTi i ri (7.21)
i =1
7. 5.4. Hµm Lagrange :
Sau khi x¸c ®Þnh ®éng n¨ng vµ thÕ n¨ng cña toµn c¬ cÊu, ta cã hµm Lagrange cña
robot cã n bËc tù do :
⎛ ∂T ⎞
∂T T
1n i i n
⎟q j qk + 1 ∑ mi gTi i ri
∑∑∑ Trace⎜ i J i i
L= && (7.22)
⎜ ∂q ⎟
∂q k
2 i =1 j =1 k =1 2 i =1
⎝j ⎠
Chóng ta chó ý r»ng, trong hµm Lagrange vÉn ch−a ®Ò cËp ®Õn ¶nh h−ëng cña
nguån truyÒn ®éng (gåm c¸c phÇn tÜnh (stator) vµ phÇn ®éng (Rotor) cña ®éng c¬ ®iÖn).
7. 5. 5. Ph−¬ng tr×nh ®éng lùc häc robot :
Ta ®· biÕt lùc tæng qu¸t ®Æt lªn kh©u thø i cña robot cã n kh©u (Ph−¬ng tr×nh
Lagrange - Euler) :
d ∂L ∂L
Fi = − (7.23)
dt ∂q i ∂q i
&
Sau khi thiÕt lËp hµm Lagrange, víi p = 1... n, ta tÝnh ®−îc :
TS. Ph¹m §¨ng Ph−íc
- 90
Robot c«ng nghiÖp
(p lµ chØ sè lÇn l−ît lÊy theo j vµ k)
⎛ ∂Ti ∂Ti T ⎞ ⎛ T⎞
⎟qk + 1 ∑∑ Tr⎜ ∂Ti J i ∂Ti ⎟q j
∂L 1 n i n i
= ∑∑ Tr ⎜ & &
Ji (7.24)
∂q p 2 i =1 k =1 ⎜ ∂q p ∂qk ⎟ 2 i =1 j =1 ⎜ ∂q j ∂q p ⎟
& ⎝ ⎠ ⎝ ⎠
Thay ®æi chØ sè gi¶ j thµnh k trong sè h¹ng thø hai ,vµ ®Ó ý r»ng :
T
⎛ ∂Ti ∂Ti T ⎞ ⎛ T⎞ ⎛ T⎞
⎟ = Tr ⎜ ∂Ti J i ∂Ti ⎟ = Tr ⎜ ∂Ti J i ∂Ti ⎟
⎜
Tr J (7.25)
⎜ ∂q i ∂q ⎟ ⎜ ∂q ∂q p ⎟ ⎜ ∂q ∂q j ⎟
⎝j p⎠ ⎝j ⎠ ⎝p ⎠
⎛ ∂Ti ∂Ti T ⎞
∂L n i
= ∑∑ Tr ⎜ ⎟qk
&
Ji
ta cã : (7.26)
∂q p i =1 k =1 ⎜ ∂qk ∂q p ⎟
& ⎝ ⎠
Còng ®Ó ý r»ng : trong Ti(q1, q2, . . . , qi), víi qi lµ c¸c biÕn khíp cña i khíp ®Çu tiªn. Do
∂Ti
= 0.
vËy, nÕu i < p th×
∂q p
⎛ ∂T ⎞
∂T T
∂L n i
= ∑∑ Tr ⎜ i J i i ⎟qk
&
Cuèi cïng ta cã : (7.27)
∂q p i = p k =1 ⎜ ∂qk ⎟
∂q p
& ⎝ ⎠
LÊy vi ph©n theo thêi gian t cña ph−¬ng tr×nh trªn :
⎛ ∂T ⎞
∂T T
d ∂L dni
= ∑∑ Tr⎜ i J i i ⎟ qk
&
dt ∂q p dt i = p k =1 ⎜ ∂qk ⎟
∂q p
& ⎝ ⎠
⎛ ∂T ⎞ ⎡2 T⎤
∂T T ⎟ qk + ∑∑∑ Tr ⎢ ∂ Ti J i ∂Ti ⎥ qk qm +
n i n i i
= ∑∑ Tr ⎜ i J i i && &&
⎜ ∂q ⎟
∂q p ⎢ ∂qk ∂qm ∂q p ⎥
⎝k ⎠ ⎣ ⎦
i = p k =1 i = p k =1 m =1
⎡ ∂ 2Ti ∂T T ⎤
n i i
+ ∑∑∑ Tr ⎢ &&
J i i ⎥ qk qm (7.28)
⎢ ∂q p ∂qm ∂qk ⎥
⎣ ⎦
i = p k =1 m =1
(BiÕn ®æi theo chó ý (7.25))
Sè h¹ng cuèi cña ph−¬ng tr×nh Lagrange Euler lµ :
⎛ ∂ 2Ti ⎞
∂T T
∂L 1 n i i
= ∑∑∑ Tr⎜ ⎟ q j qk +
&&
Ji i
∂q p 2 i = p j =1 k =1 ⎜ ∂q j ∂q p ⎟
∂qk
⎝ ⎠
⎛ ∂ 2Ti ⎞
∂T T ⎟ q j qk + ∑ mi g ∂Ti i ri
1n i i n
∑∑∑ Tr ⎜
+ &&
Ji i (7.29)
2 i =1 j =1 k =1 ⎜ ∂qk ∂q p ⎟
∂q j ∂q p
⎝ ⎠ i= p
d ∂L ∂L
Fp = −
Cuèi cïng ta cã lùc tæng qu¸t cña kh©u p :
dt ∂q p ∂q p
&
Thay thÕ c¸c chØ sè p vµ i thµnh i vµ j, ta sÏ cã :
⎛ ∂T j ⎞ ⎡ ∂ 2T j ∂T jT ⎤
∂T jT ∂T j j
j j j
n n n
Fi = ∑∑ Tr ⎜ ⎟ qk + ∑∑∑ Tr ⎢ qk qm − ∑ m j g
&& &&
⎥
Jj Jj rj
⎜ ∂q ⎟
∂qi ∂qk ∂qm ∂qi ⎦ ∂qi
⎝k ⎠ ⎣
j =i k =1 j =i k =1 m =1 j =i
(7.30)
Víi mét robot cã n bËc tù do th× :
TS. Ph¹m §¨ng Ph−íc
- 91
Robot c«ng nghiÖp
q = [q1, q2, . . . ,qn]T
q = [q 1 , q 2 , ... , q n ] T
& && &
F = F[F1, F2, . . . , Fn]T
vµ
§Ó cho gän, ta biÓu diÔn :
F = J ( q ) q + C ( q, q ) q + G ( q )
&& && (7.31)
Trong ®ã :
J thÓ hiÖn t¸c dông cña qu¸n tÝnh, lµ mét ma trËn ®èi xøng (n x n);
C thÓ hiÖn t¸c dông cña lùc ly t©m vµ Cariolis, lµ mét vect¬ (n x 1);
G thÓ hiÖn t¸c dông cña lùc träng tr−êng, còng lµ mét vect¬ (n x 1).
§©y lµ ph−¬ng tr×nh ®éng lùc häc cña robot.
NÕu thªm vµo ph−¬ng tr×nh trªn c¸c t¸c dông kh¸c nh− : FEX ®Æc tr−ng cho c¸c
ngo¹i lùc t¸c dông lªn trôc, V ®Æc tr−ng cho hiÖu øng ma s¸t, ta cã :
F = J ( q) q + C ( q, q) q + G ( q) + V ( q) + FEX
&& && & (7.32)
TS. Ph¹m §¨ng Ph−íc
nguon tai.lieu . vn
