Accuracy Analysis of General Parallel Manipulators with Joint Clearance

Size: px
Start display at page:

Download "Accuracy Analysis of General Parallel Manipulators with Joint Clearance"

Transcription

1 007 IEEE Internatonal Conference on Robotcs and Automaton Roma, Italy, Aprl 007 WeC6. Accuracy Analyss of General Parallel Manpulators wth Jont Clearance Jan Meng, ongjun Zhang, Tnghua Zhang, Hong Wang, Zexang L Abstract ue to the jont clearance, parallel manpulators always exhbt some poston and orentaton errors at the moble platform. Ths paper ams to provde a systematc framework for the error analyss problem of general parallel mechansms nfluenced by the jont clearance. A novel and effcent method s proposed to evaluate the maxmal pose errors of general spatal parallel manpulators wth jont clearance. I. INTROUCTION espte the fact that jont clearance smplfes the assembly and manufacturng of parallel mechansms, the generated pose error of lnks, however, can not be gnored when the mechansm requres hgh accuracy. Compared wth other error sources, such as assembly and manufacturng errors and motor actuaton errors, etc, jont clearance has more sgnfcant mpact on the poston accuracy for both seral and parallel manpulators. Therefore, t s qute mportant to provde an accurate model that can predct the effects of jont clearance on the mechansm s postonng performance. As ndcated n [1], contrary to the assembly and manufacturng errors, jont clearance leads to uncertan error motons at an arbtrary pose of the mechansm. Its effects are hghly non-repeatable and can not be rectfed wth any knd of calbraton. Much research has been devoted to the problem of accuracy n parallel mechansms. Some authors appled probablstc analyss to determne the pose error of clearanceaffected jonts and movng platform []. Parent-Castell and Venanz [3] used the vrtual work method to determne the poston that the movng platform reaches when a gven external load s exerted on t. P.Voglewede and I.Ebert-Uphoff [4] amed to predct precsely the moble platform error motons for 3-RRR and some specal parallel mechansms wth jont clearance. In ths paper, we wll provde a deeper nsght nto the accuracy problem of parallel mechansms affected by jont clearance. A general and yet effcent error evaluaton method s also proposed whch can handle any knd of parallel mechansms, no matter planar or spatal, overconstraned or non-overconstraned. Furthermore, ts effcency makes t possble to compute the global maxmum pose errors of a clearance-affected mechansm n a prescrbed workspace, other than just at a gven theoretcal confguraton. Fg. 1. Fg.. Pose Error of A Rgd Body Clearance-affected Revolute Jont II. POSE ERROR ANALYSIS OF PARALLEL MANIPULATORS WITH CLEARANCE-AFFECTE PAIRS A. End-effector and Jont Pose Error Fg.1 shows the end-effector of a parallel mechansm wth some confguraton (ncludng the poston and orentaton) error. Suppose that the nomnal (or deal) confguraton of the end-effector (the dashed rgd body n Fg.1) s g 0 SE(3), whch s the relatve transformaton from the nomnal body frame B to the nertal frame A. Wth some confguraton error, the relatve transformaton from the actual body frame B to the nertal frame A, however, gves rse to the real confguraton g of the end-effector. Let A be another frame attached to the end-effector, whch s chosen such that the relatve transformaton from B to A s g 0 (see Fg.1). Then t s clear that wthout confguraton error, A frame wll co- Ths project s supported by RGC Grant No. HKUST6301/03E, HKUST676/04E, and J.Meng,.Zhang and Z.L are wth epartment of Electrcal and Electronc Engneerng, Hong Kong Unversty of Scence and Technology,Clear Water Bay, Kowloon, Hong Kong eejmeng@ust.hk T.Zhang and H.Wang are wth the vson of Control and Mechatronc Engneerng, HIT Shenzhen Graduate School, P.R. Chna Fg. 3. Clearance-affected Prsmatc Jont /07/$ IEEE. 889

2 WeC6. ncde wth the nertal frame A, and B frame wll concde wth B. Hence the error rgd moton of the end-effector from the nomnal confguraton to the real confguraton can be descrbed ether by the relatve transformaton from B to B, whch s g0 1 g, or the relatve transformaton from A to A, that s g g0 1. Correspondngly, the confguraton error can be defned as (g0 1 g e) and (g g0 1 e), whch are called the confguraton error wth respect to the body frame and the nertal frame respectvely. Normally the error rgd motons g0 1 g and g g 1 0 are restrcted wthn a small neghborhood of e. Hence the two expressons of the confguraton error can be reasonably approxmated by ther frst order terms, that s, g0 1 g e = eê1 xbeê ybeê3 zbeê4 αbeê5 βbeê6 γb e (ê 1 x b + +ê 6 γ b ) where {e = 1,, 6} s the canoncal bass of R 6, and the sx components x b,, γ b (and x a,, γ a ) represent the three translaton and three rotaton errors of the end-effector about the x, y and z axs of the B (and A) frame respectvely. In the above two equatons, the operaton dentfes R 6 wth se(3) (see [5] for more detals). Hence f we let δg b =( x b, y b, z b, α b, β b, γ b ) T and δg a =( x a, y a, z a, α a, β a, γ a ) T, then t s easy to see that g0 1 g e δg ˆb and g g0 1 e δg ˆa.To avod the confuson wth the term confguraton error, δg b wll be called the pose error wth respect to body frame, and δg a the pose error wth respect to the ntal frame. From now on, the pose error of the end-effector other than the confguraton error wll be manly used and studed n ths paper. As the relatve transformaton from B to A frame s g 0, the two expressons of the pose error are related by δg a = Ad g0 δg b, where for g =(p, R), the Adjont map s defned by R ˆpR Ad g = 0 R The pose error of the end-effector s caused by many factors, such as assembly and manufacturng errors (occurrng at the adjacences between the lnks and the jonts), actuaton motor errors, jont clearance, etc. As the major source of the pose error comes from the error motons caused by the jont clearance [6], n ths paper, we wll focus on studyng the mpact of jont clearance on the devatons of the endeffector s confguratons. Hence, we wll assume that jont clearance s the unque error source. In the reman part of ths subsecton, we study the mpact of jont clearance on the pose error of two lnks connected by a clearance-affected jont. Only revolute and prsmatc jonts are consdered n ths paper, but other types of jonts can be analyzed n a smlar manner. A clearance-affected knematc par (.e., a par wth jont clearance) actually has 6-oFs. Thus at a gven nomnal confguraton (the theoretcal confguraton of the par), there are 6 twsts assocated wth a par. Some twsts result n the desred motons of the jont, called deal twsts, whereas others gve rse to the error motons caused by jont clearance. For example, for a clearance-affected revolute par n Fg., f we attach a local coordnate frame C to the bearng wth z-axs along the par s theoretcal axs, then e 6 generates the deal rotatonal motons of the jont, whereas e 1,,e 5 generate error motons of the shaft wth respect to the bearng. The pose error of the relatve confguraton of the two lnks wth respect to the local frame C thus s gven by δυ c = 5 e σ (1) where σ 1, σ, σ 3 represent 3-nfntesmal error translatons along x, y, z axs, and σ 4, σ 5 represent - and g g0 1 e = eê1 xaeê yaeê3 zaeê4 αaeê5 βaeê6 γa enfntesmal error rotatons about x, y axs. Hereafter (ê 1 x a + +ê 6 γ a ) =1 σ,=1,, 5 wll be called the error motons caused by the jont clearance, wth values restrcted by the jont geometry and the magntude of jont clearance. For a partcular desgn of a clearance-affected revolute jont, t s possble to formulate explctly the set of constrants that bound the values of σ. For example, for the journal bearng desgn (Fg.) of a revolute jont, because of the axal symmetry, the x, y axs can be chosen freely. Wth orgn chosen at the jont center, and the geometrcal dmensons L, and the magntude of jont clearance ɛ r, ɛ a gven beforehand, the values of σ are constraned by a set of quadratc and second order cone constrants : ( σ 1 + L σ 5) +( σ L σ 4) ɛ r ( σ 1 L σ 5) +( σ + L σ 4) ɛ r σ 4 + σ5 σ () 3 ɛ a σ 4 + σ5 + σ 3 ɛ a For other types of desgns of clearance-affected revolute jont, though the structure may be more complcated, one can get smlar constrant nequaltes as (). If there s an nertal reference frame A such that the relatve transformaton from C to A s g R ac, then the pose error caused by jont clearance can be wrtten wth respect to the nertal frame A by δυ a = Ad g R ac δυ c = 5 =1 (Ad g R ac e ) σ (3) Smlarly, for a clearance-affected prsmatc par(fg.3), f we attach a local coordnate frame C to the supportng carrage wth x-axs along the par s deal translatonal axs, then the pose error of the slder wth respect to the supportng carrage s gven wth respect to C by δγ c = 6 e j τ j (4) j= where τ j,j =,, 6 are also error motons caused by the jont clearance. If y, z axs are parallel to the two faces of the supportng carrage, and the orgn les on the geometrcal center of the jont (see Fg.3), then values of the 890

3 WeC6. error motons are restrcted by a set of lnear constrants ɛ a + L τ 6 τ 4 + τ ɛ a ɛ a L τ 6 + τ 4 + τ ɛ a ɛ b L τ 5 + τ (5) 4 + τ 3 ɛ b ɛ b + L τ 5 τ 4 + τ 3 ɛ b Supposng that the transformaton from C to the nertal frame A s g P ac, then the pose error s gven wth respect to frame A by δγ a = Ad g P ac δγ c = 6 j= (Ad g P ac e j) τ j (6) Fnally, we remark that the two desgns of Fg. and Fg.3 have ever been ntensvely studed n [7] and [3]. However, the constrants of error motons n those lteratures are formulated as sets of general nonlnear nequaltes. In ths paper, we have reformulated them nto sets of so-called convex constrants, as seen from () and (5). As we wll see later, ths reformulaton wll make the convex optmzaton feasble to the evaluaton of the pose errors. Fg. 4. Tsa s 3-UPU Manpulator B. Gross Pose Error of Jonts The desred motons of an actuated jont are controlled by the motors mounted on the jont. Assumng there s no motor actuaton error, then at a nomnal confguraton, the pose error of an actuated jont s solely rased by ts jont clearance. Hence the expresson of gross pose error of an actuated jont s the same as Eq.(1) and (3). However, for passve jonts n a parallel mechansm, n addton to the errors motons caused by jont clearance, the dle motons about the deal twsts also play a role n the pose errors of ther two lnks. Hence, the gross pose error of a passve revolute jont s gven wth respect to ts local frame by δυ c = 5 e σ + e 6 θ (7) =1 where θ s the dle rotaton error of the shaft about the jont s deal axs. Wth respect to the nertal frame A, the gross pose error s δυ a = Ad g R ac δυ c = 5 =1 (Ad gac R e ) σ +(Ad g R ac e 6 ) θ = 5 =1 (Ad gac R e ) σ + ξ θ where ξ s the deal twst of the jont when expressed n the nertal frame. Smlarly, two gross pose errors of a passve prsmatc jont can be wrtten as Fg. 5. A Parallelogram Jont C. End-effector Pose Error of A Parallel Mechansm In ths subsecton, we study the pose error of the endeffector of a parallel mechansm caused by the clearance of ts consttuent jonts. We use some examples to llustrate how to get the pose error of the end-effector at a general theoretcal confguraton(say, g 0 ), subject to the error motons generated by all the clearance-affected jonts. Example 1: Pose Error Analyss of Tsa s Manpulator Consder a Tsa s 3-UPU manpulator n Fg.4. Each U- jont conssts of two perpendcularly connected clearanceaffected revolute jonts modelled by Fg.. The three prsmatc jonts P 1,,P 3, beng actuated, are also clearanceaffected and modelled by Fg.3. The nertal frame A and the body frame B are located at the center of the base and the moble platform respectvely. And the local frames of R j (denoted Cj R), and that of P (denoted by C P ), =1,, 3, j =1,, 4, are located at the center of each jont, wth x, and 6 δγ c = e j τ j + e 1 µ (8) j= δγ a = Ad g P ac δγ c = 6 j= (Ad gac P e j) τ j + η µ respectvely, where µ s the dle translatonal error along the jont s deal axs, and η s the deal twst representaton of the jont n the nertal frame A. Fg. 6. Orthoglde Manpulator 891

4 WeC6. y axs chosen at one s convenence. The moble platform(or end-effector) deally can undergo three of pure translatons T (3) under some geometrcal condtons [8]. At a nomnal confguraton g 0 T (3) whch s not sngular for the sake of smplcty, the transformaton matrx from Cj R and CP to A, denoted by gj R and gp respectvely, can be obtaned through the nverse knematc analyss of the deal mechansm at g 0. From the last subsecton, the pose error of P s gven wth respect to the nertal frame A by P = 6 l= (Ad g P e l) τ l =1,, 3 and the pose error of R j s gven by R j = 5 l=1 (Ad gj R e l) σ jl + ξ j θ j where =1,, 3,j =1,, 4, and ξ j s the deal twst assocated wth R j. The pose error of th open subchan, =1,, 3, can thus be wrtten by E = P + 4 j=1 R j = 6 l= (Ad g P e l) τ l + 4 j=1 ξ j Θ j + 4 j=1 5 l=1 (Ad g R j e l) σ jl = A dγ + 4 j=1 B j dυ j + J dθ = A dγ + B dυ + J dθ where A = [Ad g P e Ad g P e 6 ], dγ = ( τ,, τ 6 ) T, B j = [Ad g R j e 1 Ad g R j e 5 ], dυ j = ( σ j1,, σ j5 ) T, B = [B 1 B 4 ], dυ = (dυ T 1,,dΥT 4 )T, J = [ξ 1 ξ 4 ], dθ =( θ 1,, θ 4 ) T. Equatng E 1 and E, E and E 3, we can get J dθ 1 dθ dθ 3 = A dγ 1 dγ dγ 3 + B dυ 1 dυ dυ 3 (9) (10) where J1 J J = (11) J J 3 s a nonsngular square matrx, and A1 A A = B1 B, B = A A 3 B B 3 From Eq.(10) and (11), we may see that the dle passve jont motons are unquely determned by the error motons τ l and σ jl caused by the jont clearance. After the dle motons dθ are determned, we can substtute t nto Eq.(9) for any subchan to obtan the pose error of the movng platform the nomnal confguraton g 0. Choosng the frst subchan as example, we can wrte the result n a concse way as follows E g0 = A 1 dγ 1 + B 1 dυ 1 + J 1 dθ 1 = A 1 dγ 1 + B 1 dυ 1 +[J 1 00] dθ 1 dθ dθ 3 = A 1 dγ 1 + B 1 dυ 1 +[J 1 00] J 1 dγ 1 dυ 1 A dγ + B dυ dγ 3 dυ 3 = H (dγ T 1,dΓ T,dΓ T 3,dΥ T 1,dΥ T,dΥ T 3 ) T In the above equaton, we get a lnear transformaton matrx H that maps the error motons of the clearanceaffected jonts, τ l and σ jl, to the pose error of the movng platform E g0. From Eq.(10), τ l and σ jl are free varables of the mechansm, the number of whch reflects the total extra of of the mechansm at a non-sngular nomnal confguraton. However, the values of τ l and σ jl are bounded wthn the sets of convex constrants gven by nequaltes () and (5) for each jont. Therefore, the maxmal pose error of the movng platform at the nomnal confguraton g 0, denoted by Eg m 0, can be obtaned by sx convex optmzatons appled to each components of E g0. If one further want to get the maxmal pose error of the movng platform n a prescrbed workspace W, for example, a cube centered by the home confguraton, he can dscretze W and calculate Eg m at any pont of W. After that, the global maxmal pose error n W can be obtaned by E m = (max g W Em g (1),, max g W Em g (6)) T where E m g () s the th component of E m g. For overconstraned parallel manpulators, however, not all error motons assocated wth jonts clearance are free varables, as seen from the next example. Example : Pose Error Analyss of A Parallelogram Jont In ths example, we study the pose error of a parallelogram (P a ) composng of four clearance-affected revolute jonts. Snce P a usually serves as an extended par n complex parallel mechansms, e.g., the elta robot, we also attach a frame C, called the local frame of the parallelogram jont, to the base lnk of the parallelogram, as shown n Fg.5. Frame C s located at the center of the base lnk, wth z axs parallel to the deal rotaton axs of R j, and x axs pontng to the geometrcal center of R 1. Thus the deal confguraton space of the parallelogram C M sacrclenthex-y plane of C. At each nomnal confguraton g 0 C M, we defne four local frames C j for R j, =1,, j =1,. Asthe choce of x, y axs of C j s free, for our convenence, C j may be chosen wth orgn lyng at the geometrcal center of R j, and x, y, z axs algned wth that of C. Hence, the transformaton from C j to C s gven by I pj g j = 0 1 where p j s the geometrcal center of R j. Assume that R 11 s the actuaton jont, then the pose error of four revolute jonts wth respect to the C frame s gven by R j = 5 (Ad gj e l ) σ jl + ξ j θ j =1,,j =1, l=1 where θ 11 =0. The pose error of th subchan, =1, thus can be easly got as E = R 1 + R = j=1 5 l=1 (Ad g j e l) σ jl + j=1 ξ j Θ j = j=1 A j dυ j + j=1 ξ j Θ j 89

5 WeC6. where A j = [Ad gj e 1 Ad gj e 5 ], dυ j = ( σ j1,, σ j5 ) T, = 1,, j = 1,. As θ 11 = 0, by equatng E 1 and E, we can get J Θ 1 Θ 1 = A Θ dυ 1 dυ 1 dυ (1) where J = [ξ 1 ξ 1 ξ ], and A = [ A 11 A 1 A 1 A ]. Clearly, J s a 6 3 full column rank matrx. J may combne wth three twsts Ad g1 e 3, Ad g1 e 4, and Ad g1 e 5 to form a 6 6 non-sngular square matrx J =[ξ 1 ξ 1 ξ Ad g1 e 3 Ad g1 e 4 Ad g1 e 5 ]. Then Eq.(1) can be re-wrtten as follows J Θ 1 Θ 1 Θ σ 13 σ 14 σ 15 = A dυ 1 σ 11 σ 1 dυ (13) where A =[ A 11 A 1 Ad g1 e 1 Ad g1 e A ]. Therefore, only 17 error motons, dυ 1, dυ and σ 11, σ 1 are free varables of the mechansm. The other three jont error motons σ 13, σ 14, σ 15, plus the dle passve jont motons Θ 1, Θ 1, Θ, are unquely determned by these free varables. enote by B = J 1 A, B() the th row of B, =1,, 6, and g ac the transformaton from the local frame C to the nertal frame A. Then pose error of the end-effector at the nomnal confguraton g 0 s gven wth respect to the nertal frame by E g0 = Ad gac ( R 11 + R ( 1 ) ) dυ11 = Ad gac [A 11 A 1 ] + Ad dυ gac ξ 1 Θ 1 ( 1 ) dυ11 = Ad gac [A 11 A 1 ] dυ 1 dυ 1 +Ad gac ξ 1 B(1) σ 11 σ 1 dυ = H (dυ T 11 dυ T 1 σ 11 σ 1 dυ T ) T subject to the constrants ( σ j1 + L σ j5) +( σ j L σ j4) ɛ r ( σ j1 L σ j5) +( σ j + L σ j4) ɛ r σj4 + σ j5 σ j3 ɛ a σj4 + σ j5 + σ j3 ɛ a σ 13 B(4) dυ 1 σ 14 = B(5) σ 11 σ 15 B(6) σ 1 dυ (15) where =1, and j =1,. Thus the calculaton of maxmal pose error of an overconstraned parallel manpulator at g 0 s a convex optmzaton problem wth lnear equalty constrants. By [9], such a problem can also be quckly solved by very effcent algorthms. P a s often used as an extended passve par n some common parallel mechansms, e.g, the elta robot. In ths case, R 11 s a passve revolute jont nstead of an actuated one, whch mples that θ Assume that the body frame B orgnally concdes wth A frame, then the deal confguraton space of the end-effector s gven by { } I (eˆωθ 11 I)v Q = θ [0, π] where ω s the deal rotaton axs of the parallelogram, and v the vector connectng R 1 and R at the home confguraton. Clearly, at the nomnal confguraton g 0 = Q(θ 11 ), the confguraton error of the end-effector caused by θ 11 s gven wth respect to the nertal frame A by F g0 = dq Q 1 0 ωeˆωθ 11 v θ 11 = θ dθ The total pose error thus s the lnear summaton of F g0 and E g0 (16) T g0 = E g0 + F g0 = H (dυ T 11 dυ T 1 σ 11 σ 1 dυ T ) T +ωeˆωθ11 v θ 11 subject to the constrants of (15). The value of θ 11,as we wll see n the next example, s determned by the loop closure equaton of the parallel mechansm that contans ths P a. Example 3: Pose Error Analyss of Orthoglde The Orthoglde manpulator [10], as shown n Fg.6, conssts of three P RP a R dentcal subchans. The actuated jonts are the three orthogonal prsmatc ones. Ideally, the manpulator s an overconstraned one, wth end-effector moton a subset of T (3). At a nomnal confguraton g 0 T (3), usngthe same approach as before, we may fnd the pose errors of (14) prsmatc jonts P and revolute jonts R j, =1,, 3, j =1,. As the end-effector s connected to the prsmatc jonts through a set of three passve parallelograms, the pose error of the parallelogram s gven by Eq.(16). Thus the pose error of the th subchan, =1,, 3, s shown to be E = P + R 1 + T + R = 6 l= (Ad g P e l) τ l + 5 l=1 (Ad g1 R e l) σ 1l + ξ 1 Θ 1 +H (dυ T 11 dυt 1 σ 11 σ 1 dυ T )T +ω eˆωθ11 v θ l=1 (Ad g R e l) σ l + ξ Θ = A dγ + j=1 B j dω j + H dυ +J dθ = A dγ + B dω + H dυ + J dθ where g P s the transformaton matrx from the local frame of the P to the nertal frame A, and gj R the one from the local frame of R j to A. Furthermore, A = [Ad g P e Ad g P e 6 ], dγ =( τ,, τ 6 ) T, B j = [Ad g R j e 1 Ad g R j e 5 ], dω j = ( σ j1,, σ j5 ) T, 893

6 WeC6. B = [B 1 B ], dω = (dω T 1,dΩT )T, dυ = (dυ T 11 dυ T 1 σ 11 σ 1 dυ T )T, J = [ξ 1 ω eˆωθ11 v ξ ], dθ = ( θ 1, θ 11, θ ) T, =1,, 3, j =1,. Equatng E 1 and E, E and E 3, we can get J dθ 1 dθ dθ 3 = A dγ 1 dγ dγ 3 + B dω 1 dω dω 3 + C dυ 1 dυ dυ 3 where J1 J J = (17) J J 3 s a 1 9 full column rank matrx, and A1 A A = B1 B, B = A A [ 3 ] B B 3 H1 H C = H H 3 Note that the rotaton axs of the twst Ad g P e 4 s along the deal translaton axs of P, =1,, 3, hence the former equaton may be re-wrtten as dθ REFERENCES 1 τ 14 J dθ τ 4 = A dγ 1 dγ + B dω 1 dω + C dυ [1] K.L. Tng, J. Zhu, and. Watkns. The effects of jont clearance 1 on poston and orentaton devaton of lnkages and manpulators. dυ Mechansm and Machne Theory, 35: , 000. dθ 3 dγ 3 dω 3 dυ [] M.Mayouran and J.Rastegar. Stochastc modelng of the mechancal 3 behavor of mechansms n the presence of jont clearances. In τ Proceedngs of 1990 ASME Mechansms Conference, pages , where [3] S. Venanz. Methods for Clearance Influence Analyss n Planar and [ J J1 Ad = g P e ] Spatal Mechansms. Ph thess, Unversty of Bologna, Bologna, 4 J 1 Ad g P e 4 J Ad g P e 4 J 3 Ad g P 3 e 4 s a 1 1 non-sngular matrx, and A A = 1 A A A, 3 A = [Ad g P e Ad g P e 3 Ad g P e 5 Ad g P e 6 ], dγ = ( τ, τ 3, τ 5, τ 6 ) T, =1,, 3. From the above transformed equaton, one may see that only dγ, dω and dυ, = 1,, 3, are free varables of the mechansm. Three jont error motons τ 4, together wth the dle passve jont motons dθ, =1,, 3, are unquely determned by totally =78 free varables. The pose error of the end-effector thus s gven by E g0 = A 1 dγ 1 + B 1 dω 1 + H 1 dυ 1 + J 1 dθ 1 = A 1 dγ 1 + B 1 dω 1 + H 1 dυ 1 +[J 1 0 0] dγ J 1 A 1 dω 1 dυ 1 dγ + B dω + C dυ dγ 3 dω 3 dυ 3 subject to the constrant of () for revolute jonts, (5) for prsmatc jonts, (15) for parallelogram, and the followng lnear equaltes τ 14 τ 4 = (dγ 1T,,dΓ T 3,dΩ T 1,,dΩ T 3,dΥ T 1,,dΥ τ T 3 ) T 34 where the matrx can be derved from Eq.(3). Hence the calculaton of the maxmal pose error of the Orthoglde manpulator at g 0 s a stll a convex optmzaton problem wth lnear equalty constrants. III. CONCLUSION In ths paper, we propose a general method to evaluate the pose(poston and orentaton) error of the end-effectors of parallel manpulators due to the jont clearance. We show that for non-overconstraned parallel manpulators, the error motons caused by jont clearance are free varables subject to some constrants defned by the jont geometry and magntude of clearance. However, for overconstraned parallel mechansms, part of these jont error motons are dependent on the remanng ones. For some partcular desgns of common lower pars, the pose error analyss for both non-overconstraned and overconstraned parallel manpulators can be formulated nto standard convex optmzaton problems, whch makes t possble to compute the maxmum pose errors n a prescrbed workspace other than just at a sngle confguraton. Italy, 004. [4] P. Voglewede and I.E. Uphoff. Applcaton of workspace generaton technques to determne the unconstraned moton of parallel manpulators. ASME Journal of Mechancal esgn, 16:83 90, 004. [5] R. Murray, Z.X. L, and S. Sastry. A Mathematcal Introducton to Robotc Manpulaton. CRC Press, [6] C.H. Han, J.W. Km, J.W. Km, and F.C. Park. Knematc senstvty analyss of the 3-UPU parallel mechansm. Mechansm and Machne Theory, 37: , 00. [7] H.H.S. Wang and B. Roth. Poston errors due to clearance n journal bearng. ASME Journal of mechansms, Transmssons, and Automaton n esgn, 111:315 30, [8] L.W. Tsa. Knematcs of a three-of platform wth three extensble lmbs. n Recent Advances n Robot Knematcs: Analyss and Control (Lenarc J., Husty M.L. Eds.), Kluwer, pages 49 58, [9] S. Boyd and L. Vandenberghe. Convex Optmzaton. Cambrdge Unversty Press, 1st edton, 001. [10].Chablat and P. Wenger. Archtecture optmzaton of a 3-dof translatonal parallel mechansm for machnng applcatons, the orthoglde. IEEE Transacton on Robotcs and Automaton, 19(3): ,

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

Iterative General Dynamic Model for Serial-Link Manipulators

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 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

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot Sensors & Transducers 214 by IFSA Publshng, S. L. http://www.sensorsportal.com Desgn and Analyss of Landng Gear Mechanc Structure for the Mne Rescue Carrer Robot We Juan, Wu Ja-Long X an Unversty of Scence

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

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

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

χ 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

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

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

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

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX 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 information

Development of a Novel 3-DoF Purely Translational Parallel Mechanism

Development of a Novel 3-DoF Purely Translational Parallel Mechanism 27 IEEE Internatonal Conference on Robotcs and Automaton Roma Ital 1-14 Aprl 27 Development of a Novel -DoF Purel Translatonal Parallel Mechansm Yunjang Lou Jangang L Jnbo Sh and Zeang L Abstract In vew

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

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

Advanced Mechanical Elements

Advanced Mechanical Elements May 3, 08 Advanced Mechancal Elements (Lecture 7) Knematc analyss and moton control of underactuated mechansms wth elastc elements - Moton control of underactuated mechansms constraned by elastc elements

More information

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI 2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The 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 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

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

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

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

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

On a direct solver for linear least squares problems

On a direct solver for linear least squares problems ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance 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 information

The Analysis of Coriolis Effect on a Robot Manipulator

The 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 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

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

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

A New Concept of Modular Parallel Mechanism for Machining Applications

A New Concept of Modular Parallel Mechanism for Machining Applications Submtted to ICRA A New Concept of Modular arallel Mechansm for Machnng Applcatons Damen Chablat and hlppe Wenger Insttut de Recherche en Communcatons et Cybernétque de Nantes, rue de la Noë, 44 Nantes,

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

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

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

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

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. 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

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

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

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

arxiv:quant-ph/ Jul 2002

arxiv:quant-ph/ Jul 2002 Lnear optcs mplementaton of general two-photon proectve measurement Andrze Grudka* and Anton Wóck** Faculty of Physcs, Adam Mckewcz Unversty, arxv:quant-ph/ 9 Jul PXOWRZVNDR]QDRODQG Abstract We wll present

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

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

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Supplemental document

Supplemental document Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

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

Dimensional Synthesis of Wristed Binary Hands

Dimensional Synthesis of Wristed Binary Hands Dmensonal Synthess of Wrsted Bnary Hands Neda Hassanzadeh Department of Mechancal Engneerng, Idaho State Unversty 921 S. 8th. Ave, Pocatello, ID 83209, USA e-mal:hassneda@su.edu Alba Perez-Graca Assocate

More information

9 Characteristic classes

9 Characteristic classes THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct

More information

A new construction of 3-separable matrices via an improved decoding of Macula s construction

A new construction of 3-separable matrices via an improved decoding of Macula s construction Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

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

Supporting Information

Supporting Information Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to

More information

Formulas for the Determinant

Formulas 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 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

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

POSITION ANALYSIS OF THE RRP-3(SS) MULTI-LOOP SPATIAL STRUCTURE (*)

POSITION ANALYSIS OF THE RRP-3(SS) MULTI-LOOP SPATIAL STRUCTURE (*) Author manuscrpt, publshed n "Journal of Mechancal Desgn 18, 1 (006) 7-78" aper No. MD-04-1161 OITION ANALYI OF THE RR-3() MLTI-LOO ATIAL TRCTRE (*) Carlo Innocent Department of Mechancal and Cvl Engneerng

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Nonlinear Effects in Stiffness Modeling of Robotic Manipulators

Nonlinear Effects in Stiffness Modeling of Robotic Manipulators A. Pashkevch, A. Klmchk D. Chablatonlnear effect n stffness modelng of robotc manpulators, In: Proceedngs of Internatonal Conference on Computer and Automaton echnology (ICCA 29), Vence, Italy, 29, World

More information

APPENDIX A Some Linear Algebra

APPENDIX 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 information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

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

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

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The 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 information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

A Multi-Axis Force Measurement System for a Space Docking Mechanism

A Multi-Axis Force Measurement System for a Space Docking Mechanism 3rd Internatonal Conference on Materal, Mechancal and Manufacturng Engneerng (IC3ME 215) A Mult-Axs orce Measurement System for a Space Dockng Mechansm Gangfeng Lu a*, Changle L b and Zenghu Xe c Buldng

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES Abdelkader Benchou, PhD Canddate Nasreddne Benmoussa, PhD Kherreddne Ghaffour, PhD Unversty of Tlemcen/Unt of Materals

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

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY

ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY CIPRIAN ACATRINEI Natonal Insttute of Nuclear Physcs and Engneerng P.O. Box MG-6, 07725-Bucharest, Romana E-mal: acatrne@theory.npne.ro. Receved March 6, 2008

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

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

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

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

Inexact Newton Methods for Inverse Eigenvalue Problems

Inexact Newton Methods for Inverse Eigenvalue Problems Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.

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

Linear Approximation with Regularization and Moving Least Squares

Linear 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 information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

SL n (F ) Equals its Own Derived Group

SL n (F ) Equals its Own Derived Group Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

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

COS 521: Advanced Algorithms Game Theory and Linear Programming

COS 521: Advanced Algorithms Game Theory and Linear Programming COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

5 The Rational Canonical Form

5 The Rational Canonical Form 5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors. SCALARS AND ECTORS All phscal uanttes n engneerng mechancs are measured usng ether scalars or vectors. Scalar. A scalar s an postve or negatve phscal uantt that can be completel specfed b ts magntude.

More information

2.3 Nilpotent endomorphisms

2.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 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

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge

More information

CHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS COURSE CODES: FFR 35, FIM 72 GU, PhD Tme: Place: Teachers: Allowed materal: Not allowed: January 2, 28, at 8 3 2 3 SB

More information

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem Appled Mathematcal Scences Vol 5 0 no 65 3 33 Interactve B-Level Mult-Objectve Integer Non-lnear Programmng Problem O E Emam Department of Informaton Systems aculty of Computer Scence and nformaton Helwan

More information

On the Multicriteria Integer Network Flow Problem

On the Multicriteria Integer Network Flow Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of

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

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information