DIFFERENTIAL KINEMATICS OF CONTEMPORARY INDUSTRIAL ROBOTS

Size: px
Start display at page:

Download "DIFFERENTIAL KINEMATICS OF CONTEMPORARY INDUSTRIAL ROBOTS"

Transcription

1 Int. J. of Appled Mechancs and Engneerng 0 vol.9 No. pp.-9 DOI: 0.78/jame-0-00 DIFFERENTIAL KINEMATICS OF CONTEMPORARY INDUSTRIAL ROBOTS T. SZKODNY Insttute of Automatc Control Slesan Unversty of Technology ul. Akademcka -00 Glwce POLAND E-mal: tadeusz.szkodny@polsl.pl The paper presents a smple method of avodng sngular confguratons of contemporary ndustral robot manpulators of such renowned companes as ABB Fanuc Mtsubsh Adept Kawasak COMAU and KUKA. To determne the sngular confguratons of these manpulators a global form of descrpton of the end-effector knematcs was prepared relatve to the other lnks. On the bass of ths descrpton the formula for the Jacoban was defned n the end-effector coordnates. Next a closed form of the determnant of the Jacoban was derved. From the formula sngular confguratons where the determnant s value equals zero were determned. Addtonally geometrc nterpretatons of these confguratons were gven and they were llustrated. For the exemplary manpulator small correctons of jont varables preventng the reducton of the Jacoban order were suggested. An analyss of postonal errors caused by these correctons was presented. Key words: knematcs manpulators mechancal system robot knematcs.. Introducton The controllers of such renowned companes as ABB Fanuc Mtsubsh Adept Kawasak COMAU and KUKA make possble movement programmng n jont space or the Cartesan space. The followng commands PTP LIN and CIRC can be appled to programmng n the Cartesan space. The mentoned commands requre a start pont and end pont. For programmng n the Cartesan space these ponts must be descrbed n Cartesan coordnates relatve to the base frame connected to the base of a manpulator. These coordnates can be obtaned by a vson system. Durng the realzaton of such programmed movement the robot happens to stop before reachng the border area and before reachng the start or the end pont. The entrapment takes place when the manpulator reaches the sngular confguratons. It s wthout doubt the major problem of modern ndustral robots whch makes t mpossble for the robots wth a vson system to cooperate wth the cameras properly. The lnear system of ordnary dfferental equatons that descrbes the knematcs can be appled to programmng the robots n the Cartesan space. In ths system a manpulator Jacoban s present. For the succesve va ponts nterpolatng trajectory programmed n the Cartesan space jont varables of these ponts can be computed. To these calculatons the standard algorthm for solvng the system of lnear equatons whch s one of the elements of a computer software lbrary can be appled. The lack of protecton n software aganst the reducton of the Jacoban rank n ths algorthm can cause the nterrupton of the calculatons and the robot can stop performng ts operatons. Ths rank decreases n sngular confguraton. The problem of the nverse knematcs solutons n a dfferental form for the sngular confguraton s presented n Chaccho (99) Kozłowsk (00) Nakamura (009) Sclano (00) Spong (997) Tchoń (997). In Chaccho (99) Kozłowsk (00) Nakamura (009) Sclano (00) for a sngular confguraton t s proposed to use: sngular value decomposton technques SVD of the Jacoban the damped least-squares nverse of the Jacoban DLS or sngularty robust nverse of the Jacoban. In the work

2 T.Szkodny (Tchoń 997) a method of avodng sngulartes by usng of dynamcal systems of the own moton s presented. The same work presents a soluton of the nverse knematcs n sngular confguratons obtaned usng: the method a normal form Jacoban attached the zero and the Jacoban space. A common feature of knematc sngularty avodance methods presented n Chaccho (99) Kozłowsk (00) Nakamura (009) Sclano (00) Tchoń (000) s ther hgh computatonal complexty. A very smple approach to the problem of determnng the knematc sngulartes based on the dfferental descrpton s presented n Spong (997) on the examples of manpulators wth three degrees of freedom. Smple because based on closed forms of the determnants of manpulators Jacoban allowng a very smple determnaton of jont varables descrbng sngular confguratons. Smple rules of lnear algebra are recommended n Spong (997) for solvng the nverse knematcs n the feld of speed. Only the methods whch do not requre a large number of calculatons can be of practcal use. Ths paper presents a very smple method of correcton allowng to avodance of the sngular confguraton of modern ndustral robots. In ths method the values of thrd and ffth jont varables are corrected. These values were determned from the closed form of the Jacoban determnant of modern ndustral robots (Spong 997). The closed form of formulas descrbng the thrd and ffth jont varables for the sngular confguratons of manpulators was determned. On the bass of these formulas one s able to determne whether n the next step of the numercal calculatons the freezng of the program caused by a reducton of the Jacoban rank wll take place. These correctons prevent the reducton of the Jacoban rank and allow the use of any method of solvng the nverse knematcs ncludng the methods descrbed n the aforementoned works (Chaccho 99; Nakamura 009; Sclano 00; Tchoń 000). In the second chapter the knematc structure and a descrpton of the end-effector knematcs n relaton to the other lnks ncludng the base of the manpulator are presented. The formulas consttutng the dfferental descrpton of the knematcs are presented n the thrd chapter. A general descrpton of sngular confguratons ncludng ther llustraton and the examples of the calculatons of jont varables correctons preventng the reducton of the Jacoban rank are presented n the fourth chapter. The ffth chapter summarzes the paper.. The global descrpton of the knematcs Fgure llustrates the robot manpulators of the majorty of modern robots wth numbered lnks. The lnk 0 s attached to the ground other lnks - are movable. The lnk 0 wll be called a base lnk the last lnk wth number wll be called an end-effector. The grpper or another tool s attached to ths lnk. Neghborng lnks are connected by revolute jonts. Fgure llustrates the manpulator knematcs schema wth the co-ordnate systems (frames) assocated wth lnks accordng to a Denavt-Hartenberg notaton (Jezersk 00; Kozłowsk 00; Szkodny 0a; 0b). The x 7 y 7 z 7 frame s assocated wth the grpper. The poston and orentaton of the lnks and tool are descrbed by homogenous transform matrces. Matrx A descrbes the poston and orentaton of the -th lnk frame n relaton to --st. T s a matrx that descrbes the poston and the orentaton of the end-effector frame n relaton to the base lnk. Matrx E descrbes the grpper frame n relaton to the end-effector frame. Matrx X descrbes the poston and the orentaton of the grpper frame n relaton to the base lnk. Parameters l are used n the Denavt-Hartenberg notaton. For further descrpton of the knematcs jont varables wll be used. Varables for = and 90. To smplfy the followng denotatons wll be used: S sn C cos j S sn C cos j j j j j. Matrces A are descrbed by Eq.(.a). C 0 S lc S 0 C ls A S C 0 ls C S 0 lc A

3 Dfferental knematcs of contemporary ndustral robots C 0 S lc C 0 S 0 S 0 C ls S 0 C 0 A A lnk lnk lnk lnk lnk lnk 0 lnk Fg.. The manpulator. O x x λ x x x lnk lnk y y z y z x A l O=O A z O A λ x y lnk z λ7 E x7 z7 O7 y7 l A y lnk lnk O y z A y O lnk x x0 T X grpper λ y0 z0 O0 l x0 z A lnk 0 Fg.. Knematc scheme of the manpulator Denavt-Hartenberg parameters jont varables and jont frames.

4 T.Szkodny C 0 S 0 S 0 C 0 A C S 0 0 S C 0 0 A. (.a) Matrx E s descrbed by Eq.(.b) E. (.b) Matrx - T descrbes the end-effector frame knematcs n relaton to --st frame. Equaton (.a) (Szkodny 0a) descrbes these matrces. T A. (.a) j From Eqs (.a) and (.a) one obtans j C S 0 0 S C 0 0 T CC CS S S SC SS C C T S C CCC SS CCS SC CS CS SCC CS SCS CC SS SS T SC SS C C C C C C S S S S C C C C S S C S S S S C C C S S C S C S C C S S C C S S T SCC CS SCS CC 0 0 CCS SC lc CS CCS SCS CC ls CC SCS SS SS 0

5 Dfferental knematcs of contemporary ndustral robots 7 S C C C S S C S C S C C S S C C S S C C C C S S S S C C C C S S C S S S T SCC CS SCS CC 0 0 SCS CC ls ls C C SCS CCS SC lc lc CS CCS SS SS 0 C C C C S S S S C C S CCC SS CSC S SCC CS S S CCC SS CSC C SCC CS T 0 CSCCS SC CSS SSCS CC SSCCS SC CSS CSCS CC C CCS SC SSS 0 CSCS CCSSS SSCS CCCSS CCS SC 0 lc SS S ls ls CC SCSC ls CSS ls ls CC SCS S. (.b) lc CS lc CS These matrces are necessary to the dfferental descrpton of the end-effector knematcs n relaton to the base frame. The descrpton s presented n the next chapter.. The dfferental descrpton of the knematcs To make the descrpton of the manpulator knematcs ndependent of the shape of the grpper the requred poston and orentaton of the grpper frame x 7 y 7 z 7 (descrbed by the matrx X ) are converted to the end-effector frame x y z. Correlaton Treq XreqE makes the converson possble. Therefore n further consderatons we wll focus on the end-effector knematcs. The movement of the x 7 y 7 z 7 frame n relaton to the x 0 y 0 z 0 frame wll be descrbed by usng dsplacement dfferentals 7 dx 7 7 dy 7 7 dz 7 of orgn O 7 and x-y-z current angle dfferentals 7 d 7x 7 d 7y 7 7z d (Crag 99; Szkodny 0). These dfferentals are descrbed n the x 7 y 7 z 7 frame. req

6 8 T.Szkodny Smlarly the movement of the x y z frame n relaton to the x 0 y 0 z 0 frame wll be descrbed by usng the dsplacement dfferentals dx dy dz of the orgn O and x-y-z current angle dfferentals d x d y d z. These dfferentals are descrbed n the x y z frame. Further we wll focus on the movement descrpton of the end-effector and therefore one t s necessary to convert the grpper 7 dfferentals to end-effector dfferentals. The dfferental equaton T E X whch connects 7 end-effector and grpper dfferentals results from the work (Szkodny 0a). In ths equaton and 7 7 are dfferental transformaton matrces respectvely of the end-effector and grpper frame. These matrces have the followng forms 0 dz dy dx dz 0 dx dy dy dx 0 dz d7z d7y dx d7z 0 d7x dy d7y d7x 0 dz7 From the above Eq.(.) results 7 7 T X 7E E 7E. (.) From Eqs (.) and (.b) results Eq.(.) whch makes t possble to compute the end-effector dfferentals from grpper dfferentals. 7 d x d 7x 7 7 dx dx7 7 d 7y 7 d y d 7y 7 7 dy dy7 7 d 7x 7 d z d 7 dz dx7. 7z (.) To descrbe the end-effector knematcs one wll apply the Cartesan dfferental matrx D and the jont dfferental matrx dx dy dz dx dy dz T T dq d d d d d d (Szkodny 0). Equaton (.) connectng these matrces s a dfferental descrpton of the end-effector knematcs. d D J q. (.) In ths equaton an end-effector Jacoban J s present descrbed n the x y z frame whch has the form of Eq.(.).

7 Dfferental knematcs of contemporary ndustral robots 9 x x x x x x ' ' ' ' ' ' y y y y y y ' ' ' ' ' ' z z z z z z ' ' ' ' ' ' J x x x x x x ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' y y y y y y z z z z z z. (.) Fgure llustrates the dfferental dsplacement and the rotaton of the end-effector caused by the dfferental dsplacement and the rotaton of the x y z frame. The vectors dr and d respectvely descrbe the dsplacement and rotaton of the x y z frame n relaton to the x - y - z - frame caused by a dfferental ncrease of the jont varable d. The x y z frame s connected wth the -th lnk and concdes wth x - y - z - frame for d 0. The dsplacement of dr of x y z frame s caused by the dsplacement of d r and rotaton of vector d by an angle d. Therefore T z- z d dr z d z d y y x- d x dr x y y- dr x d d Fg.. Dfferental dsplacements and rotatons of -th and -th frames. d d d r r d. (.) Ths movement wll be recalculated n relaton to the x y z frame usng Eq.(.). dx d a r dy b dr dz c dr. (.) The vectors a b c are the versors of the x y z frame descrbed n the x - y - z - frame. The rotaton of the x y z frame s the same as the rotaton of x y z frame n relaton to the x - y - z - frame

8 0 T.Szkodny whch s a vector d. Ths vector can be recalculated n relaton to the x y z frame by the followng Eqs (.7) x d d a dy b d dz c d. (.7) In Fg. one can see that all lnks are connected by revolute jonts. Therefore the angles d are vectors drected along the z - axs havng a length equal to the dfferental d. Therefore d kd. (.8) k s the versor of the z - axs descrbed n the x - y - z - frame. For such jonts no dsplacement of the coordnate system x y z n relaton to the x - y - z - frame takes place therefore dr 0. Thus from Eqs (.) and (.8) one obtans the followng dr d d dy j dx d. (.9) The vectors and j are versors of the x - and y - axes respectvely descrbed n x - y - z - frame. dx and dy are the x and y - coordnates of the vector d n the x - y - z - frame. After takng nto account Eqs (.9) and (.) one obtans the followng dx ax dy ay dx d dy bx dy by dx d (.0) dz cx dy cy dx d. a x ay are the x and y - coordnates of the versor a n x - y - z - frame. The coordnates of the versors b and c were marked n a smlar way. From Eqs (.7) and (.8) one obtans the followng correlatons d x az d d y bz d d z cz d. (.) The dfferental of each coordnate n the Cartesan matrx D s the sum of the correspondng coordnate dfferentals caused by dfferental changes n jont varables d. For example dx dx. In the Jacoban matrx J the dervatves of Cartesan coordnates of the x y z endeffector frame wth respect to the jont varables are present. The dervatve x x only because the

9 Dfferental knematcs of contemporary ndustral robots dfferental dx depends on the value of d. The stuaton s smlar wth other dervatves n ths matrx. Therefore the elements of matrx J can be presented n the form of Eqs (.) and (.) takng nto account the relatons (.0) and (.). x x ax dy ay dx y y z z b d b d x y y x x y y x c d c d (.) x x y y a z b z z z c z. (.) Matrces - T descrbed by Eq.(.b) can be presented n the form Eq.(.). T. (.) ax bx cx dx ay by cy dy az bz cz dz Quanttes appearng n formulas (.) and (.) can be replaced by the correspondng elements of the matrx - T resultng from the form Eq.(.b) and correspondng to the form Eq.(.). It s easy to notce that the ndex n Eqs (.)-(.) s the number of column n the matrx J. After usng Eqs (.b) (.)-(.) and smplfcatons one obtans the followng elements of the end-effector Jacoban J x l ( S C C C S )+l S (S C C +C S )+l S (S C C +C S )+ -λ C (S C C +C S )+λ ( C C S C S C C +S S S ) y l (S C S C C ) l S (S C S C C ) l S (S C S C C )+ +λ C (S C S C C )+λ (C S S +S C S C C C C )

10 T.Szkodny z l S S +l S S S +l S S S λ S S C x C S S S C S +C C C C y C S C C C S C +S S S z C C S +S C ; x l ( C C C C +C S S +C S S )S +(S C C C S S S +C C S )C )+ +l C S +λ (C C C S S ) λ (S S C C C ) y l (C C C S +C S C S S S )S (S C C S +S S C +S C S )C )+ -l S S +λ ( C C S S C ) λ (S C C +C S ) z l ((S C S C C )C (C C S +S C )S ) l C +λ C S x S C C +C S y S S C +C C z S S ; x l C S +λ (C C C S S )+λ ( S S C +C C ) y l S S +λ ( C C S S C )+λ ( C S S C C ) z l C +λ C S x S C C +C S y S S C +C C z S S ; x λ S S y λ S C z 0 x S C y S S z C ;

11 Dfferental knematcs of contemporary ndustral robots x λ C y λ S z 0 x S y C z 0; x 0 y 0 z 0 x 0 y 0 z. (.) Equaton (.) wth dependences (.) s a dfferental descrpton of the end-effector knematcs.. The sngular confguratons For the sngular confguratons of manpulators the J determnant equals zero. Thus to determne such confguratons one needs to calculate the determnant. The closed form of the determnant was obtaned by Symbolc Math Toolbox lbrary of Matlab. It s descrbed by Eq.(.). det J = l S (λ C l S )( l S l S +λ C +l ). (.) O O O Fg.. Confguraton at whch 0. The z and z axes are collnear. O * O O O O =l O O =λ O Fg.. Confguraton at whch C ls 0. Ponts O O and O are stuated on a straght lne.

12 T.Szkodny The sngular confguratons appear n zero values of the factors on the rght sde of Eq.(.). The jont varables at whch t wll take place are descrbed by Eq.(.). * * C l S 0 arc tg l S 0 0 ls ls C l 0 ** l l ls l () arctg ls l (.) O O ** () z 0 O l O O O =l O O =l O O =λ ** ** ( ) ( ) ** ** ( ) ( ) O z 0 ** ( ) O O O Fg.. Confguraton at whch ls ls C l 0. Pont O s placed on the lne passng through the z 0 axs. Fgures and llustrate these confguratons. For example let us analyze knematc sngulartes of the manpulator IRB-00. The manpulator has the followng parameters (Szkodny 009) l 0mm l 00mm l 0 mm 7 mm 70 mm 8 mm. The ranges of angle changes and are the followng: 70. The calculatons below were mode n Matlab on a PC from an Intel Pentum processor wth a frequency of GHz. Let us assume that for the confguraton from Fg.. For ths confguraton * 0 0. In ths case the det J =0 and the rank J =. The rank of the Jacoban J s less than and therefore a standard numercal soluton algorthm of Eq.(.) wll nterrupt the calculatons. One can avod zerong the determnant by ncreasng or decreasng the angle by a mnmum ncrement mn resultng from the resoluton of the encoder and the gear rato. Let us assume that a typcal resoluton of

13 Dfferental knematcs of contemporary ndustral robots encoders s the followng / rad. and the value of the gear rato equals 00. For such data mn. 0 0 rad. After the correctons mn one obtans det J. 70 and rank J =. Let us estmate the grpper poston error caused by ths correcton. The knematc scheme n Fg. results n the error not greater than. Let us assume that a typcal value of the mn 7 parameter grpper s the followng 7 0 mm. For such a grpper mn causes a change of the grpper poston not greater than. 0 mm. * For the confguraton from Fg. the jont varable The manpulator IRB-00 cannot reach ths varable because t s outsde the range of ts changes. Let us assume that for a confguraton from Fg.. ** For these jont varables one obtans (). 7 ** ( ) **. The jont varable ( ) ** s outsde the range of ts changes. For the jont varable () det J and rank J =. The rank of the Jacoban J s less than and therefore a standard numercal soluton algorthm of Eq.(.) wll nterrupt the calculatons. In order to prevent a decrease n the rank of the Jacoban J one wll ncrease and reduce the jont varable ** () by a mnmum ncrement mn equal to mn. After correctons ** mn of the jont varable () one obtans det J. 8 0 and rank J =. The change n mn causes an error of the grpper poston whch equals max mm.. Summary A dfferental descrpton of the manpulator knematcs presented n ths paper s the bass for desgn of contemporary software drvers of ndustral robots ndependent of the software manufacturers robots. Ths descrpton can be appled n IRB manpulators seres ; n Fanuc manpulators M M M70 M0 M900; n a KUKA manpulators KR KR KR KR; n a Mtsubsh manpulators RV-A RV-A RV-S RV-S RVS; n Adept manpulators s00 s0 s80 s700; n Kawasak manpulators seres M FS00N ZHE00U KF. For KUKA KRC Adept s00 and Mtsubsh RV-AJ manpulators the orgn of second and thrd lnk frame s at the same pont O =O (see Fg.). Therefore for these manpulators one has to assume that l 0. In modern ndustral robot controllers actuator varables calculaton s realzed n two stages n master-level of control and n a sngle stage n the slave-level control. In the frst step n the master-level the trajectory generator calculates the Cartesan coordnates of the trajectory va-ponts. The trajectory s specfed by a sutable Cartesan programmng command (for example LIN). The calculated va-ponts n the Cartesan space are pre-assgned tmes. The tmes resultng from the speed n the Cartesan space declared n a sutable command parameter are assgned to the ponts. In the second stage jont varables are determned (Crag 99; Szkodny 0) for the earler computed va-ponts n the Cartesan space. To determne these varables the teratve calculaton algorthm of solutons Eq.(.) s used. Next actuator varables are calculated for jont varables of these va-ponts. Such calculated actuator varables are sent from the master level nto the slave level of control wth constant frequency. The actuator varables of the va-ponts are calculated at the tme of the robot moton to the prevously calculated va-ponts. If a va pont determned by the trajectory generator n the Cartesan space s sngular calculatng jont varables usng the teratve algorthm wthout a protecton aganst the Jacoban rank reducton results n an undesrable

14 T.Szkodny nterrupton of the calculatons and stoppng the robot. The robot stops before the a sngular va-pont when approachng the prevously calculated nonsngular va-ponts. It may happen that the va-pont obtaned n the frst stage of the calculaton s near a sngular vapont. In ths case the teratve calculaton algorthm of jont varables may be convergent but the resultng actuator varables wll be very dfferent from the correspondng varables of the prevous va-pont. For large ncrements of jont varables and already pre-declared tmes of the prevous and current va-ponts t may happen that the average rate of change of actuator varables s much greater than the maxmum actuator speed. Furthermore n ths case the actuator varables acceleratons may occur above the maxmum acceleraton that can be reached by the actuators. The problem of these exceedances s solved by the trajectory generator n the slave-layer of control. Command moton parameters (e.g. LIN) are sent from the master level nto the slave level control. These are among others the declared translatonal movement speed and acceleraton. In general the acceleraton s declared n percent of the maxmum permssble values. These parameters are the bass for the generaton of actuator varables n the slave level of control. Waveforms of actuator varables are generated at a constant speed sectons connected wth fxed segments acceleratons. If the movement command e.g. LIN declares translatonal movement speed and acceleraton at the level of 0% of ther maxmum values the result s that the n slave level of control the waveforms wth speed and acceleraton rates of not more than 0% of the maxmum actuator speed and acceleraton wll be generated. In the controllers of contemporary ndustral robots one can wrte hs/her own applcatons. These applcatons can be wrtten usng Eq.(.) derved from ths work. In ths equaton dfferentals d should be replaced by the ncrements and dfferentals dx dy dz d x d y and d z by correspondng ncrements x y z x y and z. Such dscretzed Eq.(.) s the bass for the second stage of the calculaton at master level of control. Durng ths stage the jont varables of the va-pnts are teratvely calculated from the dscretzed Eq.(.). Recall that at ths level of control n the frst stage the Cartesan coordnates of the va-ponts are calculated. Assume that the jont varables of the prevous va-pont of trajectory were calculated teratvely and ( ) denote them by. For the calculaton of the jont varables of the current va-pont we wll also apply the teratve calculaton algorthm. In the next for example the k-th teraton of calculatons the ntal values of the lnk varables ( ) ( ) are equal to the end values ( ) ( ) from the prevous teraton. We start k ( ) start ( k ) calculate the Jacoban J and ncreases end k ( ) k ( ) end k from the dscretzed Eq.(.). The end values of start k the jont varables n the k-th teraton ( ) ( ) are the sum ( ) ( ) +. Subsequent teratve ( ) calculatons cause that the end values end( k) of teraton approach values representng the soluton of the nverse knematcs of the current va-pont. If the current va-pont s sngular the approachng t gves rse to the rsk of a reducton of the Jacoban ( ) start ( k ) ( ) k ( ) J rank. These protecton measures are presented n the fourth chapter by usng the example of the IRB-00 manpulator. They consst n the correcton of the angles and by mnmum ncrements mn and mn seen by the servos. ( ) ( ) * If durng the teratve calculatons satsfes the nequalty we assume: ( ) end ( k ) = ** () ( ) end ( k ) = * sgn ( ) k start( k) sgn ( ) k ( ) mn. If ( ) mn ( ) ** start( k) ( ) start k start( k) mn mn we substtute. Also for the angle ( ) ( ) satsfyng the nequalty

15 Dfferental knematcs of contemporary ndustral robots 7 ( ) * ( ) start( k) start( k) mn we assume: ( ) end = sgn ( ) ( ) k ( ) ( k ) mn. For the other varables end ( k ) = ( ) start( k). Usng such protecton measures one should check the range of the changes n the grpper poston caused by above angles correctons. For the sample IRB-00 manpulator suggested correctons can at most result n changes of the poston whch s the sum of errors.e.. 0 mm mm mm. The typcal postonng accuracy of ndustral robots are of the order 0 mm. Thus the errors caused by such correctons are acceptable. If the grpper poston errors caused by such correctons were equal to or greater than the requred postonng accuracy nverse knematcs equatons n a global form presented n Szkodny (00) would need to be used. For these equatons we can apply smple crtera for selectng solutons at the jont varables. Such crtera are much smpler and faster to mplement than the crtera proposed n the Jacoban methods (Chaccho 99; Nakamura 009; Sclano 00; Tchoń 000). The dfferental descrpton of the knematcs ncludng correctons to prevent the loss of the Jacoban rank presented n ths paper consttutes the bass for the creaton of our own software free from the basc defect namely undesred stoppng n sngular postons. Nomenclature A homogeneous matrx descrbng the -th lnk frame n relaton to the --st lnk a versor of x axs descrbed n --st frame ax ay az coordnates of versor a n --st frame b versor of y axs descrbed n --st frame bx by bz coordnates of versor b n --st frame c versor of z axs descrbed n --st frame c c c coordnates of versor c n --st frame x y z D Cartesan dfferental matrx d vector descrbes the poston of the -th lnk frame orgn n relaton to --st frame descrbed n --st frame dx dy dz the coordnates of the vectors d n --st frame d d d r dfferental dsplacement of the x y z frame orgn n relaton to the x - y - z - frame caused by a dfferental ncrease of the jont varable r dfferental dsplacement of the -th lnk frame orgn n relaton to the x - y - z - frame caused by a dfferental ncrease of the jont varable d d dfferental rotaton of the x y z frame n relaton to x - y - z - frame caused by a dfferental ncrease of the jont varable d d dfferental rotaton of the -th lnk frame n relaton to x - y - z - frame caused by a dfferental ncrease of the jont varable d dq jont dfferental matrx E homogeneous matrx descrbng the grpper frame n relaton to the end-effector frame

16 8 T.Szkodny versor of x - axs descrbed n --st frame J end-effector Jacoban J descrbed n the -th lnk frame j versor of y - axs descrbed n --st frame k versor of z - axs descrbed n --st frame l dsplacement along x axs T homogeneous matrx descrbng the end-effector frame n relaton to base frame T req requred matrx T T homogeneous matrx descrbng the end-effector frame n relaton to --st frame X homogeneous matrx descrbng the grpper frame n relaton to base frame X req requred matrx X x yz base frame x yz -th lnk frame x' yz ' ' frame connected wth -th lnk concdng wth --st frame for d 0 rotaton angle of the -th lnk frame n relaton to --st frame about the x axs Δ dfferental transformaton matrces of the end-effector descrbed n the -th frame 7 Δ 7 dfferental transformaton matrces of the grpper descrbed n the7-th frame rotaton angle of the -th lnk frame n relaton to the --st frame about z axs the -th lnk varable dsplacement along z axs References Chaccho P. Chavern S. and Sclano B. (99): Drect and nverse knematcs for coordnated moton tasks of two manpulator system. Journ. of Dynamcs Systems Measurement and Control. vol.8 No. pp Crag J.J. (989): Introducton to Robotcs. New York: Addson-Weseley Publ. Comp. charter. Jezersk E. (00): Dynamcs and Control of Robots. Warsaw: WNT 00 ch. (n Polsh). Kozłowsk K. Dutkewcz P. and Wróblewsk W. (00): Modelng and Control of Robots. Warsaw: PWN ch (n Polsh). Nakamura Y. and Hanafusa H.(009): Inverse knematc solutons wth sngularty robustness for manpulator control. Journ. of Dynamcs Systems Measurement and Control. vol.08 No. pp.-7. Sclano B. Scavcco L. Vllan L. and Orolo G. (00): Robotcs Modellng Plannng and Control. Sprnger Verlag Berln 00 chapter.. Spong M.W. and Vdyasagar M. (997): Robot Dynamcs and Control. Warsaw: WNT chapter..-. (n Polsh). Szkodny T. (009): Basc component of computatonal ntellgence for IRB-00 robots. Man-Machne Interactons. Berln Hedelberg: Sprnger-Verlag part.xi pp.7-. Szkodny T. (00): Inverse Knematcs Problem of IRB Fanuc Mtsubsh Adept and Kuka Seres Manpulators. J. of Appled Mechancs and Engneerng. Zelona Góra. Unversty Press vol. No. pp.87-8.

17 Dfferental knematcs of contemporary ndustral robots 9 Szkodny T. (0): Foundaton of Robotcs. Glwce. Slesan Unversty of Technology Publ. Company ch.. (n Polsh). Szkodny T. (0a): Knematcs of Industral Robots. Glwce. Slesan Unversty of Technology Publ. Company ch.8. (n Polsh). Szkodny T. (0b): Foundaton of Robotcs Problems Set. Glwce. Slesan Unversty of Technology Publ. Company ch. (n Polsh). Tchoń K. Mazur A. Dulęba I. Hossa R. and Muszyńsk R. (000): Moble Manpulators and Robots. Academc Publ. Company PLJ 000 chapter.. (n Polsh). Receved: May 0 Revsed: June 8 0

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

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

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

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

χ 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

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

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

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

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

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

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

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

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

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω

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

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

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

Numerical Solution of Ordinary Differential Equations

Numerical Solution of Ordinary Differential Equations Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples

More information

PHYS 705: Classical Mechanics. Canonical Transformation II

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

More information

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

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

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 Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach A Bayes Algorthm for the Multtask Pattern Recognton Problem Drect Approach Edward Puchala Wroclaw Unversty of Technology, Char of Systems and Computer etworks, Wybrzeze Wyspanskego 7, 50-370 Wroclaw, Poland

More information

6.3.4 Modified Euler s method of integration

6.3.4 Modified Euler s method of integration 6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from

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

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

This column is a continuation of our previous column

This column is a continuation of our previous column Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard

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

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent

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

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

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

More information

CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION

CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION Górnctwo Geonżynera Rok 0 Zeszyt / 006 Igor Konstantnovch Mladetskj * Petr Ivanovch Plov * Ekaterna Nkolaevna Kobets * Tasya Igorevna Markova * CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR

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

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

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

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

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

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

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

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

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

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

More information

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

DARWIN-OP HUMANOID ROBOT KINEMATICS. Robert L. Williams II, Ph.D., Mechanical Engineering, Ohio University, Athens, Ohio, USA

DARWIN-OP HUMANOID ROBOT KINEMATICS. Robert L. Williams II, Ph.D., Mechanical Engineering, Ohio University, Athens, Ohio, USA Proceedngs of the ASME 1 Internatonal Desgn Engneerng Techncal onferences & omputers and Informaton n Engneerng onference IDET/IE 1 August 1-15, 1, hcago, IL, USA DET1-765 DARWIN-OP HUMANOID ROBOT KINEMATIS

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

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

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

On the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros

On the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong

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

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

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

Modelling and Analysis of Planar Robotic Arm Dynamics Based on An Improved Transfer Matrix Method for Multi-body Systems 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

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

Chapter 3 Differentiation and Integration

Chapter 3 Differentiation and Integration MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton

More information

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra MEM 255 Introducton to Control Systems Revew: Bascs of Lnear Algebra Harry G. Kwatny Department of Mechancal Engneerng & Mechancs Drexel Unversty Outlne Vectors Matrces MATLAB Advanced Topcs Vectors A

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

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,

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

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

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

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

More information

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

Coordinate-Free Projective Geometry for Computer Vision

Coordinate-Free Projective Geometry for Computer Vision MM Research Preprnts,131 165 No. 18, Dec. 1999. Beng 131 Coordnate-Free Proectve Geometry for Computer Vson Hongbo L, Gerald Sommer 1. Introducton How to represent an mage pont algebracally? Gven a Cartesan

More information

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices Internatonal Mathematcal Forum, Vol 11, 2016, no 11, 513-520 HIKARI Ltd, wwwm-hkarcom http://dxdoorg/1012988/mf20166442 The Jacobsthal and Jacobsthal-Lucas Numbers va Square Roots of Matrces Saadet Arslan

More 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

An efficient algorithm for multivariate Maclaurin Newton transformation

An efficient algorithm for multivariate Maclaurin Newton transformation Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

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

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

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

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

Min Cut, Fast Cut, Polynomial Identities

Min Cut, Fast Cut, Polynomial Identities Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.

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

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

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

Integrals and Invariants of Euler-Lagrange Equations

Integrals and Invariants of Euler-Lagrange Equations Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,

More 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

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

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals

More 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

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 DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

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

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

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

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

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017 17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran

More 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

Robot Modeling and Kinematics Errata R. Manseur. List of errors and typos reported as of 4/1/2007:

Robot Modeling and Kinematics Errata R. Manseur. List of errors and typos reported as of 4/1/2007: Robot Modelng and Knematcs Errata R. Manseur. Lst of errors and typos reported as of 4//7:. Fle CF4.wrl mentoned on page 45 s not on the CD. Fle ZYZ_Sm.wrl mentoned on page 67 s not on the CD. 3. Pg. 43.

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

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

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to ChE Lecture Notes - D. Keer, 5/9/98 Lecture 6,7,8 - Rootndng n systems o equatons (A) Theory (B) Problems (C) MATLAB Applcatons Tet: Supplementary notes rom Instructor 6. Why s t mportant to be able to

More information

SIMPLE LINEAR REGRESSION

SIMPLE LINEAR REGRESSION Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two

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

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton

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

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

1 Convex Optimization

1 Convex Optimization Convex Optmzaton We wll consder convex optmzaton problems. Namely, mnmzaton problems where the objectve s convex (we assume no constrants for now). Such problems often arse n machne learnng. For example,

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

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

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

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