Ghostscript wrapper for D:\Digitalizacja\MTS84_t22z1_4_PDF_artyku³y\mts84_t22z3_4.pdf M E C H A N I K A  T E O R E T Y C Z N A  I  S T O S O W A N A  3 ­ 4  22  (1984)  F A L E  S P R Ę Ż O NE  M E C H A N O ­ T E R M O ­ E L E K T R Y C Z N E  W  S T A Ł Y M  O Ś R O D KU  P I E Z O E L E K T R Y C Z N Y M  K R Y S T Y N A  M A J O R K O W S K A ­ K N A P  ( W A R S Z A W A )  Politechnika  Warszawska  1.  Wprowadzenie  Ostatnie  lata  ś wiadczą  o  wciąż  wzrastają cym  zainteresowaniu  zagadnieniami  pól  sprzę ż onych.  Badania  z a r ó w n o  teoretyczne,  jak  i  d o ś w i a d c z a l ne  zachowania  się  sprę ż y­ stego  ciała  stałego,  k t ó r e  znajduje  się we wzajemnym  o d d z i a ł y w a n i u  z polem  elektromagne­ tycznym,  a  czę sto  r ó w n i e ż  z  polem  temperatury,  stały  się  bardzo  popularne  z  powodu  olbrzymiego  zapotrzebowania  ze  strony  techniki.  M o ż na  powiedzieć,  że  piezoelektrycz­ n o ś ć  — dziedzina  zajmują ca  się  sprzę ż eniem  pola  deformacji  z  polem  elektromagnetycz­ nym  oraz  t c r m o ­ p i e z o e l e k t r y c z n o ś ć,  włą czają ca  ponadto  pole  temperatury,  są  na  etapie  dynamicznego  rozwoju.  Cztery  lata  temu  m i n ę ła  setna  rocznica  odkrycia  zjawiska  piezoelektrycznego  przez  braci  Piotra  i  P a w ł a  Curie.  W  1880 r.  odkryli  oni  prosty  efekt  piezoelektryczny  w  kry­ stalicznym  kontinuum  [1] — wytwarzanie  elektrycznej  polaryzacji,  a  więc  i  pola  elek­ trycznego  przez  mechaniczne  n a p r ę ż e n i e.  Odwrotny  efekt  piezoelektryczny — pojawianie  się  n a p r ę ż eń  i  o d k s z t a ł c e ń  k r y s z t a ł u  wskutek  zastosowania  z e w n ę t r z n e go  pola  elektrycz­ nego  był przewidziany  teoretycznie  przez  G . L i p p m a n a  [2] w  1881 r.  jako  konsekwencja  termodynamiczna  prostego  efektu.  W  tym  samym  r o k u  bracia  Curie  zweryfikowali  eks­ perymentalnie  to  odkrycie  [3].  Obydwa  efekty  zostały  udowodnione  i  doś wiadczalnie  zademonstrowane  w  anizotropowych  niecentrosymetrycznych  k r y s z t a ł a c h ,  np.  [4, 5]  w  pewnych  syntetycznych  m a t e r i a ł a c h ,  np.  ceramice  piezoelektrycznej,  np.  [6, 7], w tzw.  piezoelektrycznych  teksturach  (np.  drewno  [89], s k ó r a  [9],  k o ś ć  [10), w  ceramice  ferro­ elektrycznej  [11], w  p ł y n n y c h  k r y s z t a ł a c h  [12]. O b s z e r n ą  bibliografię  dotyczą cą  znanych,  piezoelektrycznych  m a t e r i a ł ó w  m o ż na  znaleźć  w  [13]. Badania  eksperymentalne  wykazują   że  efekt  piezoelektryczny  wystę puje  r ó w n i e ż  w  centrosymetrycznych  k r y s z t a ł a c h ,  np.  [14,  15].  Najbardziej  ś cisłe  s f o r m u ł o w a n i e  zależ noś ci  mię dzy  piezoelektrycznoś cią  i  s t r u k t u r ą   k r y s z t a ł ó w  z o s t a ł o  podane  przez  Woldemara  Voigta.  O n  właś nie  ustalił,  w  k t ó r y c h  z  32  klas  krystalograficznych  m o ż e  istnieć  efekt  piezoelektryczny,  i  p o d a ł  dla  każ dej  klasy,  k t ó r e  z osiemnastu  w s p ó ł c z y n n i k ó w  piezoelektrycznych  mogą  mieć  w a r t o ś ć  róż ną  od zera.  W  roku  1894 z o s t a ł a  opracowana  przez  W .  V o i g t a  fenomenologiczna  l i n i o w a  teoria  piezoelektrycznoś ci  [16]  oparta  na  zasadach  termodynamicznych  s f o r m u ł o w a n y c h  przez  K e l v i n a  5 1 0  К .  M A J O R K O W S K A ­ K N A P  Przez  ponad  35  lat  od  o d k r y c i a  piezoefektu  nie  p r z y w i ą z y w a no  ż adnej  wagi  do  jego  praktycznego  znaczenia.  Praktyczne  zastosowania  z a p o c z ą t k o w a ła  praca  P.  Langevina,  k t ó r y  w  1917  roku  uż ył  przetwornik  kwarcowy  do  wzbudzenia  fal  akustycznych  w  wodzie.  N a s t ę p n ie  p o w s t a ł y  pierwsze  generatory,  rezonatory  i  filtry  piezoelektryczne  dzię ki  pra­ com  P.  Langevina,  W .  G .  Cady,  A .  M .  N i c o l s o n a  i  K .  S.  V a n  D y k e a .  W  r o k u  1919  W .  G .  Cady  pierwszy  z a s t o s o w a ł  kryształy  piezoelektryczne  do  kontroli  czę stotliwoś ci  genera­ tora  d r g a ń .  N a s t ę p n ie  stosowali  je  G .  W .  Pierce  i  M i l l e r .  O k o ł o  dziesię ciu  m i l i o n ó w  ge­ n e r a t o r ó w  kontrolowanych  dzię ki  piezoefektowi  w  k r y s z t a ł a c h  było  u ż y w a n y ch  podczas  II  wojny  ś wiatowej  w  łą cznoś ci  mię dzy  czołgami,  ziemią  i  samolotami.  W  tych  zastosowa­ niach  wykorzystywano  element  piezoelektryczny  drgają cy  w  zmiennym  polu  elektrycznym.  Wyjaś nienie  własnoś ci  piezoelektrycznych  kryształów  umoż liwiły  badania  ich  struktury  za  p o m o c ą  promieni  Rontgena.  Dalszy  dynamiczny  rozwój  z a s t o s o w a ń  praktycznych  (takich  j a k :  w  p r z y r z ą d a ch  re­ jestrują cych  dź wię k,  w  p r z y r z ą d a ch  pomiarowych  ciś nień  szczególnie  g w a ł t o w n y c h  eks­ plozji,  a  również  ciś nienia  k r w i ,  n a p r ę ż eń  w  inż ynierskich  elementach  konstrukcyjnych,  do  retransmisji  i  mierzenia  p r ę d k o ś ci  ś wiatła,  w  telewizji,  komputerach,  do  wykonywania  w z o r c ó w  czasu,  w  r ó ż n y ch  typach  broni  ogniowej  i  artyleryjskiej)  p o w o d o w a ł  intensy­ fikację  z a r ó w n o  b a d a ń  podstawowych,  jak  i  aplikacyjnych.  W  o m ó w i o n y c h  zastosowaniach  wykorzystywano  obję toś ciowe  fale  piezoelektryczne.  Olbrzymie  zastosowanie  techniczne  piezoelektrycznych  fal  powierzchniowych  spowo­ d o w a ł o ,  że  skupiają  one  wciąż  na  sobie  wiele'uwagi  badaczy.  Ze  wzglę du  na  ich  bardzo  niską  p r ę d k o ść  oraz  bardzo  mał ą  d ł u g o ś ć  w  p o r ó w n a n i u  z  falami  elektromagnetycznymi  tej  samej  czę stotliwoś ci  (redukcja  w  wymiarze  z a r ó w n o  prę dkoś ci  jak  i  długoś ci  o k o ł o  1 0 " s )  znajdują  zastosowanie  m.in.  w  nieniszczą cych  oszacowaniach  sejsmologicznych,  w  o b r ó b c e  s y g n a ł ó w  z  w a ż n y mi  zastosowaniami  m.in.  do  radaru,  telekomunikacji,  elek­ tronowej  techniki  wojskowej.  Akustoelektronika  p o w s t a ł a  w  cią gu  ostatnich  20  lat  j a k o  osobna  gałąź  nauki  i  techniki  bazuje  na  piezoelektrycznych  falach  powierzchniowych.  Odkryto,  zbadano  i  zastosowano  szereg  typaw  fal  powierzchniowych  od  dawno  znanych  fal  Rayleigha  (17],  L o v e ' a  [17],  Stoneleya  [18]  do  stosunkowo  niedawno  odkrytych,  np.  Bleusteina­Gulajewa  [19,  20],  typu  „ I e a k y "  [21,  22],  wystę pują cych  z a r ó w n o  na  p ł a ­ skich,  jak  i  zakrzywionych  (sferycznych  bą dź  cylindrycznych)  powierzchniach.  W  niniejszej  pracy  skoncentrowano  się  na  wybranych  zagadnieniach  dynamicznych,  a  mianowicie  interesują cych  typach  fal  (o  przebiegach  harmonicznych  w  czasie  i  w  prze­ strzeni)  w  o ś r o d ku  piezoelektrycznym.  2.  Fale  harmoniczne  w  oś rodku  piezoelektrycznym —  teoria  W .  Voigta  Klasyczne  s f o r m u ł o w a n i e  teorii  piezoelektrycznoś ci  bazują ce  na  podstawowych  r ó w ­ naniach  liniowych  wyprowadzonych  przez  W .  Voigta  zawiera  fundamentalna  literatura:  prace  z  lat  50­tych  W .  G .  Cady  [4]  (pierwsze  wydanie  1946  г .),  W .  P.  M a s o n a  [5],  z  lat  б О ­г ­70­tych  W .  P.  M a s o n a  [23],  H .  F .  Tierstena  [24,  25]  oraz  wydana  w  1983  r.  mono­ grafia  W .  Nowackiego  [26].  L i n i o w a  teoria  piezoelektrycznoś ci  W .  Voigta  w  quasi­statycznym  przybliż eniu  [26]  M E C H A N O ­ T E R M O ­ E L E K T R Y C Z N E  P A L B  511  reprezentowana  jest  przez  nastę pują ce  podstawowe  zwią zki  i  r ó w n a n i a  napisane  w  kar­ tezjań skim  u k ł a d z i e  w s p ó ł r z ę d n y ch  p r o s t o k ą t n y c h:  —  r ó w n a n i a  ruchu  i  r ó w n a n i e  pola  elektrycznego  Ojt.j+Xi  = 4 « i »  £ i . i = 0»  xeB,  / > 0  (2.1)  —  zwią zki  konstytutywne  °Ч/  =  cljklEkl~  ekijEk  Di  = eMeM  + eikEk,  Ek = ­<рл   —  definicję  tensora  odkształcenia  1  (2.2)  e tj  *=­b­0  (3.1)  oraz  r ó w n a ń  M a x w e l l a  pola  elektromagnetycznego  [65]  r o t H  =  D  +  J ,  r o t E  =  ­B,  (3.2)  d i v D  =  0,  d i v B  =  0.  (3.3)  gdzie  H ,  E ,  D ,  В ,  I  —  oznaczają  odpowiednio  wektory:  n a t ę ż e n ia  pola  magnetycznego,  n a t ę ż e n ia  pola  elektrycznego,  indukcji  elektrycznej,  indukcji  magnetycznej  i n a t ę ż e n ia  p r ą du  przewodzenia  (pomija  się  prąd  konwekcyjny).  Wektory  charakteryzują ce  pole  elektromagnetyczne  są  p o w i ą z a ne  r ó w n a n i a m i  k o n ­ stytutywnymi  D = e 0 E + P ,  В  =  fi0H,  J  =  crE  (3.4)  Tutaj  P  jest  wektorem  polaryzacji  elektrycznej,  e0  oznacza  p r z e n i k a l n o ś ć  elektryczną   M E C H A N O ­ T E R M O ­ E L E K T R Y C Z N E  F A L E  5 1 5  p r ó ż n i,  ц0  p r z e n i k a l n o ś ć  magnetyczną  p r ó ż n i,  a  jest  przewodnoś cią  elektryczną  właś ciwą.  (Założ yliś my  tu  p r o p o r c j o n a l n o ś ć  w e k t o r ó w  I  i  E  oraz  i z o t r o p o w o ś ć  materiału ,  jeś li  chodzi  o  p r z e w o d n o ś ć  elektryczną ).  W  w y n i k u  wykonania  operacji  rotacji  na  r ó w n a n i u  (3.2)2  przy  wykorzystaniu  (3.2),  i  zwią zku  (3.4)2  dochodzimy  do  r ó w n a n i a :  r o t r o t E  =  fx0b­LioaE  (3.5)  M a m y  więc  do  rozwią zania  u k ł a d  r ó w n a ń  (3.1)  i  (3.5).  Po  wstawieniu  do  niego  zwią zków  konstytutywnych  dla  zagadnienia  quasi­statycznego  (2.2)  przy  wykorzystaniu  definicji  tensora  o d k s z t a ł c e n i a  (2.3)  otrzymamy  u k ł a d  sześ ciu  r ó w n a ń ,  w  w y n i k u  rozwią zania  k t ó r e g o  wyznaczamy  s k ł a d o w e  wektora  przemieszczenia  щ  i  s k ł a d o w e  natę ż enia  pola  elektrycznego  Et,  co  pozwoli  na  obliczenie  p o z o s t a ł y c h  wielkoś ci  charakteryzują cych  sprzę ż one  pola.  N i e  podajemy  w a r u n k ó w  brzegowych  i  p o c z ą t k o w y c h,  gdyż  bę dziemy  zajmowa ć  się  w  dalszej  kolejnoś ci  jedynie  ciałem  nieograniczonym  z  warunkami  regular­ noś ci  w  n i e s k o ń c z o n o ś c i.  Jednowymiarowy  problem  propagacji  fal  w  nieograniczonym  o ś r o d ku  piezoelektrycz­ nym  był  analizowany  w  pracy  J.  J .  K a y m a  [64]  i  M .  Ahmeda  [66].  W a r t o  o d n o t o w a ć ,  że  w  ramach  tej  teorii  moż liwych  jest  pięć  fal  płaskich  dla  danego  kierunku  propagacji,  ale  dwie  z  nich  są  m a ł o  interesują ce  pod  ką tem  z a s t o s o w a ń ;  mają  bowiem  p r ę d k o ść  bliską   p r ę d k o ś ci  fal  elektromagnetycznych  w  dielektrykach.  P o z o s t a ł e  trzy  interesują ce  nas  są   właś nie  tymi,  k t ó r y c h  przybliż one  p r ę d k o ś ci  daje  uproszczenie  quasi­statyczne.  Badanie  sprzę ż eń  towarzyszą cych  propagacji  fal  sprę ż ystych  i  elektromagnetycznych  w  krysztale  tetragonalnym  klasy  4  mm  było  przedmiotem  b a d a ń  K .  Majorkowskiej­Knap  [67].  P o r ó w n a n i e  dwu  r ó ż n y ch  p r z y p a d k ó w  przy  uwzglę dnieniu  bą dź  p o m i n i ę c iu  p r ą du  przewodzenia  (jest jedynie  p r ą d  przesunię cia)  wskazuje  na  pojawienie  się  efektów  dyspersji  i  t ł u m i e n i a  fal  w  pierwszym  z  nich.  Propagacji  czystych  m o d ó w  fal  powierzchniowych  na  powierzchni  p ó ł n i e s k o ń c z o n e go  piezoelektryka  były  poś wię cone  m.in.  prace  С.  C .  Tsenga  [68,  69].  Zagadnienie  wzbudzania  fal  mechanicznych  w  niepiezoelektrycznym  o ś r o d ku  poprzez  wzbudzenie  fal  elektromagnetycznych  w  przylegają cej  warstwie  piezoelektrycznej  było  badane  teoretycznie  w  pracy  F .  Bardati,  G .  Barzilal,  G .  Gerosa  [70].  4.  Płaskie  fale  harmoniczne  w  nieograniczonym  oś rodku  termopiezoelcktrycznym  W a ż n ym  u z u p e ł n i e n i e m  fal  piezoelektrycznych  są  fale  termopiezoelektryczne.  W  ter­ mopiezoelektrycznoś ci  (podstawowe  r ó w n a n i a  R .  D .  M i n d l i n  [71],  rozwinię cie  teorii  W .  N o w a c k i  [72])  opieramy  się  na  termosprę ż ystoś ci  [73]  i  elektrodynamice  o ś r o d k ów  cią głych  [65,  74].  Ciał o  sprę ż yste,  na  k t ó r e  działają  siły  masowe  i  ź r ó d ła  ciepła,  obcią ż enia  z e w n ę t r z ne  i  pole  elektryczne  oraz  o d d z i a ł y w a n i a  termiczne,  doznaje  deformacji.  Powstaje  W  n i m  pole  przemieszczeń  u(x,  t)  i  sprzę ż one  z  nim  pole  temperatury,  które j  przyrost  0  ==  T  — To,  gdzie  T  >  0,  a  T0  jest  t e m p e r a t u r ą  stanu  naturalnego  w  skali  K e l v i n a .  W y ­ mienione  pola  powodują  z m i a n ę  pola  elektromagnetycznego.  W  quasi­statycznym  przy­ bliż eniu  liniowej  teorii  t e r m o p i e z o e l e k t r y c z n o ś ci  dysponujemy:  13»  516  К .  M A J O R K O W S K A ­ K N A P  —  r ó w n a n i a m i  ruchu  i  r ó w n a n i e m  pola  elektrycznego  tU+pt  в ,,  = 0.  R ó w n a n i a  (4.4)  i (4.5)  tworzą  komplet  s p r z ę ż o n y ch  ze  sobą  r ó w n a ń  r ó ż n i c z k o w y ch  termo­ piezoelektrycznoś ci.  Wystę pują  w nich  nieznane  wielkoś ci  ut,   Ejt.j + Et­fp.i  + E?  =  0  dlaxeB  (5.1)  cptii  = 0  dla х  G В '  gdzie  ( 5 . 1 ) ! — t o  r ó w n a n i a  ruchu,  (5.1)2  —  r ó w n a n i e  bilansu  siły  wewną trzczą steczko­ wej  —  fundamentalny  warunek  r ó w n o w a g i  p o w ł o k i  (zewnę trznych  e l e k t r o n ó w )  pod  d z i a ł a n i e m  rdzenia  ( n u k l e o n ó w  i w e w n ę t r z n y ch  e l e k t r o n ó w )  tego  samego  atomu  i  otacza­ j ą c e go  pola  elektromagnetycznego  oraz  dodatkowe  działanie  są siednich  a t o m ó w  na  po­ w ł o k ę .  (5.1)3  —  r ó w n a n i e  pola  elektromagnetycznego,  do k t ó r e g o  b e z p o ś r e d n io  prowadzi  wstawienie  do Dti  =  QC  zwią zków  D{  =  е 0 2 ? | + Р |  oraz  E% =  —,i,  —  r ó w n a n i a  konstytutywne  <*u =  c,2uk_kdu  + c^(ul,j  + Ujil)  +  al2Pktkdij  +  d^ fPJę i+Piij),  (5.2)  Eu  + ^  2 ukik  6U  + dU"t.j  + Щ . t) + bi2  Pt.к  bu  + bUPj.  < +  +*>77P.,,< ­  +  6 4 , ,  ЕГ  =  ­aP(.  5 1 8  К .  M A J O R K O W S K A ­ K N A P  Wstawienie  zwią zków  (5.2)  do r ó w n a ń  (5.1)  prowadzi do nastę pują cego  u k ł a d u  sprzę ż onych  r ó w n a ń  piezoelektrycznoś ci  C 4 4 V 2 u ­ r ( c 1 2 ­ r c 4 4 ) g r a d d i v u +  rf44V 2P  + (J 1 2 ­l­u?44)graddivP + X  = QU, dĄ ĄV 2n+(di  г  +  7 7 ) g r a d d i v P ­ f l P ­ g r a d 9 ? + E °  =  0,  —  £oV29 3 + d i v P  =  Qe,  dla  xeB,  V2c>  =  0,  dla  xe  B',  D l a  r ó w n a ń  d o c h o d z ą  warunki  brzegowe  aJnnJ  =  qt,  EJnnJ  =  0,  i­e0\