I nte rna t io na l J o urna l o f   Ro bo t ics a nd   Aut o m a t io n   ( I J RA )   Vo l.   5 ,   No .   4 Dec em b er   2 0 1 6 ,   p p .   27 7 ~ 28 3   I SS N:  2089 - 4856          277         J o ur na l   ho m ep a g e h ttp : //ia e s jo u r n a l.c o m/o n lin e/in d ex . p h p /I J RA   D y na m i c A na l y s i s   o f T w o  Link  R o bo t M a ni pu l a t o f o r C o ntro l  D es i g n   U s i n g  P I D  C o m puted T o rque C o ntr o l       J o lly   At it   Sh a h * ,   S S R a t t a n **   M e c h a n ica De p a rtme n t Na ti o n a in stit u te o f   tec h n o lo g y ,   In d ia   * *   M e c h a n ica De p a rt m e n t NIT ,   Ku ru k sh e tra ,   In d ia       Art icle  I nfo     AB ST RAC T   A r ticle  his to r y:   R ec eiv ed   Sep   24 ,   2 0 1 6   R ev i s ed   No v   10 ,   2 0 1 6   A cc ep ted   No v   26 ,   2 0 1 6       T h e   Du e   to   th e ir  a d v a n tag e   o h ig h   sp e e d ,   a c c u ra c y   a n d   re p e a tab il it y ,   ro b o m a n ip u lat o rs  h a v e   b e c o m e   m a jo c o m p o n e n o f   m a n u f a c tu rin g   in d u stries   a n d   e v e n   n o a   d a y th e y   b e c o m e   p a rt  o ro u ti n e   li f e .   Tw o   li n k   ro b o t   m a n ip u lat o is  a   v e ry   b a sic   c las si c a a n d   sim p le  e x a m p le  o f   ro b o f o ll o w e d   in   u n d e rsta n d i n g   o f   b a sic   f u n d a m e n tals  o f   r o b o ti c   m a n ip u lato r.   T h e   e q u a ti o n   o f   m o ti o n   f o tw o   li n k   ro b o t   is  a   n o n li n e a d if fe re n ti a e q u a ti o n .   F o h ig h e d e g re e o f   f re e d o m ,   a th e   c lo se d   f o rm   so lu ti o n a re   v e ry   d iff icu lt   w e   h a v e   to   u se   n u m e rica so lu ti o n .   He re   w e   f o c u se d   m a in l y   o n   c o n tr o l   o f   ro b o t   m a n ip u lat o to   g e th e   d e sire d   p o siti o n   u sin g   c o m b in a ti o n   o f   two   c las sic a l   m e th o d P ID  a n d   c o m p u ted   to rq u e   c o n tro m e th o d   a f ter  d e riv in g   th e   e q u a ti o n   o f   m o ti o n .   F o th e   sa m e   sim u latio n   is  re p re se n ted   u sin g   M AT LAB   a n d   c o m p a re d   w it h   c o m p u ted   t o r q u e   c o n tro l   m e th o d .     K ey w o r d :   C o m p u ted   to r q u co n tr o l   D y n a m ic  c o n tr o l   P I D   R o b o t m a n ip u lato r   Si m u latio n       Co p y rig h ©   201 6   In s t it u te o A d v a n c e d   E n g i n e e rin g   a n d   S c ien c e   Al rig h ts  re se rv e d .   C o r r e s p o nd ing   A uth o r :   J o lly   A ti t Sh a h ,   Natio n al  i n s ti tu te  o f   tech n o lo g y ,   I n d ia,   E m ail:  j o ll y s h a h 8 0 @ y a h o o . co m       1.   I NT RO D UCT I O N   I t is  w ell  k n o w n   th a t r o b o t m a n ip u lato r s   ar h ig h l y   n o n li n ea r ,   d y n a m icall y   co u p led   an d   ti m e - v ar y in g   s y s te m s   w h ic h   ar ex te n s i v e l y   u s ed   in   i n d u s tr ies.  T h r o b o tic  m a n ip u la to r s   ar g en er all y   s u b j ec ted   to   u n ce r tai n tie s   [ 1 ] .   B ec au s o f   th ese  u n ce r tai n tie s   an d   n o n l in e ar   b eh av io r ,   it' s   ch al len g i n g   task   to   co n tr o th e   m o ti o n   o f   r o b o m a n ip u lato r   at  ac cu r ate  p o s i tio n .   Fo r   th e   a cc u r ate  w o r k   an y   r o b o h a s   to   f o llo w   p r ed ef i n ed   tr aj ec to r ies  as  clo s as  p o s s ib le.   W h ich   m a k es  tr aj ec to r y   tr ac k in g   co n tr o m o s i m p o r tan task   i n   co n tr o o f   r o b o ti m an ip u lato r   [ 2 ] .   L o o k i n g   to   t h ab o v co n tr o r eq u ir e m en o f   r o b o tic  m an ip u lato r ,   r esear ch es   h av b ee n   r ep o r ted   o n   r o b o tic  co n tr o m e th o d s   s u ch   as   p r o p o r tio n al - in te g r atio n - d er i v ati v ( P I D)   co n tr o [ 3 ] , P co n tr o l,  P I   co n tr o l[ 3 ] ,   f ee d - f o r w ar d   co m p e n s atio n   co n tr o [ 4 ] ,   ad a p tiv co n tr o [ 5 ] ,   v ar iab le  s tr u ctu r co n tr o [ 6 ] ,   co m p u ted   to r q u co n tr o l[ 4 ] . I n   r ec en y ea r s ,   m u c h   atten tio n   h as  b ee n   p aid   to   th u s o f     P I c o n tr o an d   co m p u ted   to r q u co n tr o l in d iv id u all y   f o r   r o b o tic  m a n ip u lato r s .   T w o   d eg r ee   o f   f r ee d o m   r o b o m an ip u lato r   ac ts   l ik h u m a n   a r m s .   T h o u g h   h u m a n   ar m   is   7   DOF  b u t   t o   u n d er s ta n d   th co m p lex it y   o f   h u m a n   ar m s   it ' s   v er y   ea s y   to   tak ex a m p le  o f   t w o   d eg r e o f   f r ee d o m   r o b o t   m an ip u lato r   m o v e m e n t   [ 7 ]   in itiall y   an d   t h e n   attac h   g r ip p er   w it h   i to   m ak e   it  li k h u m an   ar m .   T h to tal   en er g y   w h ic h   i s   s u m   o f   k i n eti en er g y   a n d   p o ten tial  e n er g y   o f   th t w o   lin k   s y s te m   ar d ef i n ed   an d   u s ed   to   f o r m   L a g r an g ian   eq u atio n s .   Fi n all y   to   d ef i n t h to r q u ap p lied   o n   ea ch   li n k   th e s eq u atio n s   ar u s ed   [ 8 , 9 , 1 0 ]           Evaluation Warning : The document was created with Spire.PDF for Python.
                      I SS N :   2 0 8 9 - 4856     IJ RA    Vo l.   5 ,   No .   4 Dec em b er   2 0 1 6 :   27 7     28 3     278   2.   DYNA M I CS O F   T WO - L I N K   M ANIPULA T O RS   Her f o r   th an a l y s is   p u r p o s w h ad   u s ed   t w o   d eg r ee   o f   f r e ed o m   p la n er   r o b o t m a n ip u lato r   s h o w n   i n   Fig u r 1 L et  u s   in tr o d u ce   th f o llo w in g   n o tatio n s :   q -   T h j o in t a n g le  o f   th j o in t   i,   m i   -   T h m ass   o f   lin k   i,   l i   -   T h len g th   o f   lin k   i,   -   T h g r a v itatio n al  ac ce ler ati o n .           Fig u r e   1 .   T w o   Deg r ee   o f   Fre e d o m   R o b o t m a n ip u lato r           T h f o r w ar d   k i n e m at ics eq u at i o n s   f o r   2 - DO F r o b o tic  m a n ip u lato r   h av b ee n   d er iv ed   as   [ 8 ]                                           (           )                                     ( 1 )                                         (           )                                                     ( 2 )     I n   o r d er   to   s i m p li f y   w ass u m th at  t h ce n ter   o f   m a s s   o f   ea ch   li n k   is   r i g h at  t h m id d le  o f   t h li n k .   W w ill   also   a s s u m e   th a t t h g r av itatio n a l f o r ce   ac t s   i n   t h n eg ati v d ir ec tio n   o f   t h a x is .   L et  u s   ch o o s e   q = (   1 q 2 q ) T   as  th g e n er alize d   co o r d in ates  an d   d en o te  T   th k i n etic  en er g y ,   P   th p o ten tial  e n er g y   an d   L =T −P  th L a g r an g ian   f u n ctio n .   Usi n g   th e   n o tatio n s   ab o v e   an d   t ak in g   in to   ac co u n t   th e   ass u m p tio n   th a w h av e   m ad e,   th k in et ic  en er g y   a n d   p o ten tial e n er g y   o f   th s y s te m   c an   b ex p r ess ed   as f o llo w s :                  (           )         ̇                       ̇                     ̇     ̇        (           )              ( 3 )                                           (                             )                                 ( 4 )     T h eq u atio n s   o f   m o tio n   f o r   th is   m a n ip u lato r   ar e ,            (          ̇ )                         w h er i =   1 ,   2                                               ( 5 )     Su b s ti tu t in g   P   an d   T   in   th f o r m u la  o f   L   th d y n a m ic s   eq u ati o n   o f   th t w o - li n k   r o b o m a n i p u lato r   [ 9 ] :                     (     ̈       ̈   )                       (       ̈       ̈   )   (           )           ̈                     ̇                             ̇     ̇                        (           )                     ̇                                                                   ( 6 )                         (     ̈       ̈   )                       (     ̈ )                     ̈                              ̇                      ( 7 )     w h er c 1 co s ( q 1 )   ,   s 1   =sin ( q 1 ) ,   s 12 s in ( q 1 +q 2 ) ,   c 2 co s ( q 2 )   ,   s 2   =si n ( q 2 )     Evaluation Warning : The document was created with Spire.PDF for Python.
IJ RA   I SS N:  2089 - 4856         Dyn a mic  A n a lysi s   o f Tw o   Lin R o b o t Ma n i p u la to r   fo r   C o n tr o l D esig n   Usi n g   P I C o mp u te d   ( Jo lly  A tit S h a h )     279   W h ich   ca n   b r e w r itte n   as   [ 1 0 ] ,       (   )   ̈     (       ̇ )   ̇     (   )                                       ( 8 )     T h p ar am eter s   v i   an d   k i   ar v i s co u s   a n d   C o u lo m b   f r ictio n   co ef f icien ts ,   r esp ec ti v el y .     T h eq u atio n   ca n   b w r it ten   i n   th f o r m   o f               ̈       ̇                                                                    ( 9 )       w h er e       [                     ]                                                               (     )                                          (     )                                                [         ̇         ̇         ̇       ̇   ] ,   w it h                              (     )         [         ]   [         ̇         ̇ ]       [         ]   [               (     )                   (           )           (     )               (           ) ]             *           +   [          (     ̇ )          (     ̇ ) ]       3.   DE S I G ST E P S F O P I CO M P UT E T O RQ UE   CO NT RO L   O F   RO B O T   M ANIPULAT O R S   An y   m ec h an ical  s y s te m   h av i n g   tr an s latio n   m o tio n   i s   d escr ib ed   as b elo w   s h o w n   i n   F ig u r e   2 .           Fig u r e   2.   Me ch an ical  s y s te m   w it h   tr a n s latio n al  m o t io n       Op en   lo o p   d y n a m ic  eq u a tio n   f o r   s u ch   t y p o f   s y s te m   w ill b g iv e n   b y   f o llo w i n g   eq u a tio n ,                                                                         ( 1 0 )     ste p  1   Ass u m in g   t h at  w d o n ' t k n o w   v alu o f   m ,   b ,   k .   I n   m o d el  b ased   co n tr o m eth o d   m ,   b ,   k   ar ca lcu la ted   u s i n g   r ea ti m ca lc u lat io n s .   T h u s   s y s te m   w i ll  b e   s i m p li f ied   to   m a s s   o n l y   s y s te m .   A p p l y i n g   m o d el  b ased   co n tr o la w   to   eq u atio n   ( 1 0 )   b ec o m e,     Evaluation Warning : The document was created with Spire.PDF for Python.
                      I SS N :   2 0 8 9 - 4856     IJ RA    Vo l.   5 ,   No .   4 Dec em b er   2 0 1 6 :   27 7     28 3     280                            w h er               ,   w h ic h   is   n e w   in p u t to   s y s te m   a s   s h o w n   i n   F ig u r e   3   o n   n e x t p ag e         Fig u r e   3 .   Stru ct u r o f   co m p u t er   to r q u co n tr o l       Ste p 2   A p p l y   a n y   cla s s ical  co n tr o l la w   li k P D,   P I ,   P I D.   I f   w ap p l y   s er v o   co n tr o l la w   t o   ab o v eq u atio n ,   it  w i ll g et  c o n v er ted   in to   co m p u ted   to r q u co n tr o l s y s te m .   T o   im p r o v t h ac c u r ac y   w w il u s f ee d b ac k   co n tr o l.  C o m p u ted   to r q u co n tr o alo n g   w it h   P I f ee d b ac k   co n tr o is   k n o w n   a s   P I co m p u ted   to r q u co n tr o m et h o d .   T o tal  s y s te m   b lo ck   d iag r a m   is   s h o w n   b elo w   in   Fig u r e   4 .         Fig u r e   4 .   P I f ee d b ac k   lo o p   c o m p u ted   to r q u co n tr o ller       th u s   to tal  co n tr o l la w   is   d e f in e d   as,           (   ) (   ̈           ̇                 )     (       ̇ )     (   )     w h er   (   )   is   i n teg r al  o f   tr ac k i n g   e r r o r   e( t )   an d     k d , k p , k ar co n tr o l g ain .       4.   SI M UL AT I O A ND  RE SU L T S   Fo r   s i m u latio n   p u r p o s f o llo win g   v al u es a r tak e n   i n   M A T L A B   P r o g r a m m in g .   T o tal  Sim u latio n   T im T 1   =   1 0   s ec   an d   T 2 1 2 0   s ec ;   Sa m p li n g   T i m T   0 . 0 0 2 s e c ;   m as s     m   [ 2 ;2 ]   u n it;   l1   1   u n it;   l2   1   u n it;   g r av   9 . 8   u n it ;     Gain s :   Kd   [ 1 0 0 0   0 ;0   1 0 0 0 ] ;   Kp   [ 1 0 0   0 ;0   1 0 0 ] ;   Ki  [ 1 0   0 ; 0   1 0 ] ;     Evaluation Warning : The document was created with Spire.PDF for Python.
IJ RA   I SS N:  2089 - 4856         Dyn a mic  A n a lysi s   o f Tw o   Lin R o b o t Ma n i p u la to r   fo r   C o n tr o l D esig n   Usi n g   P I C o mp u te d   ( Jo lly  A tit S h a h )     281   S i m u latio n   is   d o n to   co n f ir m   th r esu lt s   o f   th p r o p o s ed   c o n tr o s ch e m e.   Si m u la tio n   is   d o n f o r   f ir s 1 0   s ec   an d   th e n   f o r   s ec o n d   ti m 1 2 0 s ec ,   an d   at  ev er y   0 . 0 0 2 s ec   th s a m p le  r ea d in g   w a s   tak e n   to   p lo t th g r ap h s .     Fig u r s h o w s   t h p o s itio n   er r o r s   f o r   jo in t 1   an d   j o in t 2   f o r   1 0   s ec   an d   F ig . 6   s h o w s   t h v el o cit y   er r o r s   f o r   j o in t 1   an d   j o in t 2   f o r   1 0   s ec .           Fig u r e   5 .   P o s itio n   er r o r s   f o r   j o in t 1   an d   j o in t 2   f o r   1 0   s ec     Fig u r e   6 .   V elo cit y   er r o r s   f o r   jo in t 1   an d   j o in t 2     f o r   1 0   s ec       Fig u r 7   an d   F i g u r 8   s h o ws  th e   co m p ar i s o n   f o r   d er iv e d   an d   p r ed icted   p o s itio n   f o r   lin k   1   an d   lin k   2   r esp ec tiv el y   f o r   1 0   s ec .   F ig u r e   s h o w s   t h to r q u r eq u ir ed   f o r   j o in 1   a n d   j o in 2   f o r   co n tr o p u r p o s e   f o r   1 0   s ec .         Fig u r e   9 .   T o r q u r eq u ir ed   f o r   jo in t 1   an d   j o in t 2   f o r   1 0   s ec .       Fig u r 10  s h o w s   t h p o s itio n   er r o r s   f o r   j o in 1   an d   jo in 2   f o r   1 2 0   s ec   an d   F ig u r 11  s h o w s   th v elo cit y   er r o r s   f o r   j o in t 1   an d   j o in t 2   f o r   1 2 0   s ec .           Fig u r e   10 .   P o s itio n   er r o r s   f o r   jo in t 1   an d   j o in t 2     f o r   1 2 0   s ec     Fig u r e   11 .   V elo city   er r o r s   f o r   j o in t 1   an d   j o in t 2     f o r   1 2 0   s ec   Evaluation Warning : The document was created with Spire.PDF for Python.
                      I SS N :   2 0 8 9 - 4856     IJ RA    Vo l.   5 ,   No .   4 Dec em b er   2 0 1 6 :   27 7     28 3     282   Fig u r e 1 2   an d   F i g u r e 13   s h o ws  th e   co m p ar is o n   f o r   d er i v ed   an d   p r e d icted   p o s itio n   f o r   lin k   1   an d   lin k   2   r esp ec tiv el y   f o r   1 2 0   s ec .           Fig u r e   12 .   C o m p ar is o n   o f   d er i v ed   an d   p r ed icted   p o s itio n   f o r   lin k   1   f o r   1 2 0   s ec .     Fig u r e   13 .   C o m p ar is o n   o f   d er i v ed   an d   p r ed icted   p o s itio n   f o r   lin k   2   f o r   1 2 0   s ec .       F ig u r e   14  s h o w s   t h to r q u r eq u ir ed   f o r   j o in t 1   an d   j o in t 2   f o r   co n tr o l p u r p o s e   f o r   1 2 0   s ec .         Fig u r e   14 .   T o r q u r eq u ir ed   f o r   j o in t 1   an d   j o in t 2   f o r   1 2 0   s ec .       5.   CO NCLU SI O N   I n   th i s   p ap er   w ai m ed   to   in tr o d u ce   n e w   co n tr o ller   s y s t e m   w it h   th h elp   o f   t w o   co n v en t io n al   co n tr o ller   P I an d   C o m p u ter   t o r q u co n tr o l.  T h s i m u latio n   r esu lt s   s h o w   th e   v a lid it y   o f   t h p r o p o s ed   m et h o d   an d   g i v th p o s s ib ili t y   o f   P I co m p u ted   to r q u co n t r o f o r   2   DOF  r o b o m an ip u lato r .   P o s itio n   er r o r s   an d   Velo cit y   er r o r s   ar co m in g   to   b less   th an   1 %.  A n d   ev e n   a f ter   f e w   s ec o n d s   th o s ar g et tin g   al m o s ar o u n d   ze r o   as  s h o w n   i n   s i m u latio n   g r ap h s .   Si m u latio n s   g r ap h s   als o   s h o w s   th a er r o r   b et w ee n   d e r iv ed   an d   p r ed icted   p o s itio n   f o r   lin k   1   an d   li n k   2   i s   al m o s t a r o u n d   ze r o .       RE F E R E NC E   [ 1 ]   Am o A .   K h alate ,   Go p i n ath a n   L ee n a,   Go s h aid as   R a y ,   A n   A d a p tive  F u z z C o n tr o ller   f o r   Tr a jecto r y   Tr a ck in g   o f R o b o t Ma n ip u l a to r ,   I n telli g en t Co n tr o l a n d   Au to m atio n ,   2 0 1 1 ,   3 6 4 - 370   [ 2 ]   Ko r en   W ar d ,   R a p id   S imu lta n eo u s   Lea r n in g   o Mu ltip le  B e h a vio u r s   w ith   a   Mo b ile  R o b o t ,   P r o c.   2 0 0 1 ,   Au s tr alia n   C o n f er en ce   o n   R o b o ics an d   Au to m atio n ,   S y d n e y ,   14 - 1 5   No v   2 0 0 1   [ 3 ]   J . Y. S.  L u h ,   C o n v e n tio n al  co n tr o ller   d esig n   f o r   I n d u s tr ial  r o b o ts - tu to r ial,   I E E E   T r an s .   S y s t e m s   Ma n   C y b er n et.   1 9 8 3 ,   2 9 8 3 1 6 .   Evaluation Warning : The document was created with Spire.PDF for Python.
IJ RA   I SS N:  2089 - 4856         Dyn a mic  A n a lysi s   o f Tw o   Lin R o b o t Ma n i p u la to r   fo r   C o n tr o l D esig n   Usi n g   P I C o mp u te d   ( Jo lly  A tit S h a h )     283   [ 4 ]   P r ad ee p   K.   Kh o s la,   T ak eo   Ka n ad e,   E x p er i m e n tal  o f   E v al u at io n   o f   No n l in ea r   Feed b ac k   a n d   Feed f o r w ar d   C o n tr o Sc h e m es   f o r   Ma n ip u lato r s ,   Th I n tern a tio n a Jo u r n a o R o b o tir s   R esea r ch ,   Vo l.  7 ,   No .   1 ,   Feb r u ar y   1 9 8 8 .   [ 5 ]   J . J . E .   Slo tin e,   W .   L i,  A d ap tiv Ma n ip u lato r   C o n tr o l:  A   C a s Stu d y ,   I E E E   T r an s .   Au to m a t.  C o n tr o 3 3 ,   1 9 8 8   ,   9 9 5 1 0 0 3   [ 6 ]   Ma r k   W . Sp o n g   a n d   M. Vid y a s ag ar ,   R o b o t D yn a mics a n d   co n tr o l ,   J o h n   W ile y   &   So n s , 2 0 0 4 .   [ 7 ]   J o h n   Ho ller b ac h ,   C o mp u ter,  B r a in s   a n d   th e   C o n tr o Mo ve men t ,   Ma n ch es ter   I n s tit u te  o f   T ec h n o lo g y ,   A r ti f icial  I n telli g e n ce   L ab o r at o r y ,     J u n 1 9 8 2   [ 8 ]   J o lly   S h ah ,   Dr .   S.S.Ratta n ,   P r o f .   B .   C .   Nak r a,   Kin e m atic  An al y s i s   o f   P lan er   R o b o t   Usi n g   A r ti f icial   Neu r al  Net w o r k I n tern a tio n a l J o u r n a l o R o b o tics   a n d   A u to ma tio n ,   Sep tem b e r   2 0 1 2 ,   1 4 5 - 151.   [ 9 ]   Step h en   Do d d s ,   A   t w o - li n k   m an ip u lato r s i m u latio n   a n d   co n tr o d esig n   A d v i s o r Dr .   Ma h b o u b   B ac co u c h,   Un i v er s it y   o f   Neb r ask at  O m a h O m a h a,   NE ,   6 8 1 3 4 ,   Ma r ch   2 0 1 2         [ 1 0 ]   Du y   N g u y e n - T u o n g ,   J an   P eter s ,   Lea r n in g   R o b o t   Dyn a mics  fo r   C o mp u ted   To r q u C o n tr o u s in g   Lo ca Ga u s s ia n   P r o ce s s es  R eg r es s i o n ,   E C SIS  S y m p o s i u m   o n   L ea r n in g   an d   A d ap tiv B eh a v i o r s   f o r   R o b o tic   S y s te m s ,   2 0 0 8     Evaluation Warning : The document was created with Spire.PDF for Python.