I nte rna t io na l J o urna l o f   E lect rica l a nd   Co m p ute E ng in ee ring   ( I J E CE )   Vo l.   7 ,   No .   3 J u n e   201 7 ,   p p .   1506 ~ 1 5 1 2   I SS N:  2088 - 8708 DOI : 1 0 . 1 1 5 9 1 / i j ec e . v 7 i 3 . p p 1 5 0 6 - 1512           1506       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 JE C E   M o deling   the  Vib ra tiona l Dyna m ic s o Pie z o el ectr i Actuator  by   Sy ste m  Iden tif ica tion Techniq ue       Nurul B a hia h M o hd   No o r 1 ,   M o hd   Ridzua n Ah m a d 2   1, 2   De p a rtm e n o f   Co n tro a n d   M e c h a tro n ics   E n g in e e rin g ,   F a c u lt y   o f   El e c tri c a En g in e e rin g ,     Un iv e rsiti   T e k n o lo g M a lay sia       Art icle  I nfo     AB ST RAC T   A r ticle  his to r y:   R ec eiv ed   Mar   4 ,   2 0 1 7   R ev i s ed   Ma y   7 ,   2 0 1 7   A cc ep ted   Ma y   26 ,   2 0 1 7       A c tu a to rs  b a se d   o n   sm a rt  m a teri a ls  su c h   a p iez o e lec tri c   a c tu a to rs  (P EA s)  a re   w id e l y   u se d   in   m a n y   a p p li c a ti o n to   tran sf o rm   e le c tri c a l   sig n a to   m e c h a n ica sig n a a n d   v ice   v e rs a .   Ho w e v e r,   th e   m a jo d ra w b a c k f o th e se   s m a rt  a c tu a to rs  a re   h y ste r e sis  n o n li n e a r,   c re e p   a n d   re si d u a v ib ra t io n .   In   th is   p a p e r,   P EA a re   u se d   f o a c ti v e   v ib ra ti o n   a p p l ica ti o n .   T h e re f o re ,   a   m o d e o P EA   m u st  b e   e sta b li sh e d   t o   c o n t ro t h e   v ib ra ti o n   th a t   o c c u rs  i n   t h e   sy st e m .   T h e   f re q u e n c ies   o f   1   Hz ,   2 0   Hz   a n d   5 0   Hz   w e re   tes ted   on   th e   P EA s.  T h e   re su lt s o b tai n e d   f ro m   th e   e x p e ri m e n t a w e r e   u se d   to   d e v e lo p   tran sf e f u n c ti o n   m o d e b y   e m p lo y in g   s y ste m   id e n ti c a ti o n   tec h n i q u e .   M e a n w h il e ,   th e   m o d e v a li d a ti o n   w a s   b a se d   o n   lev e o f   m o d e ls  tn e ss   to   e sti m a ti o n   d a ta,  m e a n   sq u a re d   e rro (M S E),   n a p re d i c ti o n   e rro (F P E)  a n d   c o rre latio n   tes t .   T h e   e x p e ri m e n tal  re su lt   sh o we d   th a th e   d isp lac e m e n o f   th e   a c tu a to is  in v e rse l y   p ro p o rti o n a t o   t h e   f re q u e n c y .   Th e   f o ll o w in g   c o n se q u e n c e c a u se d   th e   ti m e   re sp o n se   c rit e ria at  5 0   Hz   a c h iev e d   sm a ll e st o v e r sh o o a n d   f a ste st r e sp o n se   o f   rise   ti m e   a n d   se tt li n g   ti m e .   K ey w o r d :   A ct u ato r   P iezo elec tr ic   S y s te m   id e n ti f icatio n   Vib r atio n   Co p y rig h ©   2 0 1 7   In stit 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 :   Nu r u l B a h iah   Mo h d   No o r ,     Dep ar t m en t o f   C o n tr o l a n d   Me ch atr o n ics E n g in ee r i n g ,     Facu lt y   o f   E lectr ical  E n g in ee r in g ,   Un i v er s iti T ek n o lo g i M ala y s ia .   E m ail:  n u r u lb ah ia h m n o o r @ g m ail. co m       1.   I NT RO D UCT I O N   P iezo elec tr icity   f r o m   th Gr ee k   w o r d ,   p iezo   o r   p iezin   m ea n s   p r ess u r elec tr icit y   w a s   r s d is co v er ed   b y   J .   an d   P .   C u r ie  in   1 8 8 0   in   q u ar tz  cr y s tal.   P iezo elec tr icit y   i s   f u n d a m e n tal  p r o ce s s   in   elec tr o m ec h an ical   en er g y   co n v er s io n   i n   ce r tain   cr y s talli n m ater ials   [ 1 ] .   T h ese  cr y s talli n m ater ial s   tu r n   o u to   b elec tr ically   p o lar ized   w h e n   s u b j ec ted   to   m ec h a n ical   s tr es s   a n d   co n v er s el y   c h an g s h ap i n   t h p r ese n ce   o f   elec tr ic  eld .   So o n   a f ter ,   it  w as   d is co v er ed   th at  m ater ial s   s h o w ed   d e f o r m atio n   a s   elec tr ical   s i g n al  i s   ap p lied   p r esen ted   f o r   ac tu ato r s   m ater ia l.  P iezo elec tr ic  ac tu ato r s   ar s m ar m a ter ial  d ev ices  b ec o m co m p etiti v e   ch o i ce s   f o r   m icr o   ap p licatio n s   [ 2 ] .   P iezo elec tr i ac tu ato r s   h av b ee n   i n cr ea s in g l y   u s ed   b ec au s o f   h i g h   s tiff n es s ,   h i g h   r eso lu tio n ,   f a s r esp o n s e,   n o   m ag n etic   eld   a n d   h i g h   p r ec is io n   [ 3 ] - [ 4 ] .   A s   an   e x a m p le,   p i ez o elec tr ic  ac tu ato r   is   u s ed   in   ato m ic  f o r ce   m icr o s co p y   ( A FM)   a n d   s ca n n i n g   p r o b m icr o s co p y   ( SP M)   f o r   m i cr o   an d   n an o s ca le s   ap p licatio n s   [ 3 ] .   T h ey   ar al s o   w id el y   u s ed   i n   m ec h atr o n i ap p licatio n s   s u c h   as   ca r   f u el  in j ec to r ,   m ed ical  s u r g er y   eq u ip m e n a n d   ac tiv e   v ib r atio n   co n tr o [ 5 ] .   Nev er th eles s ,   p iezo elec tr ic  ac tu ato r s   r ev ea u n d es ir ab l e   d is p lace m e n t b eh a v io r   th at  li m it s   th ac c u r ac y   a n d   s tab ilit y   o f   th ac t u ato r s .   T h m aj o r   is s u e s   ar e:   i.   H y s ter e s is :   T h s y s te m   h as  d i f f er en o u tp u w it h   t h s a m i n p u ap p lied   to   p iezo elec tr ic  ac tu ato r .   T h er r o r   g en er ates  1 0 - 15   %   er r o r   u n d er   o p en - lo o p   v o lta g co n tr o l.  T h o u g h ,   th e   er r o r   r ea ch es  3 5   w h e n   th f r eq u e n c y   o f   t h in p u t si g n al  is   in cr ea s ed   [ 6 ] .   ii.   C r ee p : T h am p lit u d o f   d is p l ac e m en t o u tp u t i s   ch a n g o v er   ti m at  co n s tan t in p u t v o lta g e   s ig n al  [ 7 ] .     Evaluation Warning : The document was created with Spire.PDF for Python.
I J E C E     I SS N:  2088 - 8708       Mo d elin g   th V ib r a tio n a l D yn a mics o f P iezo elec tr ic  A ctu a to r   b y   S ystem  ....   ( N u r u l B a h ia h   Mo h d   N o o r )   1507   iii.   R esid u al  v ib r atio n T h ac tu ato r   v ib r ates  b ef o r r ea ch   s tead y   s tate  o r   fin al  p o s itio n .   T h is   is   d u to   m ec h a n ical  p r o p er ties   s u c h   as  m as s ,   s ti f f n e s s   a n d   d a m p in g   [ 8 ] .   T h h y s ter esis   b e h av io r   i n   p iezo elec tr ic  ac tu ato r   ca u s e s   i n ac cu r ac y   an d   i n s tab ilit y   to   its   ap p lied   s y s te m .   H y s ter e s is   is   c h ar ac t er ized   as  n o n lin ea r   b ec au s t h v o lta g -   d is p lace m e n c u r v is   n o o n e - o n m ap p in g   [ 3 ] , [ 9 ] .   T h s a m e   in p u v o ltag e   p r o d u ce s   d i f f er en t   o u tp u d is p lace m e n f o r   in cr ea s in g   a n d   d ec r ea s in g   i n p u v o ltag e.   I n   a d d itio n ,   h y s ter esis   b eh a v io r   al s o   ch ar ac ter ized   as  r ate - d ep en d en c y   b ec au s t h e   h y s ter esi s   cu r v o r ien tatio n   a n d   s h ap is   ch a n g ed   as th f r eq u en c y   o f   i n p u v o lta g is   ch a n g ed .     Fig u r e   1   s h o w ed   t h h y s ter es is   n o n l in ea r   p h e n o m e n o n   in   t h p iezo elec tr ic  ac tu ato r s   f o r   o p en   lo o p   an d   clo s lo o p   s y s te m .   T h c lo s lo o p   s y s te m   s h o w s   th e   h y s ter es is   e f f ec f o r   t h ac t u at o r   is   m in i m ized   b y   u s i n g   co n tr o ller .   Fi g u r e   2   s h o w ed   e f f ec o f   cr ee p   p h e n o m en in   p iezo elec tr ic  ac t u ato r .   T h is   p h e n o m e n o n   w il l   g iv s ig n i fican er r o r   w h en   t h P E A   ex te n d ed   p er io d s   o f   ti m e.   I ca n   b o b s er v ed   th at,   at  th s a m i n p u t   v o ltag s i g n al,   th d is p lace m en o u tp u i s   ch a n g ed   o v er ti m e.   A cc o r d in g   to   [ 1 0 ] ,   th h y s ter es is   an d   cr ee p   s tr o n g l y   af f ec t   t h ac t u ato r   at   f r eq u e n c y   lo w er   t h a n   1 0   Hz .   Mo r eo v er ,   v ib r atio n   is   r ep etitiv e   m o tio n s   th a t   o cc u r   in   ti m i n ter v al  w it h   r es p ec to   p o in o f   r ef er en ce .   T h v ib r atio n   is   ca u s e d   b y   p o te n tial  e n er g y   s to r ag e   in   s p r in g ,   k in et ic  en er g y   s t o r ag f o r   m a s s   a n d   g r ad u al  e n er g y   d a m p i n g   f o r   d a m p er .   T h ese  ef f ec ts   r ed u ce d   th e   p er f o r m a n ce   ac cu r ac y   an d   i n cr ea s e d   th u n ce r tai n t y   o f   m ea s u r e m e n ts   [ 1 1 ] .   Mo r eo v er ,   v ib r atio n   is   r ep etitiv m o tio n   t h at  o cc u r   i n   ti m i n ter v a w ith   r esp ec to   p o in o f   r ef er en ce .   T h v ib r atio n   o cc u r r ed   b ec au s o f   p o ten tial  en er g y   s t o r ag in   s p r in g ,   k in e tic  en er g y   s to r ag f o r   m a s s   an d   g r ad u al  en er g y   d a m p i n g   f o r   d a m p er .   T h is   ef f ec r ed u c ed   th p er f o r m a n ce   ac cu r ac y   an d   in cr ea s e d   t h u n ce r tai n t y   o f   m ea s u r e m en t s   [ 1 0 ] .   Fu r th er m o r e,   th v ib r at io n   f r eq u e n c y   r estra i n ed   th d ev elo p m en t   o f   h i g h   s p ee d - SPM   to   be   u s ed   in   s e m ico n d u cto r   i n d u s tr ie s ,   to   s tu d y   t h e   m ec h a n is m   o f   f a s s u r f ac p r o ce s s es,  a n d   to   d ev elo p ed   SP M - b ased   n an o f ab r icatio n   eq u ip m en ts           Fig u r e   1 .   H y s ter es is   n o n li n ea r   b eh av i o u r   i n   o p en   lo o p   an d   clo s lo o p   s y s te m           Fig u r e   2 .   C r ee p   p r o p er ties   o v er   ti m e   Evaluation Warning : The document was created with Spire.PDF for Python.
                      I SS N :   2 0 8 8 - 8708   I J E C E    Vo l.  7 ,   No .   3 J u n e   2 0 1 7   :   1 5 0 6     1 5 1 2   1508   T h m at h e m atica m o d el  o f   P E A   ca n   b d ev elo p ed   b y   u s i n g   p h y s ic s   la w   an d   p ar a m eter   esti m atio n   ap p r o ac h .   T h p h y s ics  la w   n ee d   f u l u n d er s ta n d in g   o f   th e   s y s te m   an d   p h y s ical  i n ter ac tio n s   r elate d   to   it .   Me an w h ile,   p ar a m eter   esti m atio n   ap p r o ac h   o n l y   r eq u ir ed   in p ut - o u tp u d ata  [ 1 2 ] .   T h p ar a m eter s   o f   m at h e m a tical  m o d el  th at  o b tain ed   b ased   o n   p h y s ics  la w   t h en   n ee d   to   b o p tim ized .   A cc o r d in g   to   [ 1 3 ] ,   th e   p ar am eter s   o f   m o d el  w a s   o b tain ed   b y   u s i n g   ar tif ic ial  n e u r al   n et w o r k   ( ANN) .   T h n eu r al  n et w o r k s   h ad   to   b e   tr ain ed   u n til  t h s i m p le   m o d el  w a s   ac h iev ed .   B esid es,  t h f aster   tr ai n i n g   o f   s i m p le  m o d el  r esu lti n g   le s s   co m p u tatio n   in   ac tu al  co n tr o l .   T h o u g h ,   t h n e u r al  n et w o r k   n ee d   to   b cr ea ted ,   tr ain ed ,   an d   ev al u ated   f o r   v ar io u s   f ee d f o r w ar d   b ac k p r o p ag atio n   n et w o r k s   u n til  t h e   n et w o r k   g av ac c u r ate  m o d e l.  Fu r t h er m o r e,   th e   p ar am etr ic  id en tific atio n   s u c h   as  r ec u r s i v least  s q u ar m et h o d ,   g en etic  al g o r ith m   a n d   p ar ticle  s w ar m   o p tim izatio n   also   ca n   b u s ed   f o r   p ar am eter s   o p ti m izatio n .   T h g en etic  alg o r it h m   ( G A )   is   h eu r is tic  m et h o d   to   esti m ate  t h p ar a m eter s   b as ed   o n   ev o lu tio n ar y   id ea s   o f   n a tu r al  s elec tio n   an d   g e n etic s   [ 1 4 ] .   T h p ar am eter s   th at  n ee d   to   b r ec o g n ized   k n o w n   as  c h r o m o s o m e s   t h at  f o r m   s tr i n g   ca lled   g e n e.   T h e   p r o ce s s   of   GA   wa s   s tar ted   w it h   i n itialized   th v al u es  o f   c h r o m o s o m es.  T h v al u es  we r r an d o m l y   c h o s e n   to   f o r m   a   p o p u latio n .   T h en ,   th e v al u atio n   t n es s   a n d   alg o r it h m   w er e   ca lc u lated   u n t il  t h b est  e s ti m atio n   v a lu e   wa s   ac co m p li s h ed .   T h ad v an ta g f o r   th i s   id e n ti ca t io n   o p ti m izatio n   tec h n iq u is ,   t h w r o n g   i n it ializat io n   v alu e   f o r   esti m atio n   p ar am eter s   ca n   b id en ti fied   s in ce   t h o p ti m a p ar a m eter s   lies   o u t s id in ter v al.   Hen ce ,   r esu lt s   ten d   to   b e   clu s ter ed   n ea r   t h s id o f   t h in ter v al  t h at  s h o u ld   b tu n ed .   Nev er th eles s ,   G A   m e th o d   ca n n o k ee p   th e   m e m o r y   w h e n   t h p o p u lat io n   w a s   ch a n g ed   [ 15 ].   T h er ef o r e,   r esear ch er   [ 1 6 ]   w as  i m p le m en tin g   p ar ticle  s w ar m   o p tim izatio n   ( P SO)   tech n iq u e .   P SO  is   k n o w n   a s   n e w   h e u r is tic   m et h o d   w h ic h   is   ea s y   to   u n d er s ta n d   an d   i m p le m en to   o p ti m ize  t h p ar am eter s .   I n   ad d itio n ,   P SO  p r o v id ed   ad v an tag e s   w it h   h ig h   o p ti m izat io n   p er f o r m a n ce   a n d   g o o d   m e m o r y   w h er ea c h   p ar ticle  in   P SO  wa s   p r ep ar ed   w it h   s m all  m e m o r y   co m p r is i n g   its   p r ev io u s   b es p o s itio n .   T h is   t ec h n iq u wa s   b ased   o n   th e   i d ea   o f   co llab o r ativ b eh a v io u r   an d   s w ar m i n g   i n   b io lo g ical  p o p u latio n   [ 17 ] .   T h p ar ticles  h o ld   p o s itio n   ( x j )   an d   v elo cit y   ( v j )   w it h   tr i al  s o lu tio n   o f   t h e   o p tim u m   p r o b lem .   Ho w e v er ,   th e   m ai n   d r a w b ac k   f o r   t h is   tec h n iq u wa s   t h t n ess   f u n ct io n   m u s b s elec ted   b ased   o n   s u itab ilit y   o f   m o d els   an d   PS h as  s lo w   co n v er g en ce   s p ee d   [ 1 8 ] .   Nex t,  p ar a m eter   esti m atio n   ap p r o ac h   s u ch   a s   s y s te m   id en ti f icatio n   ( SI)   is   tec h n iq u to   o b tain ed   th m at h e m a tic al  d escr ip tio n   f r o m   m ea s u r e m e n t s   o f   th e   in p u a n d   o u tp u t.  T h is   tec h n iq u o n l y   r eq u ir ed   d ata  co llectio n   o f   i n p ut - o u tp u d ata  to   f i n d   th ex ac o r   ap p r o x i m ate   m o d els  o f   d y n a m ic  s y s te m   [ 1 9 ] .     T h o b j ec tiv o f   SI  is   t o   f in d   th e x ac o r   ap p r o x im a te  m o d els o f   d y n a m ic  s y s te m   b ased   o n   k n o w led g e   o f   th o b s er v ed   in p u t a n d   o u t p u t d ata.       2.   M E T H O D O L O G Y   T h m ain   co m p o n e n o f   t h e   ex p er i m e n s y s te m   u s ed   i n   th i s   s tu d y   w as   p iezo elec tr ic  ac tu ato r   ( A E 0 5 0 5 D1 6 DF NE C   T o k in ) .   T h is   ac t u ato r   p r o v id ed   o u tp u d is p lace m en co r r esp o n d in g   o f   1 7 . 4   ±   2 . 0   µ w it h   1 5 0   VDC  m a x i m u m   v o lt ag e.   Fi g u r e   1   s h o w ed   th e x p er i m en tal  s et u p   th a t h a s   b ee n   u s ed   as a   p latf o r m   to   m o d el  t h P E A .   T h f r a m o f   tip   P E A   m o v e m e n wa s   ca p t u r ed   b y   a   s o f t w ar t h at  wa s   c o n n ec ted   to   in v er ted   m icr o s co p e.   Me an w h ile,   t h m icr o co n tr o ller   an d   s w itc h in g   cir cu it  wa s   co n n ec ted   to   t h l ap to p   to   co n tr o th e   f r eq u en c y   o f   P E A .           Fig u r e   3.    P latf o r m   f o r   co n d u c tin g   th e x p er i m e n ts       T h v ib r atio n   o f   P E A   wa s   c o n tr o l led   b y   th f r eq u e n c y   t h at  h a s   b ee n   g en er ated   b y   A r d u in o   U n o   m icr o co n tr o ller   an d   s w itc h i n g   cir cu it .   T h s q u ar w a v e   in p u t   s ig n al  wa s   s u p p lied   to   th P E A   w it h   5 0   d u t y   c y cle.   T h s w itc h i n g   cir cu i u s ed   t h tr an s is to r   C 5 4 6 B   as  s w itch .   W h e n   t h e   s q u ar w a v at   0   V,   th e   Evaluation Warning : The document was created with Spire.PDF for Python.
I J E C E     I SS N:  2088 - 8708       Mo d elin g   th V ib r a tio n a l D yn a mics o f P iezo elec tr ic  A ctu a to r   b y   S ystem  ....   ( N u r u l B a h ia h   Mo h d   N o o r )   1509   tr an s i s to r   wa s   i n   c u t - o f f   r e g io n   w h er at  t h is   p o in n o   cu r r e n at  co llecto r   ca u s e d   t h v o lt ag o f   P E A   d r o p p ed   n ea r l y   to   ze r o .   Ho w e v er ,   th h ig h   v o ltag let  t h tr a n s i s to r   to   s atu r ate  an d   r es u lti n g   i n   m ax i m u m   co llecto r   cu r r en a n d   m i n i m u m   v o lta g d r o p p ed .   T h er ef o r e,   th in p u t   v o lta g o f   P E A   at  th i s   m o m e n r ea c h ed   m ax i m u m   p ea k .   T h v ib r atio n   p h e n o m e n o n   v id e o   w a s   r ec o r d ed   b y   u s i n g   Di g i A cq u i s   s o f t w ar e.     T h en ,   t h i s   vi d eo   w a s   ad ap ted   in to   s till   i m ag e s   to   m ea s u r t h d is p lace m en o f   P E A .   T h s y s te m   id e n ti ca tio n   tec h n iq u wa s   e m p lo y ed   to   d ev elo p   m at h e m atica m o d els  o f   d y n a m ic   s y s te m   b ased   o n   s et   o f   m ea s u r ed   d ata  s a m p les.   T h e   d ev elo p ed   m o d el  w as   p r ed ictio n   m o d el  th a t   r ep r esen ted   th v ib r atio n   p h e n o m e n o n   o f   P E A .     T h s y s te m   id en ti f icatio n   ap p r o ac h   is   s h o w n   i n   Fi g u r 4   w h er th P E A   m o d el  w as  est i m ated   b ased   o n   in p u t - o u tp u d ata  co llectio n ,   f o llo w ed   b y   m o d el  s tr u ct u r an d   m o d el  v alid atio n .           Fig u r 4 .   S y s te m   I d en ti f icat io n   ap p r o ac h       T h m o d el  s tr u ct u r w a s   s elec ted   b y   tr y in g   d i f f er e n co m b i n atio n   o f   ze r o s   an d   p o les  u n t il  t h m o d el   g iv e   o p ti m u m   f it   f o r   te s d ata   an d   lo w er   m ea n   s q u ar ed   er r o r   ( MSE )   an d   f it n es s   p r ed ictio n   er r o r   ( FP E ) .   T h ex p er i m e n tal  d ata  co llectio n   w er i m p o r ted   in to   w o r k s p ac v ar iab le,   w h er e   in p u d ata  w er id en ti fied   as   v o ltag a n d   d is p lace m en d ata  w er e   id en ti ed   as  o u tp u t.   T h to o l b o x   id en ti ca tio n   i n   Ma tlab   R 2 0 1 3 a   a s h o w n   in   Fi g u r 2 . 3 ( a)   wa s   u s ed   f o r   m o d ellin g   p u r p o s e.   In   t h e   id en ti ficatio n   w in d o w ,   p r es s   th e   p o p u p   m e n u - i m p o r t   d ata   an d   s elec t   tim d o m ai n   d ata  to   im p o r th d ata.   I n   th i s   r esear ch ,   th er ar 3 0 0   d ata  o f   v o ltag an d   d is p lace m e n t.  T h er ef o r e,   f ir s 1 5 0   d ata  w er u s ed   f o r   es ti m a tio n   an d   t h r e m ain in g   1 5 0   f o r   v alid atio n .   Ne x t,   tr an s f er   f u n ctio n   m o d el  w as   o b tain ed   as  s h o w n   in   Fig u r 2 . 3 ( b ) .   I n   th is   p ap er ,   m o d el  v alid atio n   w a s   ev alu a ted   u s i n g   co r r elatio n   test   w h ic h   w er au to co r r elati o n   an d   cr o s s   co r r elatio n ,   MSE   an d   FP E .   T h co n f id e n ce   in ter v al  le v el  u s ed   f o r   r esid u al  test   w er 9 9 . 9 % f o r   all  m o d els.         ( a)                                                                                   ( b )     Fig u r 5 .   ( a)   Sy s te m   id en ti ca tio n   to o lb o x   w i n d o w   ( b )   T r an s f er   f u n c tio n   est i m a tio n   i n   s y s t e m   id en ti ficatio n   to o lb o x   Evaluation Warning : The document was created with Spire.PDF for Python.
                      I SS N :   2 0 8 8 - 8708   I J E C E    Vo l.  7 ,   No .   3 J u n e   2 0 1 7   :   1 5 0 6     1 5 1 2   1510   3.   RE SU L T A ND  AN AL Y SI S   A t   th is   p o in t,  th m at h e m atic al  m o d el  o f   P E A   a r p r ese n t ed .   T h r ee   tr an s f er   f u n ctio n s   h av e   b ee n   d ev elo p ed   b ased   o n   d ata  co lle ctio n   f o r   v ib r atio n   p h e n o m e n o n .   T h f r eq u en cie s   u s ed   f o r   v i b r atio n   ex p er i m e n t   we r 1   Hz,   2 0   Hz  an d   5 0   H z.   T h p u r p o s o f   v i b r atio n   ex p er i m e n tal  w a s   to   in v esti g at th d is p lace m e n t   b eh av io r   o f   P E A   at  v ar io u s   f r e q u en cie s .     i)   Fre q u en c y : 1   Hz   T h m o d el  s tr u ctu r o b tai n ed   at  f r eq u e n c y   o f   1   Hz  w a s   t h ir d   o r d er   s y s te m   w it h   th r ee   p o l es  an d   t w o   ze r o s .   T h tr an s f er   f u n ctio n   is   s h o w n   i n   eq u atio n   ( 1 ) .   S y s te m   id e n ti ca t io n   r esu lts   p r esen ted   th at  th e   estab lis h ed   m o d el  f it   v a lid atio n   d ata  b y   9 6 . 5 8   f i tn e s s   lev e w it h   m ea n   s q u ar ed   er r o r   ( MSE )   o f   0 . 0 4 9 4 8   an d   f i n al  p r ed ictio n   er r o r   ( FP E )   o f   0 . 0 4 2 1 .   Hig h   p er ce n tag o f   m o d el   t n es s   an d   lo w er   er r o r   lev el   r ec o g n ized   t h e   v alid it y   o f   m o d el.     T h e n,   r esid u al  an al y s is   w as  co n t in u ed   to   v alid ate  th m o d el.   T h m aj o r it y   o f   s a m p led   d ata  w er e   w i th i n   t h 9 9 . 9   co n d en ce   i n ter v al   li m it  f o r   au to co r r elatio n   f u n ctio n ,   s i m ilar l y ,   all   s a m p led   d ata  w er e   w it h i n   9 9 . 9   %   co n d e n ce   i n t er v al  li m it  f o r   cr o s s   co r r elati o n   f u n ctio n .   T h er ef o r e,   t h m o d el  h as   p as s ed   th e   v alid atio n   test .   T h s tep   r esp o n s at   1   Hz  i s   s h o w n   in   Fi g u r 6   an d   ti m r esp o n s cr iter ia  ca n   b o b s er v ed   i n   T ab le  1 .     G ( S ) = 8 . 633s 2 5427s 3523 s 3 + 80 . 36s 2 + 202 . 5s + 7651      ( 1 )           T ab le  1 .   T im r esp o n s cr iter i f o r   v ib r atio n   at  1   Hz               Fig u r 6 .   Step   R esp o n s f o r   v i b r atio n   at  1   Hz       ii)   Fre q u en c y : 2 0   Hz   T h tr an s f er   f u n c tio n   f o r   P E m o d el  at  2 0   Hz  w a s   g i v e n   in   eq u atio n   ( 2 ) .   T h b est  p r e d ictio n   m o d el  w a s   t h ir d   o r d er   w it h   t h r ee   p o les   an d   t w o   ze r o s T h s y s t e m   id e n ti ca tio n   s tated   th a t h m o d el  f it s   ti m e   d o m ai n   v a lid atio n   d ata  w it h   7 3 . 7   %,   MSE   o f     0 . 5 2 7 1   an d   FP E   o f   0 . 5 7 2 2 .   Me an w h ile,   t h r esid u al  a n al y s i s   s h o w ed   th a m o s o f   s a m p led   d ata  w er e   w it h i n   t h e   9 9 . 9   co n f id e n ce   b o u n d ar y   f o r   au to c o r r elatio n   an d   cr o s s   co r r elatio n .   T h er ef o r e,   th m o d el  h as   p ass ed   t h v alid atio n   test   f r o m   r esid u a an al y s is   p e r s p ec tiv e.   T h e   s tep   r esp o n s e   f o r   th i s   m o d el  is   s h o w n   i n   Fi g u r an d   th ti m r e s p o n s p er f o r m a n ce   is   p r ese n t ed   in   T ab le  2 .     G ( S ) = 0 . 732s 2 2 6 3 9s + 738 . 3 s 3 + 128 . 2 s 2 + 462 . 3s + 6 3 90         ( 2 )       T ab le  2 .   T im r esp o n s cr iter i f o r   v ib r atio n   at  2 0 Hz       R i se   T i me   ( s)   S e t t l i n g   T i me   ( s)   O v e r sh o o t ( %OS )   0 . 0 1 0 3   6 . 9 5 7 1   1 3 3 5   R i se   T i me   ( s)   S e t t l i n g   T i me   ( s)   O v e r sh o o t ( %OS )   0 . 0 0 8 9   3 . 5 4 0 8   9 0 1 . 5 3 8 9   Evaluation Warning : The document was created with Spire.PDF for Python.
I J E C E     I SS N:  2088 - 8708       Mo d elin g   th V ib r a tio n a l D yn a mics o f P iezo elec tr ic  A ctu a to r   b y   S ystem  ....   ( N u r u l B a h ia h   Mo h d   N o o r )   1511       Fig u r 7 .   Step   R esp o n s f o r   v i b r atio n   at  2 0 Hz       iii)   Fre q u en c y : 5 0   Hz   T h d ata  co llectio n   f r o m   e x p er i m e n t al   w er e   u s ed   to   p r ed ict  th P E A   m o d el  u s i n g   s y s te m   id en ti ca tio n   ap p r o ac h .   Hen c e,   th ir d   o r d er   tr an s f er   f u n ctio n   m o d el  w it h   t h r ee   p o les  an d   t w o   ze r o s   h a s   b ee n   estab lis h ed   f o r   v ib r atio n   e f f ec t   at  5 0   Hz   as  s h o w n   i n   ( 3 ) T h e   s tep   r esp o n s e   f o r   th i s   m o d el   is   s h o w n   i n   Fi g u r e   an d   th ti m r esp o n s cr it er ia  is   p r esen ted   in   T ab le  3 .   T h v alid atio n   o f   tr an s f er   f u n ct io n   as  g i v en   b y   eq u atio n   ( 3 )   f o r   v ib r atio n   wa s   o p ti m ized   b y   6 6 . 3 3   o f   to   v al id atio n   d ata,   M SE  o f   0 . 4 0 4 7   an d   FP E   o f   0 . 3 7 7 4 .   T h is   r esu lt  r ec o g n ize d   th v alid it y   o f   th m o d el  w i th   h i g h   p er ce n tag m o d el  t n ess   an d   lo w er   er r o r   lev el.   I n   ad d itio n ,   m aj o r it y   o f   s a m p led   d ata  w er e   w it h i n   th e   co n d en ce   i n ter v a li m it  f o r   au to co r r elatio n   an d   cr o s s   c o r r elatio n .   T h er ef o r e,   t h m o d el  h a s   p ass ed   th v alid atio n   test   a n d   ca n   b u s ed   f o r   th is   p r o j ec t.     G ( S ) = 1 . 91 s 2 1 . 5 6 2 e4 s 1 . 1 9e5 s 3 + 92 6 . 9s 2 + 4 . 3 2 6 e4 s + 3 . 2 1 8 e6                                                  ( 3 )                       T ab le  3 .   T im r esp o n s cr iter i f o r   v ib r atio n   at  5 0 Hz                       Fig u r 8 .   Step   R esp o n s f o r   v i b r atio n   at  5 0 Hz       T h s i m i lar   ex p er i m en ta s et u p   h as  b ee n   u s ed   to   es tab lis h ed   th P E A   m o d el.   T h er ef o r e,   th r ee   P E A s   m o d el   h a v e   b ee n   d ev elo p ed   at   v ar io u s   f r eq u en cie s   o f   1   Hz,   2 0   Hz  an d   5 0   Hz.     T h d is p lac e m en t b eh a v io u r   o f   P E A   ca n   b o b s er v ed   b ased   o n   ex p er i m e n tal  d ata  co lle ctio n .   A 1   Hz,   th d i s p lace m en o f   P E A   w a s   in cr ea s i n g l y   m o r t h an   d is p la ce m en at  5 0   Hz.   T h f r eq u en c y   i s   i n v er s el y   p r o p o r tio n al  to   d is p lace m e n o f   R i se   T i me   ( s)   S e t t l i n g   T i me   ( s)   O v e r sh o o t ( %OS )   0 . 0 0 2 4 5   1 . 2 0 1 6   4 4 8 . 2 4 5 6   Evaluation Warning : The document was created with Spire.PDF for Python.
                      I SS N :   2 0 8 8 - 8708   I J E C E    Vo l.  7 ,   No .   3 J u n e   2 0 1 7   :   1 5 0 6     1 5 1 2   1512   P E A .   T h u s ,   th e   m o d el  th at   h a s   b ee n   d ev elo p ed   b y   SI   tech n i q u s h o w ed   t h at  o v er s h o o w a s   h ig h   a 1   Hz   an d   d ec r ea s ed   as  th f r eq u e n c y   u p   to   5 0   Hz.   T h o v er s h o o o f   t h s y s te m   d r o p p ed   to   6 6 . 4 2 f r o m   i n itial   f r eq u en c y   a n d   r ea ch   f a s t r esp o n s as t h f r eq u en c y   at  5 0   Hz.       4.   CO NCLU SI O N   T h is   p ap er   d is cu s s ed   ab o u m o d eli n g   o f   p iezo elec tr ic  ac tu ato r   f o r   v ib r atio n   p h en o m en o n   b y   i m p le m en t in g   s y s te m   id en ti f i ca tio n   tec h n iq u e.   S y s te m   id e n tif icat io n   ap p r o ac h   h as   b ee n   u s ed   to   est i m a te  t h e   d y n a m ic  m o d el  o f   P E A   b y   i n p u t - o u tp u e x p er i m e n tal  d ata.   T h m o d el  o f   P E A   m u s b estab lis h ed   f ir s i n   o r d er   to   d esig n   th co n tr o ller   to   m i n i m ize  t h r esid u al   v ib r a tio n   i n   th e   s y s te m .   T h m o d el s   o f   P E th at   h a v e   b ee n   d ev elo p ed   b y   e m p lo y i n g   th is   tec h n iq u e   we r s elec ted   b ased   o n   h ig h   f it  to   v alid atio n   d ata,   lo w e s   MSE   an d   FP E ,   au to co r r elatio n   an d   cr o s s   co r r elatio n   f o r   v alid atio n   test .   T h en ,   t h co n tr o ller   ca n   b d esig n   b ased   o n   m o d el  r eq u ir e m e n t.       RE F E R E NC E S   [1 ]   B.   Ja ff e ,   P iez o e lec tri c   c e ra m ic s ,”   El se v ier,  v o l.   3 2 0 1 2 .   [2 ]   Q.  Xu ,   Id e n ti f ica ti o n   a n d   c o m p e n sa ti o n   o f   p iez o e lec tri c   h y ste r e sis  w it h o u m o d e li n g   h y ste re sis  in v e rse ,   IEE T ra n sa c ti o n o n   In d u stri a El e c tr o n ics ,   v o l .   6 0 ,   p p .   3 9 2 7 - 3 9 3 7 ,   2 0 1 3 .   [3 ]   B.   M .   Ba d a n d   W .   G .   A li ,   Na n o   p o siti o n in g   f u z z y   c o n tr o f o r   p iez o e lec tri c   a c tu a to rs,   In t.   J .   E n g .   T e c h . ( IJ ET - IJ ENS ),   v o l.   1 0 ,   p p .   7 0 - 7 4 ,   2 0 1 0 .   [4 ]   M .   A b d e laz iz,  e a l. A   H y b rid   F u z z y   Ro b u st  c o n tr o f o P iez o e lec tri c   a c tu a to rs,   J o u rn a o Ame ric a n   S c ien c e v o l.   8 ,   p p .   4 8 3 - 4 8 9 ,   2 0 1 2 .   [5 ]   S .   X ia o   a n d   Y.  L i,   M o d e li n g   a n d   h ig h   d y n a m ic  c o m p e n sa ti n g   th e   ra te - d e p e n d e n h y ste re sis   o f   p iez o e lec tri c   a c tu a to rs  v ia  a   n o v e m o d if ied   in v e rse   P re isa c h   m o d e l,   IEE T ra n sa c ti o n s   o n   Co n tro S y ste ms   T e c h n o l o g y ,   v o l.   2 1 ,   p p .   1 5 4 9 - 1 5 5 7 ,   2 0 1 3 .   [6 ]   Z.   Ch a n d   Q.  Xu ,   Re c e n a d v a n c e in   th e   c o n t r o o f   p iez o e lec tri c   a c tu a to rs,   In ter n a ti o n a J o u rn a o Ad v a n c e d   Ro b o ti c   sy ste ms ,   v o l.   1 1 ,   p p .   1 - 1 1 ,   2 0 1 4 .   [7 ]   Y.  L iu ,   e a l . ,   Hy ste re sis a n d   c re e p   m o d e li n g   a n d   c o m p e n sa ti o n   f o a   p iez o e lec tri c   a c tu a to u sin g   a   f ra c ti o n a l - o r d e M a x w e l re sistiv e   c a p a c it o a p p r o a c h ,   S ma rt M a ter i a ls  a n d   S tru c t u re s ,   v o l.   2 2 ,   p p .   1 1 5 0 2 0 ,   2 0 1 3 .   [8 ]   H.  S h e n ,   e t   a l. A   lo w - p o w e c ircu it   f o p iez o e lec tri c   v ib ra ti o n   c o n tro b y   sy n c h ro n ize d   sw it c h in g   o n   v o lt a g e   so u rc e s,   S e n so rs   a n d   a c t u a t o rs   A:  Ph y sic a l ,   v o l.   1 6 1 ,   p p .   2 4 5 - 2 5 5 ,   2 0 1 0 .   [9 ]   P .   L i,   e a l. A   sim p le  f u z z y   s y ste m   f o m o d e ll in g   o f   b o t h   ra t e - in d e p e n d e n a n d   ra te - d e p e n d e n h y ste re sis  in   p iez o e lec tri c   a c tu a to rs,   M e c h a n i c a S y ste ms   a n d   S ig n a l   Pro c e ss in g ,   v o l.   3 6 ,   p p .   1 8 2 - 1 9 2 ,   2 0 1 3 .   [1 0 ]   L .   L iu ,   e a l .,  De v e lo p m e n o f   a n   a p p ro a c h   to w a rd   c o m p re h e n siv e   id e n ti f ica ti o n   o f   h y ste re t ic  d y n a m ics   in   p iez o e lec tri c   a c tu a to rs,   IEE T r a n sa c ti o n o n   Co n tro S y ste ms   T e c h n o l o g y ,   v o l.   2 1 ,   p p .   1 8 3 4 - 1 8 4 5 ,   2 0 1 3 .   [1 1 ]   J.  R F o rtu n ,   e a l. F latn e ss - b a se d   a c ti v e   v ib ra ti o n   c o n t ro f o p i e z o e lec tri c   a c tu a to rs,   IEE E/ AS M T ra n sa c ti o n s   o n   M e c h a tr o n ics ,   v o l .   1 8 ,   p p .   2 2 1 - 2 2 9 ,   2 0 1 3 .   [1 2 ]   M.  Jo u a n e h ,   F u n d a m e n tals o f   m e c h a tro n ics ,”   Ce n g a g e   Lea rn in g ,   2 0 1 2 .   [1 3 ]   P .   Ro n k a n e n ,   e a l. Disp lac e m e n c o n tro o f   p iez o e lec tri c   a c tu a to rs  u si n g   c u rre n a n d   v o lt a g e ,   IEE E/ AS M T ra n sa c ti o n o n   M e c h a tro n ics ,   v o l.   1 6 ,   p p .   1 6 0 - 1 6 6 ,   2 0 1 1 .   [1 4 ]   P .   S e n g u p ta  a n d   B.   L i,   M o d if i e d   Bo u c W e n   m o d e f o h y ste r e sis  b e h a v io o f   RC  b e a m c o lu m n   jo in ts  w it h   li m it e d   tran sv e rse   re in f o rc e m e n t,   En g i n e e rin g   S tru c t u re s ,   v o l.   4 6 ,   p p .   3 9 2 - 4 0 6 ,   2 0 1 3 .   [1 5 ]   V .   Ha ss a n i,   e a l. A   s u rv e y   o n   h y st e re sis  m o d e li n g ,   id e n ti f ica ti o n   a n d   c o n tr o l ,”   M e c h a n ica sy st e ms   a n d   sig n a l   p ro c e ss in g ,   v o l.   4 9 ,   p p .   2 0 9 - 2 3 3 ,   2 0 1 4 .   [1 6 ]   J.   W .   L ian g ,   e a l. M o d e l - Ba se d   Co n tr o f o P iez o e lec tri c - A c t u a ted   S y ste m Us in g   In v e rse   P ra n d tl - Ish li n sk ii   M o d e a n d   P a rti c l e   S w a r m   Op ti m iza ti o n ,   In ter n a ti o n a J o u rn a o M e c h a n ic a l,   Aer o sp a c e ,   In d u stria a n d   M e c h a tro n ics   E n g i n e e rin g ,   v o l.   8 ,   p p .   6 0 9 ,   2 0 1 4 .   [1 7 ]   H.  Ho ss e in i,   e a l. T h e   d e sig n   o f   ro b u st  s o f se n so u sin g   A NFIS   n e t w o rk ,   J o u rn a o I n stru me n ta ti o n   T e c h n o l o g y ,   v o l.   2 ,   p p .   9 - 1 6 ,   2 0 1 4   [1 8 ]   G .   W a n g ,   e t   a l. M o d e li n g   a n d   i d e n ti f ica ti o n   o f   a s y m m e tri c   Bo u c W e n   h y ste re sis  f o p iez o e lec tri c   a c tu a to v ia  a   n o v e d if f e re n ti a e v o lu ti o n   a lg o ri th m ,   S e n so rs   a n d   Act u a t o rs   A:  P h y sic a l ,   v o l.   2 3 5 ,   p p .   1 0 5 - 1 1 8 ,   2 0 1 5 .   [1 9 ]   M .   N.  M a sla n ,   e a l. De v e lo p m e n o S y ste m   Id e n ti f ica ti o n   f o P iez o e lec tri c   P a tch   A c tu a t o r,   in   Ap p li e d   M e c h a n ics   a n d   M a ter ia ls ,   p p .   2 4 5 - 2 4 9 2 0 1 5 .   Evaluation Warning : The document was created with Spire.PDF for Python.