Modelling and Analysis of Planar Robotic Arm Dynamics Based on An Improved Transfer Matrix Method for Multi-body Systems

Size: px
Start display at page:

Download "Modelling and Analysis of Planar Robotic Arm Dynamics Based on An Improved Transfer Matrix Method for Multi-body Systems"

Transcription

1 he 4th FoMM World ongress, ape, awan, ctober 5-3, 5 D Number:.6567/FoMM.4H.W.S3.46 Modellng and Analyss of Planar Robotc Arm Dynamcs Based on An mproved ransfer Matrx Method for Mult-body Systems W. hen H.L X.Z.Zhang anjn Key Laboratory of the Desgn and ntellgent Academy of Mltary Medcal Scences ontrol of the Advanced Mechatroncal System anjn, hna anjn, hna Abstract: An mproved transfer matrx method for mult-rgd-body systems s presented based on the structural characterstcs of the manpulator and the relatons between moton transfer and force transfer. hus, a new model applcable to seral robots s obtaned by usng of dscrete tme transfer matrx method. Based on the mproved method, a unfed model for a manpulator wth n DF s derved. hs model contans both knematc and dynamc characterstcs of the system. Accordng to the moton transformaton relaton between the jont space and the operatonal space of the manpulator, the soluton and the procedure are developed for both forward and nverse knematcs of the unfed model. he smulaton results of a DF planar manpulator show the rghtness of the modellng method and exhbt ts superorty to other methods. Keywords: Mult-body system, Planar manpulator, An mproved transfer matrx method, Dynamcs modelng, Forward knematcs, nverse knematcs. ntroducton Wth the rapd progress of scence and technology, robots are frequently beng appled n areas such as ndustral automaton, mltary applcatons, semconductor manufacture and space exploraton. Seral jont robots are constructed by smulatng the human forearm. her advantages nclude flexblty of moton, a compact structure and hgh versatlty. hey are the most commonly used structures n robots, ncludng well-known products such as the Programmable Unversal Machne for Assembly (PUMA) and the Selectve omplance Artculated Robot Arm (SARA) []. A seral jont robot s a complex system wth multple degrees of freedom, multple varables, a hgh degree of non-lnearty and mult-parameter couplng. he knematcs and dynamcs modellng and analyss of a robot are the foundaton of the study. he Denavt-Hartenberg (D-H) method s a general method for the study of robot knematcs[]. Based on the D-H method, Yang W.L.[3] studed the knematcs modellng of a notched contnuum manpulator ncludng the mechancs-based forward knematcs and the curve-fttng-based nverse knematcs. n order to establsh the forward knematcs model by usng D-H procedure, the complant contnuum manpulator featurng the hyper-redundant degrees of freedom s smplfed nto fnte dscrete jonts. Al-Mashhadany Y.. [4] presents a cogntve archtecture for soluton of nverse knematcs problem (KP) of 6-DF elbow manpulator wchen@tjut.edu.cn, lhao4_989@6.com, 56787@sohu.com wth sphercal wrst by Locally Recurrent Neural Networks (LRNNs) and smulated the soluton by usng MALAB/Smulnk. Whle common dynamcs modellng methods for manpulators nclude the Lagrange method, the Newton-Euler equaton and the Kane equaton. Korayem M.H.[5] studed the path plannng dynamc model of three jonts weghtlftng robot usng Lagrange method. XE L.M and HEN L.[6] analyed the knematcs and dynamcs of the system by use of momentum conservaton. Wth the dea of an augmentaton approach, the dynamc equatons can be lnearly dependent on a group of nertal parameters. Allan J.F.[7] studed the geometrcal model of the new manpulator, an analytcal soluton for the nverse knematc equatons, the dynamc model used to compute torques/forces at the actuators, and smulatons performed wth MALAB and AA. Yaskevch N. [8] presented a new recursve algorthm based on Kane equatons matrx form for closed form dynamcs equatons calculaton for smple knematcal chan. However, n the studes of mathematcal models of mult-jont seral robotc arms, almost all methods are amed at system knematcs or dynamcs separately. And these two types of the models can hardly be unfed n the descrpton of one mathematcal model. n fact, the knematcs and dynamcs usually are nterrelated each other wthn the study of a manpulator. n ths paper, an mproved dscrete tme transfer matrx method for mult-body systems s frst presented and a mathematcal model establshed of a mult-jont seral robotc arm n the planar moton accordng to ts structural characterstcs. he equatons thus developed descrbe both the knematcs and the dynamcs. Next, based on ths unfed model, and by settng approprate boundary condtons, a soluton method s developed. Fnally, numercal smulaton s used to valdate the effectveness of the method.. he mproved Dscrete me ransfer Matrx Method for Mult-body Systems he dscrete tme transfer matrx method for mult-body systems s a knd of dynamcs method based on the classc transfer matrx method. t s applcable to general mult-body systems [9]. ts modellng prncple s to perform the lnearaton and the dscretaton of the dynamcs equatons of the components n the tme doman by usng a numercal ntegraton method. Moreover, the prncple assembles of the mult-body system transfer equatons and transfer matrces to calculate system moton by gradual ntegraton. herefore ths approach has the characterstcs of low matrx order, low computaton complexty, hgh

2 modellng flexblty, etc. A mult-jont seral manpulator s a typcal mult-body system, whch has features of tme varablty, non-lnearty and wde range of moton. n ths paper, a unfed knematcs and dynamcs model of a planar manpulator wth n rotaton jonts s developed usng the mproved dscrete tme transfer matrx method for mult-body systems. he modellng theory of ths method s explaned brefly n the followng secton. A. Modellng theory Accordng to the modellng dea of the dscrete tme transfer matrx method for mult-body systems, the state vector descrbng the planar moton of a rgd body s defned as[9] x, y,, m, q x, q,] () [ y where x and y are poston s of the connectng pont wth respect to the nertal reference system, s the angular dsplacement of the rgd body rotatng around -axs, m s the of the connectng pont nternal moment at the -axs, and q and q are s of x nternal force at the connectng pont n the nertal system. y o y o ( ) y θ Fg.. A Rgd Body wth Planar Moton As shown n Fg., a rgd body wth planar moton has one nput end and one output end; the nput end s, output end s, and the centre of mass s. n the dagram, the nerta system s oxy; n the body-fxed system x y, whose orgn s fxed at nput end, the of the output end s (b, b ), the s of the centre of mass are ( c, c ). Because the system s rgd body, the body-fxed system o x y s obtaned by rotatng the translatonal system θ degrees as follows[9] x y x y x y x y x y y () x y x y x y c c c c c c b c b s c s c s c b s b c where s sn and c cos. When the above method s used to deduce the transfer x c x x (3) (4) relaton between the state vectors of moton for the nput end and output end of the rgd body, the numercal ntegraton method should be determned frst to perform dscretaton and lnearaton on the non-lnear terms. ommon numercal ntegraton methods nclude Euler s method, the Newmark-β method, the Wlson-θ method, and the Houbolt method. n ths study, the Newmark-β method s used to perform the dscretaton on the non-lnear functons, such as sn θ and cos θ, n the model [9]. Because the effects of an external force and an external moment exerted on a rgd body are equvalent to an external force and a moment exerted at the rgd body mass centre, the rgd body s treated as f the external forces and the external moments are all exerted at the centre of mass. he moton equaton for the centre of mass s as follows mx my q q x, y, q q x, y, f f x, y, where m s the mass of the rgd body, x and y are the poston s of the centre of mass relatve to the nertal system, q x, and q y, are the nput end nternal forces, q x, and q y, are the output end nternal forces, and f x, and f y, are the external forces exerted at the centre of mass of the rgd body. he moment of momentum theorem s appled to the rgd body as follows J mx y my x m m m qx, y qy, x fx, y f y, x where J s the absolute moment of momentum of the rgd body relatve to the orgn of the body-fxed system ; J s the rotaton nerta of the rgd body relatve to ; and s the absolute angular acceleraton of the rgd body. Equaton (4) s lneared, equaton (3) s substtuted n equaton (5), and then lnearaton s performed on the temporary equaton. he temporary equaton s substtuted nto equaton (6). he lneared result s expressed n matrx format, resultng n the transfer equaton of a rgd body wth one nput end and one output end (5) (6) U (7) where and are state vectors, as defned n equaton (), and the transfer matrx U s [9] U u4 ma u 4 ma y x my mx ( t ( t u 43 ( t ) ) ( t ) ) u 45 u 46 bg bg b G bg u47 u 57 u67 (8) he above development yelds the relatons between rgd bodes wth planar moton. he matrx n equaton (8)

3 descrbes the geometrcal relaton, the moton relaton and the force relaton between rgd bodes [9]. he frst three lnes (-3) n the matrx descrbe the knematcs relatons of the rgd bodes wth planar moton, and the next three lnes (4-6) represent the force and moment stuatons durng the rgd body moton,.e., the rgd body dynamcs. Y 5 n- j n n+ B. he mproved dscrete tme transfer matrx method n the case of a mult-jont seral planar manpulator, the motor provdng the drvng force s located at jont, causng the angular dsplacements before and after jont + to dffer, as shown n Fg.. f the angular dsplacement of jont s θ, then the angular dsplacement of jont + adds θ' on the bass of the angular dsplacement θ of the prevous jont. herefore, the angle θ' s the rotaton angle of jont + relatve to jont. Hence, when calculatng the angular dsplacement of the transfer matrx, θ +=θ +θ' should be satsfed. θ + + θ' θ + Fg.. Dagram of Jont Moton Angular Dsplacement he exact dfference n ths angular dsplacement results n dfferent values of jont angular dsplacements for the nput end state vector of the (+)-th rgd body and the output end state vector of the -th rgd body. n the transfer matrx, ths relaton can be descrbed as where U U.. he mproved Dscrete me ransfer Matrx Method for Mult-body Systems f a planar robotc arm wth n rotatng jonts s consdered as the study subject, and the jonts and lnks of the manpulator are treated as rgd bodes n plane moton wth one nput end and one output end, then the entre robot s a chaned system composed of n rgd components, as shown n Fg. 3. n the dagram,, 3,, n- denote the jonts of the arm,, 4,, n represent the lnks of the arm, and n+ s end effector of the arm. (9) 3 4 Fg. 3. An n-jont Seral Planar Robotc Arm A. Jont model Under normal crcumstances, the rotatng jont of a arm s drven by a motor. n the model, the jont s smplfed as a central rgd body, whose nput end s the centre poston and the output end s the entre crcumference of the central rgd body. he central rgd body has rotaton moton only,.e., no translatonal moton. Accordng to the descrpton n secton., every jont experences a sudden change n angular dsplacement compared wth the prevous lnk connected to t; thus the calculated transfer matrx s multpled by U shown n equaton (9), n the front, and the transfer equaton s where j,3,,n. j U U () j, j j j j, j B. Lnk model he lnk of the robotc arm s a rgd body wth planar moton. n the model, t s a homogeneous lnk that s fxed to the jont and moves wth t. ts nput end, closer to the base, s the jont s output end; and ts output end, s the end farther from the base. A lnk of the robot experences rotaton moton around the jont and transport-nduced translatonal moton, resultng n a transfer equaton as follows where j,4,, n. U () j, j j j, j. verall model of the robotc arm Accordng to above analyss, a seral robotc arm wth the planar moton s composed of n jonts and n lnks. When a jont and a lnk are vewed as ndependent rgd bodes, the entre system comprses n rgd bodes. Because these rgd bodes are all connected n seres, accordng to the mproved transfer matrx method of mult-rgd-body systems, the overall transfer matrx of the robot s the product of the transfer matrces of n rgd bodes,.e., the unfed knematcs and dynamcs model of the entre robotc arm s where he n,n Uall, () U all s the overall transfer matrx of the system. can be expanded as U all X

4 U all U U UU U UU U U (3) n n n n 4 3 U where U k ( k ~ n ) s determned by equaton (8) and U k ( k ~ n) s determned by equaton (9). As descrbed above, equatons () and (3) compose the n-jont unfed knematc and dynamc model for a seral robotc arm. V. Soluton of the Model n the study of robot knematcs, f the jont varables are known, the computaton of poston and atttude of the end effector s the forward knematcs; conversely, for nverse knematcs, the end effector pose s known and the jont varables are the soluton. Whle n the study of dynamcs, f the external forces are known, the soluton for system parameters s forward dynamcs; conversely, when external forces are to be determned, t s the nverse dynamcs. n the case of the manpulator, the task of nverse dynamcs analyss s to determne the drvng moment to generate a certan movement, whch s mportant both for the dynamcs and the control of robots[]. n the case of the mult-jont seral arm shown n Fg., knematcs study mples two possble stuatons: one s to obtan the nformaton of the end effector by usng known jont moton law,.e., the forward knematc soluton; and the other s to obtan the jont varables from the known moton of the end effector, whch s an nverse knematc soluton. Usng the mult-body system dscrete tme transfer matrx modellng method n any of the above stuatons requres boundary condtons for equaton () to be set accordng to known and unknown condtons. hen, through transfer matrx teraton, unknown parameters for all the state vectors are determned. hese parameters contan all the poston and force nformaton for the jonts and lnks of the arm; thus, the soluton procedure s completely dfferent from the soluton for pure knematcs. Solutons for both stuatons above mentoned contan the nverse dynamcs problem analyss,.e., the combnaton of the knematcs problem and the nverse dynamcs problem. he two dfferent stuatons are descrbed n the followng. A. Soluton for forward knematcs usng the model f the moton law s gven for each jont of the manpulator shown n Fg., the angular dsplacement of each rgd component n the system s a known parameter due to a fxed connecton of the jont and the lnk. herefore, the system boundary condtons are,,,, m n, n x,, n,n, q, y x,, q n,n y,,,, n,n,,,, n, n (4) As shown n equaton (4), the angular dsplacements of rgd body and the end effector,.e.,, and n, n, are known parameters. he nput end dsplacements x and y of rgd body are always equal to and the nternal forces q and q, and the nternal moment m, are x, y, unknown parameters. For the end effector n+, the translatonal dsplacements xn, n and yn, n are unknown parameters, and the nternal force and nternal moment are equal to. he unknown parameters n vectors, and n, n can be calculated by substtutng equaton (4) nto equaton (). hen, wth the transfer equaton for each rgd component, the state vector of the connectng pont for each jont and lnk n the system can be obtaned. Furthermore, through lnearaton, the speed, the acceleraton, the angular velocty and the angular acceleraton of each pont can be determned. Detals of the modellng and soluton procedure are shown n Fg. 4. start Dsassemble system;ntale parametersθ,etc.; determne,,, n+,n alculate component x, y,x,y, determne each component U and =U Assemble and abtan U all, n,n+ =U all, Unknown parameter n n,n+ and, at tme t Unknown parameter n state vectors,,3 n,n- at tme t t m End m= t m < m=m+ Fg. 4. Soluton Procedure for the Unfed Model and Forward Knematcs of A Manpulator B.Soluton for nverse knematcs and model f the moton law for the end effector of the manpulator shown n Fg. s known, then the prmary focus of the soluton s to determne the moton law for the jonts. hs s because the moton and the force for every pont n the rgd body,.e., every element n the transfer matrx, are functons of the angular dsplacement parameter. Here boundary condton n equaton (4) ncludng the angular dsplacement of the nput end of jont,,, the nternal forces q x, and q y,, and the nternal moment m, are all unknown parameters. he rotaton angle n, n of the end effector s also an unknown parameter. he translatonal dsplacements xn, n and yn, n are known parameters, and the nternal force and nternal moment are both equal to. he soluton procedure s to substtute the system ntal boundary condton n equaton (4) and to calculate the ntal angular dsplacements of each jont at the moment of ntalaton. hen, accordng to the procedure shown n Fg.4, the known condton from equaton () s substtuted to determne all nformaton of the state vector for the robotc arm. n addton, the parameters of the

5 transfer matrx at ths moment and the boundary condton at the next moment are used to calculate the jont angular dsplacement and the state vectors of the unknown parameters at the next moment. he state vectors for all the jonts and lnks connectng ponts n the system at all moments can be obtaned by these teratons. Furthermore, through lnearaton, the speed, acceleraton, angular velocty and angular acceleraton for every pont are determned. he modellng and soluton procedures are shown n Fg.5. start Dsassemble system;ntale parameters; determne,,, Assemble and Uabtan all, =U all n,n+ n+,n alculate component x,y,x, y, determne each component U and =U, central rgd body (, 3) lnk (, 4) mass m (kg) =,3 radus r (m) =,3 mass centre (x,y )(m) =,3 nput (x,y ) (m) =,3 output (x,y ) (m) =,3 mass centre rotaton nerta J (kg.m ) =,3. (, ) (, ) (r, ).5 mass m j (kg) j=,4 length l j (m) j=,4 mass centre (x j,y j ) (m) j=,4 nput (m) (x j,y j )j=,4 output (x j,y j ) (m) j=,4 mass centre rotaton nerta j (kj g.m ) j=,4 (.5, ) (, ) (, ) /3 ABLE. he l parameters of a DF manpulator n ABLE, m s the mass of the rgd body. r s the radus of the central rgd body. l ndcates the length of the lnk. (x,y ), (x,y ) and (x,y ) are the mass centre, nput and output of rgd body respectvely. J s the mass centre rotaton nerta of the rgd body. Unknown parameter n and at tme n,n+, m=m+ m= Unknown parameter n state vectors,,3 n,n- at tme t t parallel to the nerta system x-axs, then the alculate rotaton angles of jonts andmoton law for jont and jont s n each state t vector ; at tme accordng to boundary condton + at tme t t m < t t m d t t sn t End d t Fg. 5. Soluton Procedure for nverse Knematcs and nverse Dynamcs (5) of A Manpulator V. Numercal Smulaton and Analyss o valdate the mproved mult-body system dscrete tme transfer matrx method, numercal smulatons are performed on a DF manpulator wth planar moton as shown n Fg.6. he structural parameters of the arm are shown n ABLE. Based on the moton stuatons of the arm, two stuatons are smulated. A. Soluton for forward knematcs model and numercal smulaton he soluton for the forward knematcs model s llustrated by the procedure shown n Fg. 4. At ths moment, the jont moton stuaton of the manpulator s known and the pose and dynamcs parameters of the end effector are unknown. f the arm ntal poston s two lnks n a straght lne, d where (=, ) s the expected value of the jont moton, d d 3, s the overall duraton of the moton, and are the moton ntalsaton and moton termnaton delay tmes, respectvely, = ms, =5s and '= ms..5 Y X Fg. 6. Dagram of A DF Manpulator wth Planar Moton m m Fg. 7. Robotc Arm Moton onfguraton Dagram

6 Dsplacement n the x-drecton x/m x x x3 x4 x_robot x_robot x3_robot x4_robot Fg. 8. Dsplacement of the Rgd omponents for the Manpulator n the x-drecton Notably, when calculatng the speed x, x, x, y, y, and y n equatons (3) and (4) at every pont, the mplct functon θ of tme t should be calculated durng the soluton because θ of every jont s a functon of tme t. Fg.7 shows the confguraton dagram of jont moton for the entre manpulator. Fg.8 and Fg.9 are dsplacements curves n the x- and y-drectons, respectvely. For manpulator rgd bodes wth -4 output ends, the data are calculated by the mproved transfer matrx method, and shown as x, x, x 3, x 4 and y, y, y 3, y 4 n the dagrams. o valdate ths method, the dagrams also show dsplacements curves for each component calculated by the D-H method, shown as x _robot, x _robot, x 3_robot, x 4_robot and y _robot, y _robot, y 3_robot, y 4_robot n the dagrams. alculaton results from the two methods match perfectly. Dsplacement n the y-drecton y/m y y y3 y4 y_robot y_robot y3_robot y4_robot Fg. 9. Dsplacement of the Rgd Bodes of the Manpulator n the y-drecton Fg. shows the angular dsplacement for the moton of each rgd body of manpulator. Jont and lnk have dentcal rotaton angular dsplacements, and jont 3 and lnk 4 have dentcal rotaton angular dsplacements. Fg. shows the nternal moments for all the rgd bodes. n ths dagram, m and m 3 are the drvng moments exactly to mplement moton n jonts and 3. Fg. and Fg.3 are nternal forces for the rgd bodes n the x- and y-drectons, respectvely. he data show that, the rgd body stress the maxmal forces q x and q y, whch are equal to the values exerted on the nput end of rgd body. Angular dsplacements / o Fg.. Angular Dsplacements of the Rgd Bodes of the Manpulator nternal moments m / N.m m m m 3 m Fg.. nternal Moments of the Rgd Bodes of the Manpulator nternal forces n the x-drecton qx/n q x q x q x3 q x Fg.. nternal Forces of the Rgd Bodes of the Manpulator n the x-drecton nternal forces n the x-drecton q y /N q y q y q y3 q y Fg. 3. nternal fforces of the Rgd Bodes of the Manpulator n the y-drecton

7 B. Soluton of nverse knematcs and Numercal smulaton he procedure for the nverse knematc numercal soluton s shown n Fg.5. At ths moment, the pose of the end effector s known, but the jont moton and ts dynamcs parameters are unknown. n ths case, the orgnal pont of the moton of the end effector s (.3, ), movement n the x-drecton s fxed,.e., x x start xend. 3, and moton n the y-drecton s defned as y.t.t.t.t t.( y start. ).8. y y y (6) where y start ; s, 4s, and 5s. Fg.4 shows the moton confguraton dagram of the manpulator. Fg.5 shows the dsplacement, speed and acceleraton of the end effector n the y-drecton. Fg.6 and Fg.7 show dsplacements of each rgd component n the x- and y-drectons, respectvely. he results agree wth the calculatons performed usng the robotcs method, whch valdates the rghtness of the mult-body system transfer matrx method. Fg.8 shows the angular dsplacements for every rgd body moton. he data show that the jont and lnk have consstent angular dsplacement moton; the jont 3 and lnk 4 have consstent angular dsplacement moton. Fg.9 shows the nternal moments for all the rgd bodes. n the Fg.9, m and m 3 are two jonts drvng moments. Fg. and Fg. show the nternal forces of all the rgd bodes n the x- and y-drectons, respectvely. he data show that rgd body are exerted on the maxmal forces q and q. m Dsplacement, speed and acceleraton of the robotc arm end n the y-drecton m Fg. 4. Manpulator Moton onfguraton Dagram y_end v_end a_end Fg. 5. Dsplacement, Speed and Acceleraton of the Manpulator End n the y-drecton x y Dsplacement n the x-drecton x/m x x x3 x4 x_robot x_robot x3_robot x4_robot Fg. 6. Dsplacement of the Rgd Bodes of the Manpulator n the x-drecton Dsplacement n the y-drecton y/m y y y3 y4 y_robot y_robot y3_robot y4_robot Fg. 7. Dsplacement of the Rgd Bodes of the Manpulator n the y-drecton Angular dsplacements / o Fg. 8. Angular Dsplacements of the Rgd Bodes n the Manpulator nternal moments m / N.m Fg. 9. nternal Moments of the Rgd Bodes n the Manpulator m m m 3 m 4

8 nternal forces n the x-drecton q x /N 3 - q x q x q x3 q x Fg.. nternal Forces on the Rgd Bodes of the Manpulator n the x-drecton nternal forces n the x-drecton q y /N.4. q y.6 q y. q y3.8.4 q y Fg.. nternal Forces on the Rgd Bodes of the Manpulator n the y-drecton V. onclusons n hs study, a seral manpulator wth multple degrees of freedom s consdered as the subject. Based on an mproved dscrete tme transfer matrx method for mult-body systems, a new unfed model s developed, ncorporatng both knematcs and dynamcs. A soluton method s also provded for ths model. he soluton s tested on a DF arm through a numercal smulaton. he followng conclusons are reached: () Accordng to the structural characterstcs of a seral robot and the relaton between moton transfer and force transfer, the mult-body system transfer matrx method s mproved to obtan a mult-rgd-body dscrete tme transfer matrx equaton applcable to a robotc arm. By usng ths equaton, an n-jont seral planar mult-body system dynamcs model for a manpulator s derved. hs model covers both the knematcs and dynamc characterstcs of the system; () Usng ths unfed model, and applyng the two dfferent moton stuatons,.e., jont moton and end effector moton, the correspondng boundary condtons are set. A soluton method and procedure for the forward knematcs problem and the nverse knematcs problem are proposed. n addton, a soluton method and procedure for the nverse dynamcs of the model are proposed; (3) he knematcs and dynamcs numercal smulaton of a DF planar manpulator s completed. he results show that dsplacements n the x- and y-drectons and the angular dsplacements of the rgd bodes n the arm agree wth calculated results performed by a classc robotcs method. hs proves the valdty and feasblty of the model. he model soluton also yelds characterstcs of the dynamcs, whch llustrates the superorty of the proposed model to prevous methods. n summary, the mproved dscrete tme transfer matrx method for mult-body systems proposed n ths paper provdes a new method and a new approach for robot modellng. A mult-body system model of a multple degrees-of-freedom robotc arm based on ths method unfes the study of system knematcs and dynamcs, and provdes a theoretcal foundaton for the further study of robot knematcs and dynamcs characterstcs. he results from the numerc smulaton of the DF planar manpulator llustrated the valdty and feasblty of the proposed modellng method and the superorty of the model over exstng methods. Acknowledgment he authors would lke to acknowledge the fnancal support of the natonal scence foundaton of chna (347, ), hna Postdoctoral Scence Foundaton (Mltary System: M57), and anjn Key Laboratory of the Desgn and ntellgent ontrol of the Advanced Mechatroncal System. References [] Wang X.B. Research on Structure System Dynamcs and haracterstcs of Engneerng Mechancal Arm. Zhejang(PR hna): Zhejang Unversty, 4. [] Xong Y.L., Dng H. And Lu E.. Robotcs. Bejng(PR hna): Machnery ndustry Press, 993. [3] Yang W.L, Dong W. and Du Z.J. Knematcs Modelng for a Knematc-mechancs ouplng ontnuum Manpulator. n nternatonal onference on Manpulaton, Manufacturng and Measurement on the Nanoscale (3M-NAN), ape(w), pp , 4. [4] Al-Mashhadany Y.. nverse Knematcs Problem (KP) of 6-DF Manpulator by Locally Recurrent Neural Networks (LRNNs). n nternatonal onference on Management and Servce Scence (MASS), Wuhan(PR hna), pp.-5,. [5] Korayem M.H., Abbas Esfeden R. and Nekoo S.R. Path Plannng Algorthm n Wheeled Moble Manpulators Based on Moton of Arms. Journal of Mechancal Scence and echnology, 9(4): , 5. [6] Xe L.M. and hen L. Robust and adaptve composte control of space manpulator system wth bounded torque nputs. Engneerng Mechancs, 3(3): , 3 (n hnese). [7] Allan J.F., Lavoe S., Reher S. and Lambert G. Knematc and Dynamc Analyss of A Novel 6-DF Seral Manpulator for Underground Dstrbuton Power Lnes. n EEE/RSJ nternatonal onference on ntellgent Robots and Systems. San Francsco, A(USA), pp ,. [8] Yaskevch N. Recursve Algorthm for Smulaton of Sngle han Manpulator Dynamcs. n EEE/ASME nternatonal onference on Mechatronc and Embedded Systems and Applcatons,Sengalla, Ancona(taly), pp.-5, 4. [9] Ru X.., Yun L.F. and Lu Y.Q., etal. ransfer matrx method of mult-body system and ts applcaton. Scence Press (PR hna): 8. [] Hong J.Z. omputatonal dynamcs of mult-body systems. Hgher Educaton Press (PR hna): 999.

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

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

Modeling of Dynamic Systems

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

More information

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

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

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

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

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

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

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

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

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

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

Physics 111: Mechanics Lecture 11

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

More information

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

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

More information

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

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

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

More information

Lecture 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

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

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

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

So far: simple (planar) geometries

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

More information

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,

More information

DYNAMIC MODELING AND PATH PLANNING OF FLEXIBLE-LINK MANIPULATORS

DYNAMIC MODELING AND PATH PLANNING OF FLEXIBLE-LINK MANIPULATORS Int. J. Mech. Eng. & Rob. Res. M Nazemzadeh et al., Research Paper ISSN 78 9 www.merr.com ol., No., Aprl IJMERR. All Rghts Reserved DYNAMIC MODELING AND PAH PLANNING OF FLEXIBLE-LINK MANIPULAORS M Nazemzadeh

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

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

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

χ 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

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites 7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT

More information

Notes on Analytical Dynamics

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

Rigid body simulation

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

More information

MEEM 3700 Mechanical Vibrations

MEEM 3700 Mechanical Vibrations MEEM 700 Mechancal Vbratons Mohan D. Rao Chuck Van Karsen Mechancal Engneerng-Engneerng Mechancs Mchgan echnologcal Unversty Copyrght 00 Lecture & MEEM 700 Multple Degree of Freedom Systems (ext: S.S.

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

A Fast Computer Aided Design Method for Filters

A Fast Computer Aided Design Method for Filters 2017 Asa-Pacfc Engneerng and Technology Conference (APETC 2017) ISBN: 978-1-60595-443-1 A Fast Computer Aded Desgn Method for Flters Gang L ABSTRACT *Ths paper presents a fast computer aded desgn method

More information

Chapter 11: Angular Momentum

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

More information

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

1 axis. F Input force for i th joint. f Force exerted on linki by linki 1at the coordinate frame ( i 1

1 axis. F Input force for i th joint. f Force exerted on linki by linki 1at the coordinate frame ( i 1 Internatonal Journal of Robotcs and Automaton (IJRA) Vol. 3, No., December, pp. ~33 ISSN: 89-856 Dynamc Behavor of a SCARA Robot by usng N-E Method for a Straght Lne and Smulaton of Moton by usng Soldworks

More information

Parameter Estimation for Dynamic System using Unscented Kalman filter

Parameter Estimation for Dynamic System using Unscented Kalman filter Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

Tensor Smooth Length for SPH Modelling of High Speed Impact

Tensor 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

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

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

Conservation of Angular Momentum = "Spin"

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

More information

Chapter 11 Angular Momentum

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

More information

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

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

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

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

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

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

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

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

Markov Chain Monte Carlo Lecture 6

Markov Chain Monte Carlo Lecture 6 where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways

More information

Classical Mechanics ( Particles and Biparticles )

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

More information

Physics 181. Particle Systems

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

More information

High resolution entropy stable scheme for shallow water equations

High resolution entropy stable scheme for shallow water equations Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal

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

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

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

AP Physics 1 & 2 Summer Assignment

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

More information

Report on Image warping

Report 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

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:

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

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST Appled Computer Scence, vol. 13, no. 4, pp. 56 64 do: 10.23743/acs-2017-29 Submtted: 2017-10-30 Revsed: 2017-11-15 Accepted: 2017-12-06 Abaqus Fnte Elements, Plane Stress, Orthotropc Materal Bartosz KAWECKI

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

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

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros

More information

Lesson 5: Kinematics and Dynamics of Particles

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

More information

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model Identfcaton of Instantaneous Modal Parameters of A Nonlnear Structure Va Ampltude-Dependent ARX Model We Chh Su(NCHC), Chung Shann Huang(NCU), Chng Yu Lu(NCU) Outlne INRODUCION MEHODOLOGY NUMERICAL VERIFICAION

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

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

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 539 Homeworks Sprng 08 Updated: Tuesday, Aprl 7, 08 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. For full credt, show all work. Some problems requre hand calculatons.

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

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE October 1-17, 008, Bejng, Chna NON LINEAR ANALYSIS OF SRUCURES ACCORDING O NEW EUROPEAN DESIGN CODE D. Mestrovc 1, D. Czmar and M. Pende 3 1 Professor, Dept. of Structural Engneerng, Faculty of Cvl Engneerng,

More information

Adaptive Consensus Control of Multi-Agent Systems with Large Uncertainty and Time Delays *

Adaptive Consensus Control of Multi-Agent Systems with Large Uncertainty and Time Delays * Journal of Robotcs, etworkng and Artfcal Lfe, Vol., o. (September 04), 5-9 Adaptve Consensus Control of Mult-Agent Systems wth Large Uncertanty and me Delays * L Lu School of Mechancal Engneerng Unversty

More information

12. The Hamilton-Jacobi Equation Michael Fowler

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

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

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

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

More information

Operating conditions of a mine fan under conditions of variable resistance

Operating conditions of a mine fan under conditions of variable resistance Paper No. 11 ISMS 216 Operatng condtons of a mne fan under condtons of varable resstance Zhang Ynghua a, Chen L a, b, Huang Zhan a, *, Gao Yukun a a State Key Laboratory of Hgh-Effcent Mnng and Safety

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

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

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016 ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have

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

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

VIBRATION FATIGUE DESIGN METHODOLOGY OF A LARGE SCALE HEAVY DUTY ROBOT

VIBRATION FATIGUE DESIGN METHODOLOGY OF A LARGE SCALE HEAVY DUTY ROBOT ICSV14 Carns Australa 9-12 July, 2007 VIBRATION FATIGUE DESIGN METHODOLOGY OF A LARGE SCALE HEAVY DUTY ROBOT Jong Hw Seo 1, Jae Chul Hwang 1, Yong Won Cho 1, Dong Il Km 1, Hong Jae Ym 2 1 Robotcs Technology

More information

Lecture 12: Discrete Laplacian

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

More information

Modeling and Simulation of a Hexapod Machine Tool for the Dynamic Stability Analysis of Milling Processes. C. Henninger, P.

Modeling and Simulation of a Hexapod Machine Tool for the Dynamic Stability Analysis of Milling Processes. C. Henninger, P. Smpack User Meetng 27 Modelng and Smulaton of a Heapod Machne Tool for the Dynamc Stablty Analyss of Mllng Processes C. Hennnger, P. Eberhard Insttute of Engneerng project funded by the DFG wthn the framework

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

An Improved multiple fractal algorithm

An Improved multiple fractal algorithm Advanced Scence and Technology Letters Vol.31 (MulGraB 213), pp.184-188 http://dx.do.org/1.1427/astl.213.31.41 An Improved multple fractal algorthm Yun Ln, Xaochu Xu, Jnfeng Pang College of Informaton

More information

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

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

More information

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

Chapter - 2. Distribution System Power Flow Analysis

Chapter - 2. Distribution System Power Flow Analysis Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load

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

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order

More information

The Finite Element Method

The Finite Element Method The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of

More information

Implicit Integration Henyey Method

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

More information

Physics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1

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

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

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

More information