DIFFERENTIAL KINEMATICS OF CONTEMPORARY INDUSTRIAL ROBOTS
|
|
- Valentine Lang
- 6 years ago
- Views:
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 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 informationMEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut
MEV442 Introducton to Robotcs Module 2 Dr. Santhakumar Mohan Assstant Professor Mechancal Engneerng Natonal Insttute of Technology Calcut Jacobans: Veloctes and statc forces Introducton Notaton for tme-varyng
More informationChapter 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 informationIterative 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
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationCHAPTER 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 informationCSci 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 informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationStructure 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 informationSecond 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 informationTHE 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 informationNovember 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 informationIndeterminate 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 informationNMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582
NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 4. Moton Knematcs 4.2 Angular Velocty Knematcs Summary From the last lecture we concluded that: If the jonts
More informationNote 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 informationNumerical 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 informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationOn 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 informationA 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 information6.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 informationFinite 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 informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More informationThis 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 informationErrors 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 informationAPPENDIX 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 informationIrregular 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 informationSolution 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 informationCHARACTERISTICS 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 informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationOne-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 informationCOMPOSITE 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 information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationUncertainty 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 informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
More informationEEE 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 informationDARWIN-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 informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationOn 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 informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationA 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 informationModelling 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 informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationChapter 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 informationMEM 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 informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationA 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 informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationModule 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 informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationCoordinate-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 informationThe Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices
Internatonal Mathematcal Forum, Vol 11, 2016, no 11, 513-520 HIKARI Ltd, wwwm-hkarcom http://dxdoorg/1012988/mf20166442 The Jacobsthal and Jacobsthal-Lucas Numbers va Square Roots of Matrces Saadet Arslan
More informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not
More informationAn efficient algorithm for multivariate Maclaurin Newton transformation
Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,
More informationSome 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 informationMA 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 informationA 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 informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationMin 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 informationThe Analysis of Coriolis Effect on a Robot Manipulator
Internatonal Journal of Innovatons n Engneerng and echnology (IJIE) he Analyss of Corols Effect on a Robot Manpulator Pratap P homas Assstant Professor Department of Mechancal Engneerng K G Reddy college
More informationImplicit 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 informationChapter 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 informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationModeling of Dynamic Systems
Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how
More informationPhysics 111: Mechanics Lecture 11
Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton
More informationTensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q
For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationLectures - 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 informationON 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 informationLecture 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 informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationA 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 informationChapter - 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 informationKinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017
17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran
More informationStudy 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 informationRobot 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 informationNumerical 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 informationEN40: 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 information36.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 informationSIMPLE 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 informationOn 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 informationAdvanced 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 informationModule 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 informationLab 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 information1 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 informationNMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582
NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationChapter 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 informationNotes on Analytical Dynamics
Notes on Analytcal Dynamcs Jan Peters & Mchael Mstry October 7, 004 Newtonan Mechancs Basc Asssumptons and Newtons Laws Lonely pontmasses wth postve mass Newtons st: Constant velocty v n an nertal frame
More information