Iterative General Dynamic Model for Serial-Link Manipulators

Size: px
Start display at page:

Download "Iterative General Dynamic Model for Serial-Link Manipulators"

Transcription

1 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general dynamc model for seral-lnk manpulators. Ths dynamc model wll relate the set of torques or forces τ requred at the jonts of the manpulator (torques for revolute jonts, forces for prsmatc jonts) to acheve a partcular set of jont postons, veloctes and acceleratons ( ΘΘ,, Θ ): τ Note that n equaton (1), τ, Θ, Θ and Θ are all n 1 vectors, where n s the number of jonts n the manpulator, and h s some nonlnear functon. In order to derve such a dynamc model, we wll do the followng: 1. Relate lnear and angular acceleratons of coordnate frames wth respect to one another. 2. Generalze the basc laws of motons to three dmensons and apply those laws of moton to the seral confguraton of manpulators. 3. Extend the concept of moments of nerta to three dmensons (nerta tensor). [Ths part of the dscusson wll not be part of ths set of notes, but wll be handled elsewhere. 2. cceleratons between coordnate frames. asc defntons The lnear acceleraton Q of a pont Q wth respect to some coordnate frame { } s defned as the tme dervatve of the lnear velocty V Q : Note that ths defntons s smlar to that of lnear velocty tself (from Chapter 5): Smlarly, the angular acceleraton Ω of coordnate frame { } wth respect to coordnate frame { } s defned as the tme dervatve of the angular acceleraton Ω :. Notaton h( ΘΘ,, Θ ) d V Q Q ( V Q ) lm ( t+ t) V Q () t t t d V Q ( Q) lm Q( t+ t) Q() t t t d Ω ( Ω ) lm ( t+ t) Ω () t t t We defne the followng short-hand notaton for the lnear and angular acceleratons of some coordnate frame { } wth respect to a fxed unversal reference frame { } : (1) (2) (3) (4) v Ȧ ω ORG (smlar to v VORG from Chapter 5) (5) Ω (smlar to ω Ω from Chapter 5) (6) C. Lnear acceleratons between coordnate frames In Chapter 5, we derved the followng mportant equaton for two coordnate frames { } and { } wth concdent orgns (.e. no translaton between { } and { } ): VQ R( VQ) + [ R( Q) ({ } and { } orgns concdent). (7) - 1 -

2 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators Let us rewrte equaton (7) to establsh an mportant formula for the tme dervatve of a rotaton matrx tmes a vector: VQ R ( VQ) + [ R( Q) d [ R ( Q) R( Q ) + [ R( Q) (8) (9) d [ R ( Q) R( Q ) + [ R( Q) Note that equaton (10) gves us a general formula for the tme dervatve of RQ where R s some rotaton matrx and Q s some vector. Now, let us dfferentate equaton (8) wth respect to tme: (10) VQ R ( VQ) + [ R( Q) d Q R [ ( VQ) Ω [ R( Q) d + + Ω [ R ( Q) In equaton (12) we made use of the followng dentty for any 3-space vectors Q and P : (11) (12) d dq dp ( Q P) P + Q Q P + Q P. (13) Let us now substtute equaton (10) nto equaton (12): d Q R [ ( VQ) Ω [ R( Q) d + + Ω [ R ( Q) Q R( Q ) R + [ ( VQ) + Ω [ R( Q) + Ω R( Q ) R + [ ( Q) In equaton (15), let Q VQ so that: Q R( Q ) R + [ ( VQ) + Ω [ R( Q) + Ω R( VQ) R( + [ Q) (14) (15) (16) Note that we can combne two cross-product terms n equaton (16), Q R( Q ) R + [ ( VQ) + Ω [ R( Q) + Ω R( VQ) R( + [ Q) Q R( Q ) 2 R + [ ( VQ) + Ω [ R( Q) + [ Ω R( Q) (17) (18) - 2 -

3 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators Q R( Q ) 2 R + [ ( VQ) + Ω [ R( Q) + [ R( Q) (19) Equaton (19) gves the angular acceleraton Q for a vector Q defned wth respect to coordnate frame { }, when the orgns of coordnate frames { } and { } are concdent wth one another. If the orgn of coordnate frame { } s acceleratng wth respect to { } equaton (19) s easly modfed to nclude that adonal lnear acceleraton: Q ORG R( Q ) 2 R + + [ ( VQ) + Ω [ R( Q) + [ Ω R( Q) (20) Q ORG R( Q ) 2 R + + [ ( VQ) + Ω [ R( Q) + [ R( Q) (21) Equaton (21) allows for the possblty that vector Q s movng wth respect to coordnate frame { }. Let us now consder a more restrctve case namely, that vector Q s fxed wth respect to coordnate frame { }. In other words, we wll assume that there s movement only between coordnate frames { } and { } (and not wthn coordnate frame { } ), so that: VQ 0 Q 0 (22) (23) Ths assumpton smplfes equaton (21) substantally: Q ORG R( Q ) 2 R + + [ ( VQ) + Ω [ R( Q) + [ Ω R( Q) (24) Q ORG + Ω [ R( Q) + [ R( Q), VQ Q 0 (25) D. ngular acceleratons between coordnate frames Now, let us consder angular acceleratons between dfferent coordnate frames. Let us begn wth the relatonshp of angular veloctes between three coordnate frames { }, { } and { C} : ΩC + R( ΩC) Dfferentatng (26) and agan substtutng equaton (10), we get: Ω Ω d C + [ R ( ΩC) Ω Ω C + R( Ω C ) + [ R( ΩC) Ω Ω C + R( Ω C ) + [ R( ΩC) (26) (27) (28) (29) E. Summary of results elow, we summarze our results on lnear and angular acceleratons: Q ORG R( Q ) 2 R + + [ ( VQ) + Ω [ R( Q) + [ Ω R( Q) (30) Q ORG + Ω [ R( Q) + [ Ω R( Q), VQ Q 0 (31) Ω Ω C + R( Ω C ) + [ R( ΩC) (32) - 3 -

4 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 3. asc equatons of moton. Newton s law Consder Fgure 1 below, whch depcts a rgd body, whose center of mass s acceleratng wth acceleraton v Ċ under a net force F actng on the body. (Crag, Fg. 6.3) Fgure 1: Rgd body, whose center of mass s acceleratng under the acton of a net force F. Newton s second law of moton relates F and v Ċ : F mv Ċ (33) where, F m v Ċ net force actng on the body, (34) mass of the body, and, (35) acceleraton of the center of mass of the body. (36). Euler s law Consder Fgure 2 below, whch depcts a rgd body, whch s rotatng wth angular velocty ω and angular acceleraton ω under a net moment N actng on the body. (Crag, Fg. 6.4) Fgure 2: Rgd body, whch s rotatng under the acton of a net moment N

5 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators Euler s law of rotatonal moton for rgd bodes relates N and ( ωω, ): N where, C I ω ω C + I ω (37) N net moment actng on the body, (38) C I ω ω 3 3 nerta tensor, wrtten wth respect to coordnate frame { C} (at center of mass), (39) angular velocty of the body, and, (40) angular acceleraton of the body. (41) C. Dynamc modelng Gven the results of Secton 2 on acceleratons, and the basc laws of moton n the prevous two sub-sectons, we wll now proceed as follows n dervng the dynamc model for a seral-lnk manpulator: τ h( ΘΘ,, Θ ) We wll assume that the jont postons Θ, jont veloctes Θ and jont acceleratons Θ are known, so that we can compute the correspondng jont torques/forces τ requred to acheve the known jont trajectory. We then break down the development of the dynamc model n equaton (42) nto three man tasks: 1. Outward propagaton of veloctes and acceleraton from the base coordnate frame { 0} to the end-effector coordnate frame { N}. 2. Newton s and Euler s equatons of moton from the base coordnate frame { 0} to the end-effector coordnate frame { N}. 3. Inward propagaton of force balance and moment balance equatons from coordnate frame { N} to coordnate frame { 1}. 4. Propagaton of veloctes and acceleratons. ngular veloctes and acceleratons In Chapter 5, we developed the followng relatonshp for angular veloctes between consecutve lnks: ω + 1 R + 1 ω + θ + 1 Ẑ + 1 (42) (velocty propagaton) (43) Now, we want to develop an analogous relatonshp for angular acceleratons: + 1 ω + 1 g ( ω, ω ) where g () represents some functonal mappng. To do ths, we wll apply the general relatonshp that we developed n Secton 2 for the propagaton of angular acceleratons: (44) ΩC + R( ΩC) (45) Ω Ω C + R( Ω C ) + [ R( ΩC). (46) + 1 Frst, we wll get equaton (43) nto the same form as equaton (45). Recall that ω and + 1 ω are shorthand notaton, and can be wrtten less compactly as, ω lso note that, R ( Ω) and ω + 1 ( Ω + 1). (47) - 5 -

6 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators + 1 θ + 1 Z +ˆ R( Ω + 1). (48) Equaton (48) requres some explanaton. The rght-hand sde of (48) denotes the angular velocty of coordnate frame { + 1} wth respect to {} expressed n the { + 1} coordnate frame. Now, thnk about the left-hand sde of (48). The angular velocty between frames {} and { + 1} s gven exactly by the jont rate θ + 1 orented along the Ẑ + 1 axs, and that the left-hand sde of (48) s expressed n terms of the { + 1} coordnate frame. Substtutng (47) and (48) nto equaton (43), ω + 1 R + 1 ω + θ + 1 Ẑ ( Ω + 1) R + 1 ( Ω) + R( Ω + 1) ( Ω + 1) R + 1 ( Ω) + R( Ω + 1) Now, let and + 1 C so that equaton (51) becomes: C C ( ΩC) C R ( Ω) + C R( ΩC) ( ΩC) C ( Ω) + C R( ΩC) Multplyng equaton (53) by CR and lettng, CR C ( ΩC) CR C ( Ω) + CR C R( ΩC) (49) (50) (51) (52) (53) (54) ΩC Ω + R( ΩC) (55) ΩC Note that we have now transformed equaton (43) nto (45), whch means that equaton (46) can be used to propagate angular acceleratons from one lnk to the next for seral-lnk manpulators. In summary, + R( ΩC) (56) ω + 1 R + 1 ω + θ + 1 Ẑ + 1 C Ω + R( ΩC) wth substtutons: ω + 1 ( Ω + 1) ω ( Ω) (57) (58) (59) + 1 θ + 1 Z +ˆ R( Ω + 1), (60), + 1 C and. (61) From equaton (46), angular acceleratons are propagated by, Ω Ω C + R( Ω C ) + [ R( ΩC) Let us now transform ths relatonshp nto lnk notaton by reverse substtuton. Frst let, C + 1 and, so that: Ω + 1 Ω + R( Ω + 1) + Ω [ R( Ω + 1). (63) (62) - 6 -

7 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators + 1 Multply equaton (63) by : Ω + 1 Ω + 1 R( Ω + 1) Ω R + + [ ( Ω + 1) (64) In order to further develop equaton (64), we wll need the followng dentty: for any 3-space vectors Q, P and rotaton matrx R, RQ ( P) ( RQ) ( RP). (65) Thus, we can modfy the last term of (64) usng (65): Ω + 1 Ω R + 1 R( Ω + 1) Ω R + + [ ( Ω + 1) Ω + 1 Ω R + 1 R( Ω + 1) R [ Ω [ R( Ω + 1) sng, (66) (67) + 1 R R and R, (68) equaton (67) further smplfes to: Ω + 1 Ω R Ω + 1 R [ Ω ( R Ω + 1) Ω + 1 R ( Ω ) R Ω + 1 R R [ ( Ω) ( R Ω + 1) Fnally, we make the followng substtutons nto equaton (70): (69) (70) Ω ω (by defnton) (71) Ω ω and ( Ω + 1 ) ω + 1 (by defnton) (72) + 1 R + 1 Ω + 1 θ + 1 Ẑ R Ω + 1 θ + 1 Ẑ + 1 These substtutons result n: [prevously explaned n dscusson followng (48) (73) [analogous to (73) above (74) Ω + 1 R ( Ω ) R Ω + 1 R [ ( Ω) ( R Ω + 1) ω + 1 R ω θ + 1 Ẑ + 1 R ( ω) ( θ + 1 Ẑ + 1) ω + 1 R ω R ( ω) θ + 1 Ẑ θ + 1 Ẑ + 1 (75) (76) (77). ngular veloctes and acceleratons: summary Summarzng the results developed n the prevous secton, angular veloctes and acceleratons are propagated from lnk to lnk + 1 usng the followng two equatons (for revolute jonts): ω + 1 R + 1 ω + θ + 1 Ẑ + 1 [from (43) (78) - 7 -

8 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators ω + 1 R ω R ( ω) θ + 1 Ẑ θ + 1 Ẑ + 1 [from (77) (79) + 1 Note that equatons (78) and (79) allow us to compute the angular velocty ω + 1 of lnk frame { + 1} + 1 n terms of the angular velocty ω of lnk frame {} and θ + 1 ; and the angular acceleraton ω + 1 of lnk frame { + 1} n terms of the angular acceleraton ω and angular velocty ω of lnk frame {} and θ 1. lso, note that for a prsmatc jont, + θ 1 + θ so that (78) and (79) smplfy (for prsmatc jonts) to: ω + 1 R ω ω + 1 R ω (80) (81) (82) C. Lnear acceleratons of lnk frame orgns In Secton 2, we derved the followng relatonshp for the propagaton of lnear acceleratons: Q ORG R( Q ) 2 R + + [ ( VQ) + Ω [ R( Q) + [ Ω R( Q) We wll now convert equaton (83) nto lnk-specfc form as we dd above for angular acceleratons. Frst, let us make the followng substtutons: (83) (84) (85) Q + 1, ORG (for subscrpts) (86) Gven these substtutons, note that Q now denotes the orgn of lnk frame { + 1} n terms of lnk frame {} ; n short, Q P + 1, ORG Thus, equaton (83) becomes, (87) + 1, ORG, ORG + R ( + 1, ORG ) + 2 Ω [ R ( V + 1, ORG ) + Ω [ R ( P + 1, ORG ) + Ω [ Ω R ( P + 1, ORG ) + 1 Let us now pre-multply equaton (88) by : , ORG, ORG R ( + 1, ORG ) Ω [ R ( V + 1, ORG ) Ω [ R ( P + 1, ORG ) Ω + 1 [ Ω R ( P + 1, ORG ) + 1 Note that we used vector dentty (65) n dstrbutng over the cross product of vectors n (89). Let us now consder each of the terms n equaton (89) one by one. The left-hand sde of (89) can be wrtten n short-hand notaton as: (88) (89) , ORG + 1 (by defnton) (90) - 8 -

9 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators Smlarly for the frst term on the rght-hand sde of equaton (89): ORG, R ( ORG, ) + 1 R Equatons (90) and (91) smplfy equaton (89) to: (by defnton) (91) R R ( + 1, ORG ) Ω [ R ( V + 1, ORG ) Ω [ R ( P + 1, ORG ) Ω + 1 [ Ω R ( P + 1, ORG ) Let us now consder the second term on the rght-hand sde of equaton (92). The notaton + 1, ORG ndcates the lnear acceleraton of the orgn of lnk frame { + 1} wth respect to (and n terms of) lnk frame {}. Thus, for a seral-lnk manpulator, we can wrte: (92) + 1 R R ( + 1, ORG ) + 1 R( + 1, ORG ) ḋ Ẑ + 1 Smlarly, (93) + 1 R ( V + 1, ORG ) + 1 R( V + 1, ORG ) d Ẑ + 1 Note that lnear relatonshps n (93) and (94) are analogous to the manpulator-specfc angular relatonshps n (73) and (74). Equaton (92) now smplfes to: (94) R + ḋ + 1 Ẑ Ω d + 1 Ẑ Ω [ R ( P + 1, ORG ) Ω + 1 [ Ω R ( P + 1, ORG ) Let us now consder the angular velocty and acceleraton terms n (95). From pror dscusson, Ω R + 1 ( Ω) R ω + 1 Ω + 1 R + 1 ( Ω ) R ω so that equaton (95) further reduces to: R + ḋ + 1 Ẑ R + 1 ( ω) d + 1 Ẑ R ω [ R R ( P + 1, ORG ) R + 1 ω R + 1 [ ω R ( P + 1, ORG ) Next, let us make the followng substtuton n (98): + 1 R ( P + 1, ORG ) + 1 R( P + 1, ORG ) + 1 R ( P + 1) (95) (96) (97) (98) (99) - 9 -

10 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators where P s smply short-hand notaton for + 1 P + 1, ORG. Thus, equaton (98) reduces to: R + 1 ḋ + 1 Ẑ R ( ω) d + 1 Ẑ Fnally, note that the relatonshp n (65) can be rewrtten as, so that we can rearrange and group terms n equaton (100): Thus (for a prsmatc jonts): R ω R + 1 [ ( P + 1) R + 1 ω R [ ω R( P + 1) ( RQ) ( RP) RQ ( P) R + 1 ḋ + 1 Ẑ R ( ω) d + 1 Ẑ R ω R + 1 [ ( P + 1) R + 1 ω R [ ω R( P + 1) R + 1 ḋ + 1 Ẑ R ( ω) d + 1 Ẑ R( ω P + 1) + R[ ω ( ω P + 1) (100) (101) (102) (103) R ω P + 1 ω ω P [ + ( + 1) R ( ω) d + 1 Ẑ ḋ + 1 Ẑ Note that equaton (104) allow us to compute the lnear acceleraton + 1 of lnk frame { + 1} n terms of the angular acceleraton ω, lnear acceleraton and angular velocty ω of lnk frame {}, and d + 1 and ḋ + 1. lso, note that for a revolute jont, (104) ḋ + 1 d so that (104) smplfes (for revolute jonts) to: (105) R ω P. (106) + 1 ω ω P [ + ( + 1) + D. Lnear acceleraton of a lnk s center of mass In order to apply Newton s second law of moton n equaton (33), we need to know not just the lnear acceleraton of lnk frame {} but of the center of mass C of lnk as well. Once agan, we wll begn wth the relatonshp for the propagaton of lnear acceleratons derved n Secton 2: Q ORG R( Q ) 2 R + + [ ( VQ) + Ω [ R( Q) + [ Ω R( Q) (107) Let us make the followng substtutons: (108) (109) Q C (for subscrpts) (110) Gven these substtutons, note that Q now denotes the orgn of the center of mass of lnk frame { + 1} n terms of lnk frame {} ; n short, Q PC (111)

11 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators Thus, equaton (107) becomes, C, ORG + R( ) + 2 Ω [ R( V ) + C Ω [ R( P ) + Ω [ Ω R( P ) C C C (112) In equaton (112), note that for rgd lnks, C VC 0 (113) so that (112) reduces from, C, ORG + R( ) + 2 Ω [ R( V ) + C Ω [ R( P ) + Ω [ Ω R( P ) C C C (114) to the smplfed form, C, ORG + Ω [ R( P ) + Ω [ Ω R( P ) C C (115) Let us now pre-multply equaton (115) by : C, ORG + Ω [ R( P ) + Ω [ Ω R( P ) C C (116) Note that we used vector dentty (65) n dstrbutng over the cross product of vectors n (116). Smlar to earler dervaton, we now make the followng substtutons: C C (by defnton) (117) ORG, (by defnton) (118) R Ω ω and R Ω ω (by defnton) (119) Equaton (116) consequently reduces to: C + ω [ R R( P ) + ω [ ω R( P ) C C (120) Notng that, R( P ) P C C (121) equaton (120) now reduces to: C C + ω P + ω [ ω P C C + ω P + ω [ ω P C C (122) (123)

12 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 5. Lnk-specfc equatons of moton. Manpulator-specfc equatons of motons We can rewrte the basc equatons of moton n (33) and (37) n more lnk specfc notaton. Specfcally, F m C (124) N C I ω + ω C I ω (125) where, F m C N net force actng on lnk, (126) mass of lnk, (127) acceleraton of center of mass of lnk, (128) net moment actng on lnk, (129) C I nerta tensor, wrtten n frame { C } located at the center of mass of lnk, (130) ω ω angular velocty of lnk, and, (131) angular acceleraton of lnk. (132) We can now completely and succnctly wrte the outward propagaton of veloctes and net forces/moments. In the followng subsectons, we do so for revolute and prsmatc jonts, respectvely.. Summary of outward teraton 1. Revolute jonts: ω + 1 R + 1 ω + θ + 1 Ẑ ω + 1 R ω R ( ω) θ + 1 Ẑ θ + 1 Ẑ R ω P + 1 ω ω P [ + ( + 1) + 2. Prsmatc jonts: ω + 1 R ω ω + 1 R ω (133) (134) (135) (136) (137) R ω P + 1 ω ω P [ + ( + 1) R ( ω) d + 1 Ẑ ḋ + 1 Ẑ + 1 (138) 3. oth jont types: + 1 C ω PC ω [ ω + 1 PC + 1 (139)

13 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators F + 1 m + 1 v C + 1 C N I + 1 ω + 1 ω I + 1 ω Force/moment balance equatons C 1 (140) (141). Introducton Equatons (133) through (141) gve us the net forces/moments at a gven lnk requred to cause the desred/ known jont moton ( ΘΘ,, Θ ). We now must fgure out what part of those net forces/moments must be suppled by the jont actuators. To do ths, we wll wrte force/moment balance equatons about the center of mass of each lnk. Consder Fgure 3 below, whch llustrates the forces and moments actng on lnk. In Fgure 3 we use the followng notaton: P + 1 V 2 V 1 Fgure 3: Forces and moments acton on lnk. (Crag, Fg. 6.5) f force exerted on lnk by lnk 1, (142) f + 1 force exerted on lnk + 1 by lnk, (143) n moment exerted on lnk by lnk 1, (144) n + 1 moment exerted on lnk + 1 by lnk, (145) and, as before, F N net force actng on lnk, and, (146) net moment actng on lnk. (147) lso, P + 1 vector from the orgn of frame {} to the orgn of frame { + 1}, (148) V 1 vector from the center of mass of lnk to the orgn of coordnate frame {}, and, (149)

14 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators V 2 vector from the center of mass of lnk to the orgn of coordnate frame { + 1}. (150). Force balance equaton Gven the above notaton, we can wrte the force balance equaton for lnk : F f f + 1 (151) We can, of course, express equaton (151) wth respect to any coordnate frame. Let us rewrte (151) n terms of coordnate frame {} : F F F f f + 1 f ( + 1R ) + 1 f + 1 f ( + 1R ) + 1 f + 1 (152) (153) (154) C. Moment balance equaton Now, let us wrte the moment balance equaton about the center of mass of lnk : N n n V 1 f V 2 f + 1 (155) Note that V 1 f and V 2 f gve the moments nduced by forces f + 1 and f, respectvely, about the + 1 center of mass of lnk {}. Let us now wrte expressons for V 1 and V 2 n lnk-specfc notaton: V 1 ( P ) C (156) V 2 ( P + 1 P ). (157) C Substtutng (156) and (157) nto (155) and expressng wth respect to frame {} : N n n ( P ) f ( P + 1 P ) f + 1 C C (158) Rearrangng terms and keepng equaton (152) n mnd, N N N N n n + 1 PC ( f f + 1) P + 1 f + 1 n n + 1 PC F P + 1 f + 1 n n + 1 PC F P [ ( + 1R) f + 1 n n + 1 PC F P + 1 ( + 1R ) + 1 [ f + 1 (159) (160) (161) (162) D. Inward teraton of lnk forces and moments Thus the force and moment balance equatons at lnk are gven by, F f ( + 1R ) + 1 f + 1 and (163) N n n + 1 PC F. (164) P + 1 ( + 1R ) + 1 [ f

15 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators We can rewrte equatons (163) and (164) as teratons that propagate f and n from the end-effector to the lnk frame { 1} : f n ( + 1R ) + 1 f F N ( + 1R ) + 1 n + 1 PC F P + 1 ( + 1R ) [ f + 1 (165) (166) Note that equatons (165) and (166) allow us to recursvely compute the forces and moments that each lnk exerts on ts neghborng lnks by nwardly propagatng from coordnate frame { N} to frame { 1}. The last remanng queston s, once equatons (165) and (166) are computed, what should be the torques/forces for the actuators to acheve the desred jont moton? ll components of the force and moment vectors f and n are ressted by the structure of the mechansm tself, except for the torque/force about/along the jont axs. Therefore the requred torque for a revolute jont s gven by, τ n Ẑ (167) whle the requred force for a prsmatc jont s gven by, τ f Ẑ. (168) 7. Complete formulaton of the teratve Newton-Euler dynamcs Ths secton summarzes the complete formulaton of the teratve Newton-Euler dynamcs model. It conssts of (1) the outward propagaton of angular veloctes and lnear and angular acceleratons, (2) the outward propagaton of net moments and forces actng on the lnks, and (3) the nward propagaton of forces and torques between lnks. Collectvely, equatons (170) through (182) mplctly defne the relatonshp we were lookng for at the begnnng of ths dscusson namely, τ h( ΘΘ,, Θ ). Outward teraton 1. Revolute jonts: (169) ω + 1 R + 1 ω + θ + 1 Ẑ ω + 1 R ω R ( ω) θ + 1 Ẑ θ + 1 Ẑ R ω P + 1 ω ω P [ + ( + 1) + 2. Prsmatc jonts: ω + 1 R ω ω + 1 R ω oth jont types: + 1 C R ω P + 1 ω ω P [ + ( + 1) R ( ω) d + 1 Ẑ ḋ + 1 Ẑ ω PC ω [ ω + 1 PC + 1 (170) (171) (172) (173) (174) (175) (176)

16 EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators F + 1 m + 1 v C + 1 C 1 C N I + 1 ω + 1 ω I + 1 ω + 1 (177) (178). Inward teraton 1. oth jont types: f n ( + 1R ) + 1 f F N ( + 1R ) + 1 n + 1 PC F P + 1 ( + 1R ) [ f + 1 (179) (180) 2. Revolute jonts: τ n Ẑ (181) 3. Prsmatc jonts: τ f Ẑ. (182) C. Intalzaton of propagatons In order to compute equatons (170) through (175), we need to know 0ω, ω 0 and 0 that s, the angular velocty, and lnear and angular acceleraton of the base coordnate frame { 0}. For a fxed-base manpulator, 0ω T, (183) 0ω T, and, (184) G, (185) 0 where G denotes the gravty vector. Note that (185) s equvalent to sayng that the base of the robot s acceleratng upward wth acceleraton g, and therefore easly ncorporates the effects of gravty loadng on the lnks wthout any adonal effort. N + 1 N + 1 In order to compute equatons (179) and (180), we need to know fn + 1 and nn + 1 that s, the forces and moments from the envronment actng on the end-effector of the manpulator. When the manpulator end-effector s not n contact wth any object or obstacle, these are smply gven by, N + 1 fn T, and, (186) N + 1 nn T. (187)

The classical spin-rotation coupling

The classical spin-rotation coupling LOUAI H. ELZEIN 2018 All Rghts Reserved The classcal spn-rotaton couplng Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 louaelzen@gmal.com Abstract Ths paper s prepared to show that a rgd

More information

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is. Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these

More information

Spin-rotation coupling of the angularly accelerated rigid body

Spin-rotation coupling of the angularly accelerated rigid body Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s

More information

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut MEV442 Introducton to Robotcs Module 2 Dr. Santhakumar Mohan Assstant Professor Mechancal Engneerng Natonal Insttute of Technology Calcut Jacobans: Veloctes and statc forces Introducton Notaton for tme-varyng

More information

Lecture 23: Newton-Euler Formulation. Vaibhav Srivastava

Lecture 23: Newton-Euler Formulation. Vaibhav Srivastava Lecture 23: Newton-Euler Formulaton Based on Chapter 7, Spong, Hutchnson, and Vdyasagar Vabhav Srvastava Department of Electrcal & Computer Engneerng Mchgan State Unversty Aprl 10, 2017 ECE 818: Robotcs

More information

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 4. Moton Knematcs 4.2 Angular Velocty Knematcs Summary From the last lecture we concluded that: If the jonts

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

Conservation of Angular Momentum = "Spin"

Conservation of Angular Momentum = Spin Page 1 of 6 Conservaton of Angular Momentum = "Spn" We can assgn a drecton to the angular velocty: drecton of = drecton of axs + rght hand rule (wth rght hand, curl fngers n drecton of rotaton, thumb ponts

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Physics 181. Particle Systems

Physics 181. Particle Systems Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

8.6 The Complex Number System

8.6 The Complex Number System 8.6 The Complex Number System Earler n the chapter, we mentoned that we cannot have a negatve under a square root, snce the square of any postve or negatve number s always postve. In ths secton we want

More information

So far: simple (planar) geometries

So far: simple (planar) geometries Physcs 06 ecture 5 Torque and Angular Momentum as Vectors SJ 7thEd.: Chap. to 3 Rotatonal quanttes as vectors Cross product Torque epressed as a vector Angular momentum defned Angular momentum as a vector

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

Chapter 11 Angular Momentum

Chapter 11 Angular Momentum Chapter 11 Angular Momentum Analyss Model: Nonsolated System (Angular Momentum) Angular Momentum of a Rotatng Rgd Object Analyss Model: Isolated System (Angular Momentum) Angular Momentum of a Partcle

More information

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force. Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act

More information

Physics 111: Mechanics Lecture 11

Physics 111: Mechanics Lecture 11 Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton

More information

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding. Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays

More information

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15 NGN 40 ynamcs and Vbratons Homework # 7 ue: Frday, Aprl 15 1. Consder a concal hostng drum used n the mnng ndustry to host a mass up/down. A cable of dameter d has the mass connected at one end and s wound/unwound

More information

Lecture 20: Noether s Theorem

Lecture 20: Noether s Theorem Lecture 20: Noether s Theorem In our revew of Newtonan Mechancs, we were remnded that some quanttes (energy, lnear momentum, and angular momentum) are conserved That s, they are constant f no external

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the

More information

Classical Mechanics ( Particles and Biparticles )

Classical Mechanics ( Particles and Biparticles ) Classcal Mechancs ( Partcles and Bpartcles ) Alejandro A. Torassa Creatve Commons Attrbuton 3.0 Lcense (0) Buenos Ares, Argentna atorassa@gmal.com Abstract Ths paper consders the exstence of bpartcles

More information

Lie Group Formulation of Articulated Rigid Body Dynamics

Lie Group Formulation of Articulated Rigid Body Dynamics Le Group Formulaton of Artculated Rgd Body Dynamcs Junggon Km 11/9/2012, Ver 1.01 Abstract It has been usual n most old-style text books for dynamcs to treat the formulas descrbng lnearor translatonal

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force. The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,

More information

Introduction To Robotics (Kinematics, Dynamics, and Design)

Introduction To Robotics (Kinematics, Dynamics, and Design) ntroducton To obotcs Kneatcs, Dynacs, and Desgn SESSON # 6: l Meghdar, Professor School of Mechancal Engneerng Sharf Unersty of Technology Tehran, N 365-9567 Hoepage: http://eghdar.sharf.edu So far we

More information

10. Canonical Transformations Michael Fowler

10. Canonical Transformations Michael Fowler 10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst

More information

Translational Equations of Motion for A Body Translational equations of motion (centroidal) for a body are m r = f.

Translational Equations of Motion for A Body Translational equations of motion (centroidal) for a body are m r = f. Lesson 12: Equatons o Moton Newton s Laws Frst Law: A artcle remans at rest or contnues to move n a straght lne wth constant seed there s no orce actng on t Second Law: The acceleraton o a artcle s roortonal

More information

Physics 207: Lecture 27. Announcements

Physics 207: Lecture 27. Announcements Physcs 07: ecture 7 Announcements ake-up labs are ths week Fnal hwk assgned ths week, fnal quz next week Revew sesson on Thursday ay 9, :30 4:00pm, Here Today s Agenda Statcs recap Beam & Strngs» What

More information

σ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review

σ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review Chapter Four States of Stress Part Three When makng your choce n lfe, do not neglect to lve. Samuel Johnson Revew When we use matrx notaton to show the stresses on an element The rows represent the axs

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

Classical Mechanics Virtual Work & d Alembert s Principle

Classical Mechanics Virtual Work & d Alembert s Principle Classcal Mechancs Vrtual Work & d Alembert s Prncple Dpan Kumar Ghosh UM-DAE Centre for Excellence n Basc Scences Kalna, Mumba 400098 August 15, 2016 1 Constrants Moton of a system of partcles s often

More information

Rigid body simulation

Rigid body simulation Rgd bod smulaton Rgd bod smulaton Once we consder an object wth spacal etent, partcle sstem smulaton s no longer suffcent Problems Problems Unconstraned sstem rotatonal moton torques and angular momentum

More information

Rotational Dynamics. Physics 1425 Lecture 19. Michael Fowler, UVa

Rotational Dynamics. Physics 1425 Lecture 19. Michael Fowler, UVa Rotatonal Dynamcs Physcs 1425 Lecture 19 Mchael Fowler, UVa Rotatonal Dynamcs Newton s Frst Law: a rotatng body wll contnue to rotate at constant angular velocty as long as there s no torque actng on t.

More information

Modeling of Dynamic Systems

Modeling of Dynamic Systems Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how

More information

Chapter 11: Angular Momentum

Chapter 11: Angular Momentum Chapter 11: ngular Momentum Statc Equlbrum In Chap. 4 we studed the equlbrum of pontobjects (mass m) wth the applcaton of Newton s aws F 0 F x y, 0 Therefore, no lnear (translatonal) acceleraton, a0 For

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 151 Lecture 3 Lagrange s Equatons (Goldsten Chapter 1) Hamlton s Prncple (Chapter 2) What We Dd Last Tme! Dscussed mult-partcle systems! Internal and external forces! Laws of acton and

More information

coordinates. Then, the position vectors are described by

coordinates. Then, the position vectors are described by Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Three views of mechanics

Three views of mechanics Three vews of mechancs John Hubbard, n L. Gross s course February 1, 211 1 Introducton A mechancal system s manfold wth a Remannan metrc K : T M R called knetc energy and a functon V : M R called potental

More information

ENGI9496 Lecture Notes Multiport Models in Mechanics

ENGI9496 Lecture Notes Multiport Models in Mechanics ENGI9496 Moellng an Smulaton of Dynamc Systems Mechancs an Mechansms ENGI9496 Lecture Notes Multport Moels n Mechancs (New text Secton 4..3; Secton 9.1 generalzes to 3D moton) Defntons Generalze coornates

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

T f. Geometry. R f. R i. Homogeneous transformation. y x. P f. f 000. Homogeneous transformation matrix. R (A): Orientation P : Position

T f. Geometry. R f. R i. Homogeneous transformation. y x. P f. f 000. Homogeneous transformation matrix. R (A): Orientation P : Position Homogeneous transformaton Geometr T f R f R T f Homogeneous transformaton matr Unverst of Genova T f Phlppe Martnet = R f 000 P f 1 R (A): Orentaton P : Poston 123 Modelng and Control of Manpulator robots

More information

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics) CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O

More information

Linear Momentum. Center of Mass.

Linear Momentum. Center of Mass. Lecture 6 Chapter 9 Physcs I 03.3.04 Lnear omentum. Center of ass. Course webste: http://faculty.uml.edu/ndry_danylov/teachng/physcsi Lecture Capture: http://echo360.uml.edu/danylov03/physcssprng.html

More information

Angular Momentum and Fixed Axis Rotation. 8.01t Nov 10, 2004

Angular Momentum and Fixed Axis Rotation. 8.01t Nov 10, 2004 Angular Momentum and Fxed Axs Rotaton 8.01t Nov 10, 2004 Dynamcs: Translatonal and Rotatonal Moton Translatonal Dynamcs Total Force Torque Angular Momentum about Dynamcs of Rotaton F ext Momentum of a

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

MAE140 - Linear Circuits - Winter 16 Midterm, February 5

MAE140 - Linear Circuits - Winter 16 Midterm, February 5 Instructons ME140 - Lnear Crcuts - Wnter 16 Mdterm, February 5 () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator

More information

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

PHYS 1443 Section 003 Lecture #17

PHYS 1443 Section 003 Lecture #17 PHYS 144 Secton 00 ecture #17 Wednesda, Oct. 9, 00 1. Rollng oton of a Rgd od. Torque. oment of Inerta 4. Rotatonal Knetc Energ 5. Torque and Vector Products Remember the nd term eam (ch 6 11), onda, Nov.!

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics N40: ynamcs and Vbratons Homewor 7: Rgd Body Knematcs School of ngneerng Brown Unversty 1. In the fgure below, bar AB rotates counterclocwse at 4 rad/s. What are the angular veloctes of bars BC and C?.

More information

Part C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis

Part C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis Part C Dynamcs and Statcs of Rgd Body Chapter 5 Rotaton of a Rgd Body About a Fxed Axs 5.. Rotatonal Varables 5.. Rotaton wth Constant Angular Acceleraton 5.3. Knetc Energy of Rotaton, Rotatonal Inerta

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation 1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve

More information

Study Guide For Exam Two

Study Guide For Exam Two Study Gude For Exam Two Physcs 2210 Albretsen Updated: 08/02/2018 All Other Prevous Study Gudes Modules 01-06 Module 07 Work Work done by a constant force F over a dstance s : Work done by varyng force

More information

On the Hessian of Shape Matching Energy

On the Hessian of Shape Matching Energy On the Hessan of Shape Matchng Energy Yun Fe 1 Introducton In ths techncal report we derve the analytc form of the Hessan matrx for shape matchng energy. Shape matchng (Fg. 1) s a useful technque for meshless

More information

A NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT

A NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON

More information

CHAPTER 10 ROTATIONAL MOTION

CHAPTER 10 ROTATIONAL MOTION CHAPTER 0 ROTATONAL MOTON 0. ANGULAR VELOCTY Consder argd body rotates about a fxed axs through pont O n x-y plane as shown. Any partcle at pont P n ths rgd body rotates n a crcle of radus r about O. The

More information

11. Dynamics in Rotating Frames of Reference

11. Dynamics in Rotating Frames of Reference Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons

More information

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George

More information

Dynamics of a Spatial Multibody System using Equimomental System of Point-masses

Dynamics of a Spatial Multibody System using Equimomental System of Point-masses Dynamcs of a Spatal Multbody System usng Equmomental System of Pont-masses Vnay Gupta Dept. of Mechancal Engg. Indan Insttute of echnology Delh New Delh, Inda vnayguptaec@gmal.com Hmanshu Chaudhary Dept.

More information

Spring 2002 Lecture #13

Spring 2002 Lecture #13 44-50 Sprng 00 ecture # Dr. Jaehoon Yu. Rotatonal Energy. Computaton of oments of nerta. Parallel-as Theorem 4. Torque & Angular Acceleraton 5. Work, Power, & Energy of Rotatonal otons Remember the md-term

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Lesson 5: Kinematics and Dynamics of Particles

Lesson 5: Kinematics and Dynamics of Particles Lesson 5: Knematcs and Dynamcs of Partcles hs set of notes descrbes the basc methodology for formulatng the knematc and knetc equatons for multbody dynamcs. In order to concentrate on the methodology and

More information

LAB # 4 - Torque. d (1)

LAB # 4 - Torque. d (1) LAB # 4 - Torque. Introducton Through the use of Newton's three laws of moton, t s possble (n prncple, f not n fact) to predct the moton of any set of partcles. That s, n order to descrbe the moton of

More information

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have

More information

Physics 106a, Caltech 11 October, Lecture 4: Constraints, Virtual Work, etc. Constraints

Physics 106a, Caltech 11 October, Lecture 4: Constraints, Virtual Work, etc. Constraints Physcs 106a, Caltech 11 October, 2018 Lecture 4: Constrants, Vrtual Work, etc. Many, f not all, dynamcal problems we want to solve are constraned: not all of the possble 3 coordnates for M partcles (or

More information

PHYS 705: Classical Mechanics. Canonical Transformation II

PHYS 705: Classical Mechanics. Canonical Transformation II 1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m

More information

From Biot-Savart Law to Divergence of B (1)

From Biot-Savart Law to Divergence of B (1) From Bot-Savart Law to Dvergence of B (1) Let s prove that Bot-Savart gves us B (r ) = 0 for an arbtrary current densty. Frst take the dvergence of both sdes of Bot-Savart. The dervatve s wth respect to

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Implicit Integration Henyey Method

Implicit Integration Henyey Method Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure

More information

AP Physics 1 & 2 Summer Assignment

AP Physics 1 & 2 Summer Assignment AP Physcs 1 & 2 Summer Assgnment AP Physcs 1 requres an exceptonal profcency n algebra, trgonometry, and geometry. It was desgned by a select group of college professors and hgh school scence teachers

More information

NEWTON S LAWS. These laws only apply when viewed from an inertial coordinate system (unaccelerated system).

NEWTON S LAWS. These laws only apply when viewed from an inertial coordinate system (unaccelerated system). EWTO S LAWS Consder two partcles. 1 1. If 1 0 then 0 wth p 1 m1v. 1 1 2. 1.. 3. 11 These laws only apply when vewed from an nertal coordnate system (unaccelerated system). consder a collecton of partcles

More information

How Differential Equations Arise. Newton s Second Law of Motion

How Differential Equations Arise. Newton s Second Law of Motion page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons

More information

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso Supplement: Proofs and Techncal Detals for The Soluton Path of the Generalzed Lasso Ryan J. Tbshran Jonathan Taylor In ths document we gve supplementary detals to the paper The Soluton Path of the Generalzed

More information

Poisson brackets and canonical transformations

Poisson brackets and canonical transformations rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order

More information

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Professor Terje Haukaas University of British Columbia, Vancouver  The Q4 Element Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to

More information

Grover s Algorithm + Quantum Zeno Effect + Vaidman

Grover s Algorithm + Quantum Zeno Effect + Vaidman Grover s Algorthm + Quantum Zeno Effect + Vadman CS 294-2 Bomb 10/12/04 Fall 2004 Lecture 11 Grover s algorthm Recall that Grover s algorthm for searchng over a space of sze wors as follows: consder the

More information

Technical Report TR05

Technical Report TR05 Techncal Report TR05 An Introducton to the Floatng Frame of Reference Formulaton for Small Deformaton n Flexble Multbody Dynamcs Antono Recuero and Dan Negrut May 11, 2016 Abstract Ths techncal report

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

More information

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be

More information