Lie Group Formulation of Articulated Rigid Body Dynamics
|
|
- Henry Chambers
- 5 years ago
- Views:
Transcription
1 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 and angularor rotatonal moton of a rgd body separately. For example, the famous Newton s 2nd law, f = ma, for the translatonal moton of a rgd body has ts partner, so-called the Euler s equaton whch descrbes the rotatonal moton of the body. In ths artcle a smple and elegant equatons of moton of a rgd body s presented usng generalzed notatons of knematc and dynamc quanttes such as velocty, force and momentum Secton 1. Then effcent recursve algorthms for nverse and forward dynamcs of artculated rgd bodes are derved usng the notatons. A hybrd dynamcs algorthm, or a generalzed verson of the nverse and forward dynamcs algorthms, s also presented Secton 2. Fnally, analytc dervatves of the dynamcs algorthms are derved Secton Dynamcs of a Rgd Body Ths secton descrbes the equatons of moton of a sngle rgd body n a geometrc manner. 1.1 Rgd Body Moton Fgure 1: Coordnates frames for a rgd body To descrbe the moton of a rgd body, we need to represent both the poston and orentaton of the body. Let B} be a coordnate frame attached to the rgd body and A} be an arbtrary coordnate frame, and all coordnate frames wll be rght-handed Cartesan from now on. We can defne a 3 3 matrx R = [x ab, y ab, z ab ] 1 1 GEAR Geometrc Engne for Artculated Rgd-body smulaton s a C++ mplementaton of the algorthms presented n ths artcle. 1
2 where x ab, y ab, z ab R 3 are the coordnates of the coordnate axes of B} wth respect to A}. A matrx of ths form s called a rotaton matrx as t can be used to descrbe the orentatonor rotaton of a rgd body, relatve to a reference frame. Snce the columns of a rotaton matrx are mutually orthonormal and the coordnate frame s rght-handed, the rotaton matrx has two propertes, RR T = R T R = I, det R = 1, 2 and t s denoted by SO3 2. Let p R 3 be the poston vector of the orgn of B} from the orgn of A}, and R SO3 be the rotaton matrx of B} relatve to A}. The confguraton space of the rgd body moton can be represented wth the par R, p, whch s denoted as SE3. A 4 4 matrx, [ ] R p T = s called the homogeneous representaton of T = R, p SE3, and ts nverse can be obtaned wth [ ] T R = T R T p From now on, a smple declaraton, T SE3 : A} B}, wll be used to notfy that T SE3 represents the orentaton and poston of a coordnate frame B} wth respect to another coordnate frame A}. The Le algebra of SE3, denoted as se3, s dentfed as a 6-dmensonal vector space w, v R 6 where w so3, the Le algebra of SO3. ξ = w, v se3 can also be represented as a 4 4 matrx, [ ] [w] v ξ = ] where [w] = [ 0 w3 w 2 w 3 0 w 1 w 2 w 1 0 R 3 3 s a skew-symmetrc matrx. The adjont acton of T SE3 on ξ se3, Ad : SE3 se3 se3, s defned as Ad T ξ = T ξ T. 6 From 3, 4 and 5, Ad T can be regarded as a lnear transformaton, Ad T : se3 se3, whch s defned by a 6 6 matrx [ ] R 0 Ad T = 7 [p]r R where T = R, p SE3. The coadjont acton of T on ξ dse3 whch s the dual of ξ, Ad T : dse3 dse3, s defned by a 6 6 matrx 1.2 Generalzed Velocty and Force Ad T = Ad T T. 8 Let T t = Rt, pt SE3 be a moton trajectory of a coordnate frame attached to a rgd body wth respect to an nertal frame. The generalzed velocty of the rgd body s defned as [ ] V = T T [w] v = where [w] = R T Ṙ and v = R T ṗ. The physcal meanng of w R 3 s the rotatonalor angular velocty of the coordnate frame attached to the body relatve to the nertal frame, but expressed n 2 The notaton SO abbrevates specal orthogonal and specal refers to the fact that det R = +1 rather than ±1. See [2] for more detals. 2
3 the body coordnate frame. Smlarly, v R 3 represents the velocty of the orgn of the coordnate frame relatve to the nertal frame, and stll expressed n the body frame. The generalzed velocty s an element of se3, and can be smply regarded as a 6-dmensonal vector,.e., w V =. 10 v As the generalzed velocty s an nstance of se3, t follows the adjont transformaton rule defned n 7. Let A}, B} be two dfferent coordnate frames attached to the same rgd body, and T a, T b SE3 represent the orentaton and poston of the two frames wth respect to an nertal frame. Then, from 6 and 9, the generalzed veloctes of A} and B} have the followng relaton: V b = Ad Tba V a 11 where T ba SE3 : B} A}. Wth a coordnate frame attached to a rgd body, the generalzed force actng on the body can be defned as m F = 12 f where m R 3 and f R 3 represent a moment and force actng on the body respectvely, vewed n the body frame. The generalzed force s known as the member of dse3, and has the followng transformaton rule, F b = Ad T ab F a 13 where F a and F b denote a generalzed force vewed from dfferent body frames A} and B}, and T ab SE3 : A} B}. 1.3 Generalzed Inerta and Momentum The knetc energy of a rgd body s gven by the followng volume ntegral 1 e = 2 v 2 dm 14 vol whch means the sum of the knetc energes of all the mass partcles consttutng the body. By ntroducng a coordnate frame attached to the body, 14 can be restructured as the followng smple quadratc form, e = 1 2 V T IV 15 where V R 6 s the generalzed velocty of the body and I R 6 6, whch s known as generalzed nerta, represents the mass and mass dstrbuton wth respect to the body frame. To obtan an explct form of the generalzed nerta of a rgd body, let r R 3 be the poston of a body pont relatve to the body frame and R, p SE3 represents orentaton and poston of the body frame wth respect to an nertal frame respectvely. Usng v 2 = ṗ+ṙr 2, R T Ṙ = [w] and R T ṗ = v, 14 can be rewrtten as e = 1 ṗ 2 + 2ṗ T Ṙr + 2 Ṙr 2 dm 16 vol = 1 } ṗ T ṗ dm 2ṗ T R [r]dm w + w T [r] T [r]dm w 17 2 vol vol vol } mv T v 2v T w + w T = 1 2 [r]dm vol [r] T [r]dm vol = 1 2 V T IV 19 w 18 3
4 where V = w, v s the generalzed velocty of the body and the generalzed nerta, I, has the followng explct matrx form: [ vol I = [r]t [r]dm vol [r]dm ]. 20 vol [r]t dm m1 The generalzed nerta s symmetrc postve defnte, and ts upper dagonal term, I = vol [r]t [r]dm, s the defnton of the well-known 3 3 nerta matrx of the rgd body wth respect to the body frame. If the orgn of the body frame s located on the center of mass, then the generalzed nerta becomes a block dagonal matrx because [r]dm = 0. In addton, f the orentaton of the body vol frame also concdes wth the prncple axes of the body, then the generalzed nerta becomes a dagonal matrx. Let A} and B} be coordnate frames attached to a rgd body, I a and I b be the generalzed nertas of the body correspondng to the two frames. Usng 11, 15, and the fact that the knetc energy of the body should reman under change of coordnate frame, the followng transformaton rule between the generalzed nertas can be obtaned: I b = Ad T ab I a Ad Tab 21 where T ab SE3 : A} B}. If the mass, m, and the nerta matrx n a center of mass frame 3, I c R 3 3, are gven, one can get the generalzed nerta n an arbtrary body frame from 21, rather than usng 20, as [ ] RI I = c R T +m[p] T [p] m[p] m[p] T 22 m1 where R, p SE3 represents the orentaton and poston of the center of mass frame wth respect to the body frame. The generalzed momentum of a rgd body s defned as L = IV 23 where I and V are the generalzed nerta and velocty of the body expressed n a coordnate frame attached to the body. L a and L b be the generalzed momentum of a rgd body expressed n dfferent body frames, A} and B}, respectvely. Usng 11 and 21 one can derve the followng transformaton rule for generalzed momentums: L b = Ad T ab L a 24 whch s same to that of generalzed forces, and ndeed, the generalzed momentum s also known as dse Tme Dervatves of se3 and dse3 Recall that the tme dervatve of a 3-dmensonal vector x = 3 =1 x ê, expressed n a movng coordnate frame ê }, can be obtaned as d dt x = 3 =1 } d d dt x ê + x = ẋ + w x 25 dtê where ẋ = 3 d =1 dt x ê = dx 1 dt, dx2 dt, dx3 dt s the component-wse tme dervatve of x and w s the angular velocty of the movng frame. 3 A coordnate frame whose orgn s located on the center of mass of the body. 4
5 The tme dervatves of se3 and dse3 whch are expressed n a movng frame have more generalzed form. A smple approach to obtan the dervatves s to transform se3or dse3 to a statonary coordnate frame, dfferentate t there, and transform the dervatve back to the orgnal movng coordnate frame wth the correct assumpton that the dervatve of se3or dse3 can be transformed wth the rule of se3or dse3. 4 Note that dfferentatng a vector n a statonary coordnate frame wth respect to tme s just gettng the component-wse dervatve of t. Lemma 1. Let X se3 be expressed n a movng frame attached to a rgd body. Then the tme dervatve of X can be obtaned by d dt X = Ẋ + ad V X 26 where Ẋ s the component-wse tme dervatve of X and ad V : se3 se3 s a lnear transformaton defned as [ ] [w] 0 ad V = 27 [v] [w] where V = w, v se3 s the generalzed velocty of the body. Proof. Let T = R, p SE3 denote the orentaton and poston of the movng frame wth respect to an nertal frame whch s statonary n the space. By transformng X to the nertal frame, dfferentatng t there, and then transformng the result back to the orgnal body frame, one can get d dt X = Ad d T dt Ad T X = Ẋ + Ad d T dt Ad T X where d dt d represents the component-wse dfferentaton and Ẋ = dtx. Usng 7 and V = w, v = as follows: R T Ṙ, R T ṗ, one can show Ad T d dt Ad T = ad V d [ ] [ ] R Ad T dt Ad T 0 Ṙ 0 T = R T [p] R T [ṗ] R + [p] [ ] Ṙ Ṙ [w] 0 = [v] [w] Lemma 2. Let Y dse3 be expressed n a movng frame attached to a rgd body. Then the tme dervatve of Y can be obtaned by d dt Y = Ẏ ad V Y 28 where Ẏ s the component-wse tme dervatve of Y and ad V : dse3 dse3 s a lnear transformaton defned as [ ] T ad V = ad T [w] 0 V = 29 [v] [w] where V = w, v se3 s the generalzed velocty of the body. Proof. Smlarly to the proof of Lemma 1, the dervatve of Y dse3 can be obtaned by d dt Y = d Ad T dt Ad T Y = Ẏ + d Ad T dt Ad T Y. 4 In [1] the tme dervatve of a spatal velocty whch s smlar to the generalzed velocty n ths artcle was obtaned wth ths approach. 5
6 Usng 8 and V = w, v = R T Ṙ, R T ṗ, one can show Ad T d dt Ad T = ad V as follows: Ad d [ ] [ ] R T dt Ad T = T R T [p] Ṙ [ṗ] R + [p] Ṙ 0 R T 0 Ṙ [ ] [w] [v] = = ad T 0 [w] V = ad V 1.5 Geometrc Dynamcs of a Rgd Body Equatons of moton of a rgd body can be wrtten as F = d dt L 30 where F represents the net sum of the generalzed forces actng on the rgd body and L s the generalzed momentum of the body. Usng 23 and 28, the equatons of moton of the rgd body can be wrtten as F = I V ad V IV 31 where V s the component-wse tme dervatve of the generalzed velocty of the body. Note that L = İV + I V = I V because the components of I doesn t vary. The dynamcs equaton of a rgd body s coordnate nvarant,.e., the structure of the equaton stll remans under change of coordnate frame. Let A} and B} be coordnate frames attached to a body, V a and V b be the generalzed veloctes, F a and F b be the generalzed forces, and I a and I b be the generalzed nertas correspondng to A} and B} respectvely. Usng the transformaton rules, 11, 13, and 21, one can easly transform the dynamcs equatons wth respect to A}, F a = I a V a ad V a I a V a, to the equatons n B}, F b = I b V b ad V b I b V b, and ths shows the coordnate nvarance of 31 under change of coordnate frame. 2 Dynamcs of Open Chan Systems Let q R n denote the set of coordnates of all jonts n a system, and for open chan systems, n s equal to the degree-of-freedom of the system. The dynamcs equatons of the system can be wrtten as Mq q + bq, q = τ 32 where Mq R n n s a symmetrc mass matrx of the system, bq, q R n represents Corols, centrfugal, and gravty terms, and τ R n denotes torqueor force vector correspondng to the system coordnates q. Calculatng the jont torqueor force τ wth gven q, q, q s called nverse dynamcs. It s typcally used to obtan the requred jont torques whch make the system move along a gven jont trajectory. On the other hand, calculatng q wth gven q, q, τ s called forward dynamcs, and t s frequently used to smulate the system evoluton n tme by ntegratng q to get q, q at the next tme step. In general, the command nput on a jont can be ether torque or acceleraton durng the smulaton and the command type does not have to be same for all jonts. The equatons of moton can be rewrtten as M qa q p + b = τa τ p where the subscrpt a s for the acceleraton-prescrbed jonts and the subscrpt p s for the jonts wth gven, or known, torques. We can compute τ a, q p wth known q a, τ p from the equatons and 33 6
7 Inverse dynamcs q, q, q τ Forward dynamcs q, q, τ q Hybrd dynamcs q a, q a, q a, q p, q p, τ p τ a, q p Table 1: Input and output of dynamcs we call ths hybrd dynamcs 5 One possble soluton for hybrd dynamcs s to rearrange 33 and solve t wth a drect matrx nverson. For example, the acceleratons [ of the unprescrbed jonts can be obtaned by q p = Mpp Maa M τ p b p M pa q a where M = ap M pa M pp ], b = ba bp, and q = qa q p. The method, however, s not effcent for a complex system because t requres buldng the mass matrx and nvertng the submatrx correspondng to the unprescrbed jonts, whch leads to an On 2 + On 3 p algorthm where n and n p denote the number of all coordnates and the number of unprescrbed coordnates respectvely. In Table 1 the nput and output of nverse, forward, and hybrd dynamcs are summarzed. Note that hybrd dynamcs s a generalzaton of tradtonal forward and nverse dynamcs,.e., they can be regarded as the extreme cases of hybrd dynamcs when all of the jonts have gven, or known, torques and when all of the jonts are acceleraton-prescrbed respectvely. 2.1 Recursve Inverse Dynamcs A recursve Newton-Euler nverse dynamcs algorthm usng the geometrc notatons shown n Secton 1 was presented n [3]. Here the algorthm s slghtly modfed to support mult-degree-offreedom jonts n the formulaton. Let 0} be an nertal frame whch s statonary n the space, I} be a coordnate frame attached to -th rgd body of the open chan system, and λ} be a coordnate frame attached to the parent body of the -th rgd body. Also, let T SE3 : 0} }, T λ SE3 : 0} λ}, and T SE3 : λ} }. From T = T λ T and 6, the generalzed velocty of the -th body can be rewrtten as V = T T 34 = T T λ T λ T + T λ T 35 = Ad T V λ + S q 36 where S q = T T se3 represents the relatve velocty of the -th body wth respect to ts parent. S = S q se3 n s called the Jacoban of the jont connectng the -th body and ts parent and q R n represents the coordnate vector of the jont. As shown n 31, component-wse tme dervatves of the generalzed veloctes of all bodes n the system are needed to buld the dynamcs equatons for each body. Recallng that A = A AA for an arbtrary matrx A, and ad ξ1 ξ 2 = ξ 1 ξ 2 ξ 2 ξ 1 for arbtrary ξ 1, ξ 2 se3, one can derve the followng formula for V, the component-wse tme dervatve of V : V = d T dt V λt + Ṡ q + S q 37 = T T T V λt + T V λ T + T V λ T + Ṡ q + S q 38 = Ad T V λ + ad AdT V λ S q + Ṡ q + S q 39 = Ad T 5 We follow [1] for the termnology. V λ + ad V S q + Ṡ q + S q 40 7
8 36 and 40 are well suted to calculate the generalzed velocty and ts component-wse tme dervatve of each body n a open chan system from the ground to the end of the system recursvely, as the velocty and acceleraton of the ground are known n most cases. 6 Note that Ṡ 0 as the jont Jacoban S s a functon of q n general. Let F dse3 be the generalzed force transmtted to -th body from ts parent through the connectng jont, and F ext Both F and F ext be wrtten as dse3 be a generalzed force actng on the -th body from envronment. are expressed n }. From 31, the equatons of moton for the -th body can F + F ext Ad F T k = I V ad V I V 41,k k µ where the left hand sde of the equatons represents the net force actng on the body, µ s the set of chld bodes of the -th body, and Ad F T k s the generalzed force, transmtted from k-th,k chld body, expressed n }. It should be noted that the generalzed force, F, for each body can be calculated by 41 from the ends to the ground recursvely, as the end bodes have no chld. A recursve nverse dynamcs algorthm for open chan systems s shown n Table 2, and the followng s a lst of symbols for the geometrc nverse dynamcs: = ndex of the -th body. λ = ndex of the parent body of the -th body. µ = set of ndexes of the chld bodes of the -th body. q R n τ R n = coordnates of the -th jont whch connects the -th body wth ts parent body. = torqueor force exerted by the -th jont. T SE3 : λ} }, a functon of q. V se3 = the generalzed velocty of the -th body, vewed n the body frame }. V se3 = component-wse tme dervatve of V. S se3 [ n = Jacoban of T ] vewed n }. S = T T,, T T q 1 q n, where q k R denotes the k-th coordnate of the -th jont,.e., q = q 1,, qn. I = the generalzed nerta of the -th body, vewed n }. F dse3 = the generalzed force transmtted to the -th body from ts parent through the connectng jont, vewed n }. F ext }. dse3 = the generalzed force actng on the -th body from envronment, vewed n Table 2: Recursve Inverse Dynamcs 6 One can assume that V 0 = 0 and V 0 = 0, g where g R 3 denote the gravty vector, vewed n the nertal frame, wth approprate drecton and magntude. 8
9 whle forward recurson do T = functon of q V = Ad T V λ + S q V = Ad T V λ + ad V S q + Ṡ q + S q end whle whle backward recurson do F = I V ad V I V F ext τ = S T F end whle 2.2 Recursve Forward Dynamcs + k µ Ad F T k,k Featherstone [1] found that the dynamcs equatons of the -th body can be reformulated to have the followng form, F = Î V + ˆB, 42 where Î s called as the artculated body nerta of the body and ˆB s an assocated bas force. He also showed that the artculated body nerta and bas force correspondng to each body n open chan systems can be calculated recursvely, and by usng these new quanttes, forward dynamcs can be solved wth an On algorthm. A Le group formulaton of the artculated body nerta method was reported n [4]. Here, a more general form of the geometrc formulaton supportng mult-degree-of-freedom jont models s presented. Startng from 42, we ll show that the same form of dynamcs equatons stll holds for the parent of the -th bodyλ-th body, and Îλ and ˆB λ can be calculated from Î and ˆB, whch leads a backward recurson process for them. Let s assume that the equatons of moton for -th body can be wrtten as 42. By substtutng 40 for V n 42, one can get F = Î Ad T V λ + ad V S q + S q + Ṡ q + ˆB, 43 and from S T F = τ, the unknown q can be wrtten as q = S T ÎS τ S T Î Ad T V λ + ad V S q + Ṡ q S T ˆB }. 44 From 41 the dynamcs equatons for the λ-th body becomes F λ = I λ Vλ ad V λ I λ V λ Fλ ext + and by substtutng 43 for F k n 45, one can get k µλ Ad F T k, 45 λ,k F λ = I λ Vλ ad V λ I λ V λ Fλ ext + Ad Î T k Ad λ,k T V λ + ad Vk S k q k + S k q k + Ṡk q k + ˆB k }. λ,k k µλ 46 By substtutng 44 for the unknown q k n 46 and arrangng the equatons n terms of V λ, one can have the dynamcs equatons for the λ-th body wth the followng desred form: F λ = Îλ V λ + ˆB λ 47 9
10 where Î λ = I λ + k µλ k µλ } Ad Î T k ÎkS k S Tk ÎkS k S Tk Î k Ad λ,k T 48 λ,k ˆB λ = ad V λ I λ V λ Fλ ext + Ad ˆB T k + ad Îk Vk S k q k + Ṡk q k λ,k + ÎkS k S T k ÎkS k τ k S T k Îkad Vk S k q k + Ṡk q k S T k ˆB k }. 49 In summary, the On algorthm for forward dynamcs of open chan systems conssts of the followng three man recurson process: 1. Forward recurson: recursvely calculates V for each body wth Backward recurson: recursvely calculates Î and ˆB for each body wth 48 and Forward recurson: recursvely calculates q and V for each body wth 44 and 40. Table 3 shows the forward dynamcs algorthm for open chan systems wth a few addtonal ntermedate varables, such as η, Ψ, Π, and β, for the smplcty and effcency of the equatons. Table 3: Recursve Forward Dynamcs whle forward recurson do T = functon of q V = Ad T V λ + S q η = ad V S q + Ṡ q end whle whle backward recurson do Î = I + k µ Ad T,k ˆB = ad V I V F ext Π k Ad T,k + k µ Ad T,k Ψ = S TÎS Π = Î ÎS Ψ S TÎ β = ˆB + η Î + S Ψ τ S T Î η + ˆB } end whle whle forward recurson do } q = Ψ τ S TÎ Ad T V λ + η S T ˆB V = Ad T V λ + S q + η F = Î V + ˆB end whle 2.3 Recursve Hybrd Dynamcs A geometrc recursve algorthm for hybrd dynamcs of open chan systems was reported n [5], but the formulaton was based on 1-dof jonts. Here, more general form of the geometrc hybrd dynamcs s presented to support mult-degree-of-freedom jonts more convenently. One can derve the hybrd dynamcs for open chan systems wth a smlar way as n 2.2 for forward dynamcs. Let s go back to 46. Unlke the case of forward dynamcs, as q k n 46 β k 10
11 correspondng to an actve jont s already known, one can arrange 46 n terms of Vλ by substtutng 44 nto 46 for q k of passve jonts only as follows: F λ = Îλ V λ + ˆB λ 50 where Î λ = I λ + + k µλ a} k µλ p} k µλ p} Ad Î T k Ad λ,k T λ,k } Ad Î T k ÎkS k S Tk ÎkS k S Tk Î k Ad λ,k T λ,k ˆB λ = ad V λ I λ V λ Fλ ext + Ad ˆB T k + ad } Îk Vk S k q k + Ṡk q k + S k q k λ,k k µλ a} + Ad ˆB T k + ad Îk Vk S k q k + Ṡk q k λ,k + ÎkS k S T k ÎkS k τ k S T k Îkad Vk S k q k + Ṡk q k S T k ˆB k } where a and p denote the sets of the acceleraton-prescrbed and unprescrbed jonts n the system respectvely. Table 4 shows the hybrd dynamcs algorthm for open chan systems. Table 4: Recursve Hybrd Dynamcs whle forward recurson do T = functon of q V = Ad T V λ + S q η = ad V S q + Ṡ q end whle whle backward recurson do Î = I + k µ Ad T,k ˆB = ad V I V F ext f a then Π = Î β = ˆB + Î η + S q else Ψ = S TÎS Π = Î ÎS Ψ S TÎ Π k Ad T,k + k µ Ad T,k β = ˆB + η Î + S Ψ τ S T Î η + ˆB } end f end whle whle forward recurson do f a then V = Ad T V λ + S q + η F = Î V + ˆB τ = S TF else q = Ψ τ S TÎ Ad T V λ + η S T β k ˆB } 11
12 V = Ad T V λ + S q + η F = Î V + ˆB end f end whle 3 Dfferentaton of the Geometrc Dynamcs 3.1 Basc Dervatves It s useful to have the followng dervatves for dfferentatng the recursve dynamcs algorthms wth the chan rule. Lemma 3. Let p R be an arbtrary scalar varable and T SE3 be a functon of p. Then Ad T = ad T T Ad T. 53 Proof. Let ξ be an arbtrary se3. If Ad T s regarded as a 6 6 matrx, then Ad T ξ = Ad T ξ + Ad T ξ. 54 Ad T can also be regarded as a lnear mappng, Ad T : ξ T ξt, and n ths case, Ad T ξ = T ξt 55 = T ξt + T ξ T T + T ξ 56 = T T T ξt T T ξt T + T ξ T 57 = ad T T Ad T ξ + Ad T ξ. 58 As ξ se3 s arbtrary, one can see Ad T = ad T T Ad T from 54 and 58. Corollary 1. Let p R be an arbtrary scalar varable and T SE3 : λ} }. Then Ad T = ad h Ad T Ad = Ad ad T T h where h se3 s defned as h = T T. 61 Corollary 2. If p = q k where q k denotes the k-th coordnate of the -th jont, then h = Sk 62 where S k se3 denotes the k-th column of the -th jont Jacoban, S se3 n. If p / q = q 1,, qn }, then h =
13 3.2 Dervatves of the Dynamcs By applyng chan rule wth 59 and 60, the recursve algorthm for nverse, forward, and hybrd dynamcs can be dfferentated wth respect to an arbtrary scalar varable p R. Table 5 shows the dervatve of the recursve nverse dynamcs, and t can solve τ wth gven q, q, q. In Table 6, the dervatve of the recursve forward dynamcs, whch calculates q τ wth gven, s gven. The dervatve of the recursve hybrd dynamcs whch solves s presented n Table 7. τ a, qp Table 5: Dervatve of the Recursve Inverse Dynamcs whle forward recurson do T T h V V = Ad T = Ad T + ad V V λ V λ ad h Ad T V λ + S q q + S V λ + S ad h Ad T S S q + ad V q q + S end whle whle backward recurson do F = I V + I V ad V τ end whle + k µ Ad T,k = S T F + S T F Fk q q + S + Ṡ q q + Ṡ I V ad I V V V + I ad h k F k wth gven qa, qa, qa F ext Table 6: Dervatve of the Recursve Forward Dynamcs whle forward recurson do T T h V η = Ad T V λ ad h Ad T = ad S V S q + ad V q q + S end whle whle backward recurson do Î ˆB Ψ Π β Π k V λ + S q q + S + Ṡ q q + Ṡ = I + k µ Ad T Π kad h k Π k ad h k,k = ad V I V ad I V V ext V + I F + k µ Ad βk T ad h k β k,k = Ψ S T Î S + TÎ S S S + TÎ } T S Ψ = Î Ψ Î S ST Î + Î S Ψ S TÎ Î + S Ψ S TÎ } T S + Î Ψ S TÎ S + Î Ψ S T Î = ˆB + Î η + S Ψ τ S T Î η + ˆB } + Î η + S Ψ Ψ + S τ S T Î η + ˆB T } Ad T,k T 13
14 + S Ψ τ S end whle whle forward recurson do q V F = Ψ τ S Î T AdT S τ + Ψ S TÎ Ad T S T ˆB S T ˆB = Ad V λ T = Î end whle V + Î V TÎ η + ˆB S T Î η η + Î + ˆB } V λ + η S T ˆB } TÎ + S T Î AdT V λ } ad h Ad T V λ + η V λ + η ad h Ad T V λ + S q q + S + η + ˆB Table 7: Dervatve of the Recursve Hybrd Dynamcs whle forward recurson do T T h V η = Ad T V λ ad h Ad T = ad S V S q + ad V q q + S end whle whle backward recurson do Π k Î = I + k µ Ad T Π kad h k,k ˆB = ad V I V ad I V V V + I + k µ Ad βk T ad h k β k,k f a then else Π β Ψ Π β = Î = ˆB V λ + S q q + S + Ṡ q q + Ṡ F ext Π k ad h k + Î η + S q + Î η + S q q + S S T Î TÎ S + = Ψ S + S = Î Ψ Î S ST Î + Î S TÎ } T S Ψ T } Ad T,k S Ψ S TÎ Î + S Ψ S TÎ } T S + Î Ψ S TÎ S + Î Ψ S T Î + Î η + S Ψ τ S T Î η + ˆB } + Î η + S Ψ Ψ + S τ S T Î η + ˆB = ˆB + S Ψ τ end f end whle whle forward recurson do f a then S T TÎ η + ˆB S T Î η η + Î + ˆB } 14
15 V = Ad V λ T ad h Ad T V λ + S q q + S + η F V + Î V τ = Î = S T F + S T F + ˆB else q = Ψ τ S Î T AdT τ + Ψ S V F end f end whle S T Î Ad T S T ˆB S T ˆB = Ad V λ T = Î V + Î V TÎ + S T Î V λ } V λ + η S T ˆB } AdT ad h Ad T V λ + η V λ + η ad h Ad T V λ + S q q + S + η + ˆB To obtan the dervatves of the dynamcs, t s needed to run the assocated recursve dynamcs n advance, and n addton, the followng quanttes S, Ṡ ext F,, I, and h 64 for each body should be set properly n the ntalzaton step. The algorthms for the dervatves of the dynamcs are so general that they can be appled to dfferentatng the equatons of moton wth respect to any arbtrary scalar varable. It should be noted that, by mplementng the algorthms cleverly, the calculaton speed can be much faster than ts nave mplementaton, because, n some cases, many of the quanttes n 64 are zero. For example, f p = q k where q k R denotes the k-th coordnate of system, then I = 0 S = Ṡ f p = q k, then S Ṡ = I f p = q k, then Ṡ References = h = h = 0 = 0 when qk / q = S = I = 0 when qk / q = h = 0. [1] Roy Featherstone, Robot Dynamcs Algorthms, Kluver Academc Publshers, [2] Rchard M. Murray, Zexang L and and S. Shankar Sastry, A Mathematcal Introducton to Robotc Manpulaton, CRC Press,
16 [3] F. C. Park, J. E. Bobrow and S. R. Ploen, A Le group formulaton of robot dynamcs, Internatonal Journal of Robotcs Research, vol. 14, no. 6, pp , [4] S. R. Ploen and F. C. Park, Coordnate-nvarant algorthms for robot dynamcs, IEEE Transactons on Robotcs and Automaton, vol. 15, no. 6, pp , [5] Garett A. Sohl and James E. Bobrow, A recursve multbody dynamcs and senstvty algorthm for branched knematc chans, Journal of Dynamc Systems, Measurement, and Control, vol. 123, pp ,
Iterative General Dynamic Model for Serial-Link Manipulators
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
More informationMathematical 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 informationRigid 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 informationPhysics 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 informationPHYS 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 informationLecture 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 informationNMT 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 informationWeek 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 informationMEV442 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 informationInner 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 informationPhysics 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χ 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 informationWeek3, 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 informationClassical 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 informationSo 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 informationLesson 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 informationMoments 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 informationLinear 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 informationA 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 informationSpin-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 information11. 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 informationAn 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 informationCHAPTER 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 informationLecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2
P470 Lecture 6/7 (February 10/1, 014) DIRAC EQUATION The non-relatvstc Schrödnger equaton was obtaned by notng that the Hamltonan H = P (1) m can be transformed nto an operator form wth the substtutons
More informationThe 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 informationWeek 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 informationLAGRANGIAN MECHANICS
LAGRANGIAN MECHANICS Generalzed Coordnates State of system of N partcles (Newtonan vew): PE, KE, Momentum, L calculated from m, r, ṙ Subscrpt covers: 1) partcles N 2) dmensons 2, 3, etc. PE U r = U x 1,
More information1 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 informationPhysics 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 informationPhysics 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 informationPHYS 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 informationELASTIC 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 informationPhysics 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 information12. 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 informationCanonical 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 informationFirst 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 informationPHYS 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 informationcoordinates. 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 informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationNotes on Analytical Dynamics
Notes on Analytcal Dynamcs Jan Peters & Mchael Mstry October 7, 004 Newtonan Mechancs Basc Asssumptons and Newtons Laws Lonely pontmasses wth postve mass Newtons st: Constant velocty v n an nertal frame
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More informationTranslational 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 informationTransfer 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 informationThe 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 informationChapter 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 informationConservation 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 informationSpring 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 informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationTHE 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 informationSalmon: 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 informationBezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0
Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the
More informationDifference 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 informationRotational 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 informationModule 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 informationNMT 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 informationPart 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 informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationNEWTON 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 informationMechanics 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 information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More informationA Tale of Friction Basic Rollercoaster Physics. Fahrenheit Rollercoaster, Hershey, PA max height = 121 ft max speed = 58 mph
A Tale o Frcton Basc Rollercoaster Physcs Fahrenhet Rollercoaster, Hershey, PA max heght = 11 t max speed = 58 mph PLAY PLAY PLAY PLAY Rotatonal Movement Knematcs Smlar to how lnear velocty s dened, angular
More informationTechnical 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 informationPoisson 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 informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationPhysics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1
Physcs 141. Lecture 14. Frank L. H. Wolfs Department of Physcs and Astronomy, Unversty of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Course Informaton: Lab report # 3. Exam # 2. Mult-Partcle
More informationENGI9496 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 informationModeling 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 informationLagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013
Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned
More informationT 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 informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not
More informationA 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 informationGeometry and Screw Theory for Robotics
Tutoral (T9) Geometry and Screw Theory for Robotcs Stefano Stramgol and Herman Bruynnckx March 15, 2001 2 Contents 1 Moton of a Rgd Body 5 1.1 Introducton... 5 1.2 The Eucldean Space... 6 1.2.1 Ponts versus
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationLecture 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 informationTime-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 informationClassical 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 informationThree 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 informationProfessor 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 informationThe Analysis of Coriolis Effect on a Robot Manipulator
Internatonal Journal of Innovatons n Engneerng and echnology (IJIE) he Analyss of Corols Effect on a Robot Manpulator Pratap P homas Assstant Professor Department of Mechancal Engneerng K G Reddy college
More informationStudy 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 informationThe Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices
Internatonal Mathematcal Forum, Vol 11, 2016, no 11, 513-520 HIKARI Ltd, wwwm-hkarcom http://dxdoorg/1012988/mf20166442 The Jacobsthal and Jacobsthal-Lucas Numbers va Square Roots of Matrces Saadet Arslan
More informationAn efficient algorithm for multivariate Maclaurin Newton transformation
Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationCHEMICAL REACTIONS AND DIFFUSION
CHEMICAL REACTIONS AND DIFFUSION A.K.A. NETWORK THERMODYNAMICS BACKGROUND Classcal thermodynamcs descrbes equlbrum states. Non-equlbrum thermodynamcs descrbes steady states. Network thermodynamcs descrbes
More informationChapter 8. Potential Energy and Conservation of Energy
Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal
More informationKinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017
17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationRepresentation theory and quantum mechanics tutorial Representation theory and quantum conservation laws
Representaton theory and quantum mechancs tutoral Representaton theory and quantum conservaton laws Justn Campbell August 1, 2017 1 Generaltes on representaton theory 1.1 Let G GL m (R) be a real algebrac
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationWeek 8: Chapter 9. Linear Momentum. Newton Law and Momentum. Linear Momentum, cont. Conservation of Linear Momentum. Conservation of Momentum, 2
Lnear omentum Week 8: Chapter 9 Lnear omentum and Collsons The lnear momentum of a partcle, or an object that can be modeled as a partcle, of mass m movng wth a velocty v s defned to be the product of
More informationTensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q
For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals
More informationPY2101 Classical Mechanics Dr. Síle Nic Chormaic, Room 215 D Kane Bldg
PY2101 Classcal Mechancs Dr. Síle Nc Chormac, Room 215 D Kane Bldg s.ncchormac@ucc.e Lectures stll some ssues to resolve. Slots shared between PY2101 and PY2104. Hope to have t fnalsed by tomorrow. Mondays
More informationC/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 informationFeb 14: Spatial analysis of data fields
Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s
More information10. 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 informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More information