Ghostscript wrapper for D:\BBB-ARCH\ARCHIWUM-lata-78-71\MTS73_t11z1_4\MTS72_t10z1_4\mts72_t10z1.pdf M E C H AN I K A TEORETYCZNA I  STOSOWANA 1,  10  (1972) AN ALIZ A  R U C H U   P E WN E G O  U KŁAD U   WI BR O- U D E R Z E N I OWE GO O  D WÓ C H   ST O P N I AC H   SWOBOD Y Z B I G N I E W  W I Ś N I E W S KI  ( G D A Ń S K) Oznaczenia Xi  bezwymiarowe  przemieszczenie  masy  1 x 2   bezwymiarowe  przemieszczenie  masy  2 Xi  bezwymiarowa  prę dkość  masy  1 X2  bezwymiarowa  prę dkość  masy  2, Xi  bezwymiarowe  przyspieszenie  masy  1, x 2   bezwymiarowe  przyspieszenie  masy  2, t[s]  czas, T  czas  bezwymiarowy, n  stosun ek  okresu  ruchu  do  okresu  sił y  wymuszają cej, cuts"1]  czę stość  zmian  sił y  wymuszają cej, (p[rd]  ką t  przesunię cia  fazowego  pomię dzy sił ą  wymuszają cą   a przemieszczeniem  masy 1, R  współ czynnik  restytucji, v  bezwymiarowa  prę dkość  uderzenia, i5  bezwymiarowe  przemieszczenie ś rodka  masy  ukł adu w cią gu jednego  cyklu  ruchu, sztywność  sprę ż ystego  zawieszenia  masy  2, hi  bezwymiarowy  współ czynnik  oporu  tł umienia  w  pneumatycznym  elemencie sprę ż ystym, h 2   bezwymiarowy  współ czynnik  oporu  tł umienia  drgań  masy 2, Qi[N ]  cię ż ar  masy  1, Q 2   [N ]  cię ż ar  masy  2, P[N / m 2 ]  ciś nienie  w  pneumatycznym  elemencie  sprę ż ystym, S[m 2 ]  pole  powierzchni  dn a  tł oka  elementu  sprę ż ystego, F 0 [N ]  am plituda  sił y  wymuszają cej. 1.  Wprowadzenie D oś wiadczalne  i  teoretyczn e  badan ia  ukł adu  wibro- uderzeniowego  z  elementem  sprę - ż ystym,  umoż liwiają cym  cią głą   regulację   czę stoś ci  uderzeń  wykazał y  [4], że  ukł ad  ten  po- siada  szereg  cech,  odróż niają cych  go  od  innych  ukł adów  wibro- uderzeniowych. 54 Z.  WIŚ N IEWSKI Korzyś ci,  jakie  mogł oby  przynieść  zastosowanie  takiego  ukł adu  n p.  w  wibro- uderze- niowych  urzą dzeniach  do pogrą ż ania  pali,  był y  przyczyną   podję cia  przez  autora  dalszych badań  teoretycznych  [5],  dla gł ę bszego  poznania  zjawisk  towarzyszą cych  pracy  ukł adu. W  pracy  [5] badany był  model  dynamiczny  ukł adu o dwóch  stopniach swobody,  w któ- rym  ogranicznik  posiadał   skoń czoną   masę   i  liniową   sprę ż ystoś ć,  przy  czym  przyję to,  że ś rodek  masy  ukł adu porusza  się   ruchem jednostajnym. Rys.  1.  Schemat  ukł adu  badanego  w  pracy  [5] Zał oż ono  tam  również,  że sił a  oporu  towarzyszą ca  ruchowi  obu mas  drgają cych  jest proporcjonalna  do prę dkoś ci  ruchu. M odel  dynamiczny  tego  ukł adu przedstawiono  schematycznie n a rys.  1. Ruch ukł adu opisany  był   równaniami  róż niczkowymi (1.1)  x 1 +h i (x 1 - x 2 )+q=- cos (1.2) x2-   - ±-  ( Xl- x2)+ 2h2x2+ s 2x2  =  —, przy  czym  warunki  graniczne1)  miał y  postać Xj(0)  =  x2(0)  =  x0,  xx(2nń )  = x2(2mi)  =  xo— 5,  Xi(0)—x2(Q)==Rv, Równania  (1.1.) i  (1.2)  wyraż ono  we  współ rzę dnych bezwymiarowych,  przy  pomocy pod- stawień : (1.4) x  =   cot; 1J 2   = ln 2 , g  g Równania  (1.1) i  (1.2)  stanowią   ogólny  opis  ruchu  badanego  ukł adu. Z uwagi  jedn ak na zł oż oną   formę   rozwią zań,  wystę pują   trudnoś ci  n atury  matematycznej  w  przedstawieniu x)  Termin  «warunki  graniczne» stosuje  się  w teorii  ukł adów  wibro- uderzeniowych  dla  odróż nienia od typowej postaci warunków brzegowych, niezależ nych od parametrów  ukł adu. AN ALIZA  RU CH U   PEWNEGO  UKŁADU   WIBRO- UDERZENIOWEGO 55 wyników  analizy  w  postaci  dogodn ej  dla  zastosowań  praktycznych.  D latego  też  celowe jest  wprowadzenie  pewnych  zał oż eń  upraszczają cych. W  niniejszej  pracy  zajmiemy  się   badan iem  ukł adu wibro- uderzeniowego  o dwóch  stop- n iach  swobody,  zawierają cego  elem ent  sprę ż ysty  o  charakterystyce  niezależ nej  od  poł o- ż enia  m asy  drgają cej  [4],  w  przypadku  gdy  m oż na  pom iną ć  tł umienie w  elemencie  sprę - ż ystym, ja ko  zn ikom e w  porówn an iu  z  sił ą   oporu  ruchu masy  2  (rys. 2). Z godn ie  z  przyję tymi  wyż ej  oznaczeniami  mamy  wię c  li 1   g  0.  D la  takiego  przypadku okreś limy  okresowe  rozwią zan ia  równ ań  ru ch u ;  warun ki  istnienia  rozwią zań  okresowych, wynikają ce  z  m atem atyczn ych  i  fizycznych  ograniczeń  wartoś ci  param etrów  ukł adu, jak również  kryteria  stabilnoś ci  strukturaln ej. Oddzielnie  rozpatrzon y  zostanie  przypadek  ruchu,  gdy  n a  masę   2  dział a  sił a  oporu o  stał ej  wartoś ci  przył oż ona  skokowo  w  dowolnej  chwili,  w  przedziale  czasu  pomię dzy uderzeniam i. 2.  Sformułowanie problemu R ozważa  się   ukł ad wibro- uderzeniowy  z pneum atycznym elementem sprę ż ystym,  w któ- rym  wpł yw  zmiany poł oż en ia t ł o ka wzglę dem  cylin dra2)  n a ciś nienie w przestrzeni sprę ż ania jest  znikomy.  Ogran iczn ik  ruch u  masy  1  (rys.  2)  stanowi  masę   skupioną   2,  zawieszoną F o -  cos  cot Rys.  2.  Model  dynamiczny  ukł adu  wibro- uderzeniowego  z  pneumatycznym  elementem  sprę ż ystym, bez oporów  ruchu  w  elemencie  sprę ż ystym n a  sprę ż ynie  o  liniowej  charakterystyce.  R uchowi  ogranicznika  towarzyszy  sił a  oporu proporcjon aln a  do  prę dkoś ci. D o  analizy  przyjmiemy  nastę pują ce  zał oż enia. 1.  M oż liwy  jest  ruch  okresowy  masy  1,  przy  czym  okres  ruchu  równy  jest  okresowi sił y  wymuszają cej  lu b  jego  cał kowitej  krotn oś ci; 2.  sił a  wymuszają ca  zm ien ia  się   harm on iczn ie i  dział a w  kierunku  zgodnym  z  kierun- kiem  osi  symetrii  cylin dra  w  elemencie  sprę ż ystym,  zaś  czę stość  zmian  sił y  wymuszają cej jest  stał a w  cią gu  cyklu  ruch u u kł a d u ; 3.  nie wystę puje  sprzę ż enie  zwrotn e  pomię dzy  ukł adem a  ź ródł em energii; 2 )  Budowę  pneumatycznego  elementu sprę ż ystego  opisano szczegół owo w pracy [4]. 56  Z .  WIŚ N IEWSKI 4.  czas  trwania  zderzeń pomię dzy  masą  1 i  ogranicznikiem  2 jest  krótki  w  porówn an iu z  okresem  ruchu  ukł adu  (odstę pem  czasu  pomię dzy  dwom a  kolejnymi  zderzen iam i); 5.  przekazywanie  energii  kinetycznej  przy  zderzeniu  okreś la  współ czynnik  restytucji R,  przybierają cy  wartoś ci  z przedział u  [0;  1] oraz speł niają cy  zał oż enia podan e w pracy  [4]; 6. 'opory  ruchu  tł oka  wzglę dem  cylindra  w  elemencie  sprę ż ystym  są  nieznaczne  w  po- równaniu  z  oporam i  ruchu  masy  2; 7.  ruch ogranicznika  2 traktujemy  jako  drgan ia  masy  skupionej  zawieszonej  n a  liniowej sprę ż ynie; 8.  ś rodek  masy  cał ego ukł adu przesuwa  się ruchem jedn ostajn ym ,  w  kierun ku  ujemnego zwrotu  osi  x; 9.  poł oż enie ś rodka  masy  1 nie  zależy  od  chwilowego  poł oż en ia m as  niewyważ onych wibratora  bezwł adnoś ciowego,  stanowią cego  ź ródło sił y wymuszają cej  drgan ia; 10.  wszystkie  param etry  ukł adu  są  zdeterm in owan e. U wzglę dniając  powyż sze  zał oż enia,  ruch  ukł adu  opisać  m oż na  dwom a  równ an iam i róż niczkowymi: (2.1)  Jći+ ? =  cos( (2.2)  x 2 +2h 2 x 2 +s 2 x 2   =   1 z  warunkam i  okresowoś ci  ruchu  ,(1.3). 3.  Cał kowanie  równań  ruchu Ponieważ  w  równaniu  (2.1)  nie  wystę pują  czł ony  zależ ne  od  x 2)   zaś  równ an ie  (2.2) nie  zawiera  czł onów  zależ nych  od  x,  więc  zwią zek  pom ię dzy  param etram i  ruch u  obu mas  wynika  jedynie  z  warun ków  okresowoś ci  ruchu  (1.3) 3). Cał ka  równania  (2.1)  ma  postać (3.1)  x t   =   -   ^ r~  +  C 1 Korzystając  z  warunków  (1.3)  otrzymamy (3.2)  Xy =  — —— +   [mig—  - - - Z  \   Inn (3.3)  x l   =  —qx+nnq—- - Cał kując  równanie  (2.2), należy  rozważ yć  trzy  przypadki: 1)  h 2   <  s, 2)  h 2   =   s, 3)  h 2   >  s. D la  h 2   <  s  otrzymamy (3.4)  x 2   =  fc - 'l 2 l ( c o sAT + £ 1 si n AT ) + - -T 3 )  Cechą  charakterystyczną  ukł adów  wibro- uderzeniowych  jest  zależ ność  warun ków okresowoś ci  od param etrów  ukł adu  [1], AN AL I Z A  R U C H U   P E WN E G O  U K Ł ADU   WI BR O- U D E R Z E N I OWE GO  57 oraz (3.5)  x 2  = przy  czym  ozn aczon o 3 — t/ ^ —hi  •   £ — v  ?- -• fS  _  e t - 2 l I "^ c |  '  ś mlnnl P odstawiając  warun ki  (1.3) do  (3.3) i (3.5), otrzym am y  p o wykonaniu  koniecznych prze- kształ ceń  wyraż enia  okreś lają ce  prę dkość  uderzenia  oraz  kąt  przesunię cia  fazowego po- mię dzy  sił ą  wymuszają cą  i  przem ieszczeniem : ( 3 . 7 )  v  =   ^ l  +   , £ = (3.8)  sirup — Rv—nnq- \ -  -= D la  h 2  = s  otrzym am y  odpowiedn io (3.9)  x 2  = 2 m  '"  ' "J "  us2' B(\ —e~ 2nnh2 )—d -  4-  — -   I (3.12)  siny  -   J f e , _ ; W ? + ^ r - / J a £ +   ^ - ^ r —i e 2 ™'* i  wreszcie  dla  / ;2 > s: (3.13)  x 2 = -   - - ^ gd zie (3.15) (3- 16)  f , - IA1S 58 Z .  WIŚ N IEWSKI (3.17)  s i n ,  = Rozpatrzym y  z  kolei  przypadek,  gdy  opory  ruch u  ogran iczn ika  2  są   d o  pom inię cia w  porówn an iu  z  wartoś cią   sił y  wymuszają cej. M am y wię c h x   — h 2   =   0. M odel dynamiczny takiego  ukł adu przedstawia schematycznie rys.  3.  R ówn an ia  ruchu  przybiorą   postać (3.18)  ,  x x +q  =  cos(T +q>), (3.19)  x 2 +s 2 x 2   = F o cos  cot• i x • 1 J P 2 Rys.  3.  Schemat  ukł adu  bez  tł umienia Rozwią zaniem  tego  ukł adu  równ ań ,  przy  s  #   1  i  s ( fc =   1, 2, . . . , )  bę dą zwią zki: (3.20) (3.21) (3.22) qt 2 — Znn x+x 0 —cos0, oraz aby  sił a wywoł ana  ciś nieniem powietrza  n a dno tł oka w elemencie sprę ż ystym  zrówno- waż yła cię ż ar masy  uderzają cej,  czyli (3.28)  q>0. Z  (3.26)  otrzymamy  po  wykonaniu  koniecznych  przekształ ceń  nierówność gdzie (3.29) r  SR  .  L 2 R i  D — l - 1,0 1,5 1,0  S Rys.  4.  Zakresy  wartoś ci  parametrów  ukł adu,  odpowiadają cych  rozwią zaniom  okresowym,  przy  n  =   1; Xo«-   - 1;  , u - l;  ^  ==   1;  dla  J?  =   0  i  ii  =   0,5 0,2 0  0,5  1,0  1,5  2,0  s R y s .  5 .  O b s z a r y  i s t n i e n i a  r z e c z y w i s t y c h  w a r t o ś ci  s i n ? )  d l a  «  =   1 ; . v 0 =   — l ; / i =  1 0 ;  <5 =   0 , 0 1 ;  q  =   1 0,4 0,2 0  0,5   - 1,0   1,6   2,0   s R ys.  6.  Obszary  istn ien ia rzeczywistych  wartoś ci  sinip  dla  n—  1;  x 0   =   —  1;  /'•   =   1;  d  —  0, 01;  )0—sin 95, (in deks  0  ozn acza  wielkość  liczon ą   z  uwzglę dn ien iem  docią ż en ia). P o  wyko n an iu  ko n ieczn ych  d ział ań  o t rzym am y (4.6)  / U 2  =   - 5-   { (1 —  T 0)  1  r- =   +   C O S [ $ ( T -   T 0)]sm2nns  I  sm2jins (4.7)  z|5c2  =  0 j i _ l Ł  _  cos [ J ( 2H »-   TO)]}  J ° 2 ^ 7  ~ 7  sin |> ( T -   T 0)] , (4.8)  Zlw =   — ,  ,  _  s i n 2 7 i / w—j——{[ 1 —T 0 + J C O S J ( 2 T O T —  T 0 ) ] t gJ t M + 1+XV  1+ J v 2—  T O ) } +   -  -   •   (2tguins—ctg2nns), l—R  \ f l  — T 0  , .  ,1 (4.9)  A(smw)  =  - x-   ctg7r«i'+  - ——  tgnns  coss(2nn— T 0 )  - 2  \   1+ jt  /  I  s /  T tV%  T 1  ł   I  I  — I —  I T CT 'T f'nP  I 1  I  —i—  "i / r  1  f*' ^ ^ ^ u^ nn  T 0 )  lift  \ i W  przypadku  szczególnym,  gdy  docią ż enie  pojawia  się   na  począ tku  cyklu  ruchu, czyli T 0  ==  0,  zwią zki  (4.6)  —  (4.9)  przyjmą   postać i.  ,/w  ,  0  [\   ,  sinóT  ,, (4.10)  Ax 2   =  —=-   1+ c o s^ T  .-—  (l+scos2nns) (4.11)  Ax 2   =   —  ( 1—J C O S 2 7 I / M ) - T - 4 , - --   -  —sin (^T) (4.12)  Av  =   ' l+R  s(l+R) (2tg7tns—i (4.13)  / I (sin  93)  =  - -  I c t g JT/7J +   I  cos  2J7;«5 - 2sR  I  si: k  1  '  1  1  D l+R 5.  Stabilność strukturalna  układu [(l+ 3i?)ctgra2i'+ i?tgjrw]. Z  punktu  widzenia  zastosowań  technicznych badanie  stabilnoś ci  strukturalnej  przy- nosi  szczególne  korzyś ci,  gdyż  pozwala  okreś lić  wpływ  bł ę dów realizacji  ukł adu  na cha- rakter  ruchu. 64  Z .  WIŚ N IEWSKI D la  okreś lenia  warun ków  stabilnoś ci  strukturaln ej  posł uż ymy  się   m etodą   «dopaso- wywania»  [1],  polegają cą   n a  n adan iu  stał ym  w  cał kach  ogólnych  równ ań  opisują cych ruch  ukł adu, pewnych  mał ych  zaburzeń  (przyrostów),  a  n astę pn ie  n a  porówn an iu  ruchu zaburzonego  i niezaburzonego, przy  liczbie  przedział ów  ruchu  dą ż ą cej  do  nieskoń czonoś ci. N iech  cał ki  ogólne  równań  (3.18)  i  (3.19)  mają   postać (5.1)  '  x,  = X J ( T ;  Q i  C 2 ;cp), (5.2)  x 2   =  x 2 (r;C 3 ;C Ą ). N adajmy  stał ym  C t ;  C 2 ;  C 3 ;  C 4  oraz  ką towi  przesunię cia  fazowego  q>  —  mał e przyrosty  — odpowiedn io:  a v;  ^ v ;  yv; ff„ ; J y .  R uch ukł adu zostanie  zaburzon y,  a cał ki  ogólne  dla  ruchu zaburzonego  w  i'- tym  przedziale  przyjmą   post ać: (5.3)  4 V )  =   4 V ) ( T ;  C , . + «, _ I  ; C a+ / 9, _ iJ 2roi+ +4 Ą ­t). (5.4)  4 V )  ­  *2V)(TS  C 3 + y , _ i ;  C++or,_i;  2 w i + i , ) . Czas  trwania  v­go  cyklu  ruchu  wynosić  bę dzie  ( 2JC T J+ ŚV) zam iast 2n?i. Okreś lmy  róż nicę  pomię dzy  ruchem  zaburzon ym  i  niezaburzonym  n a  koń cu  v- go  prze- dział u  ruchu (5.5)  Axty  = (5.6)  Ax[l>  = oraz  n a  począ tku  ( v+ l) - go  przedział u (5.7)  M (5.8)   p P odobnie  okreś limy  przyrosty  pochodn ych n a  koń cu  v- go  przedział u (5.9)  t  V (5.10) oraz  na  począ tku  (j»+ l)- go  przedział u (5.11)  4 i i f p + 1 ) - (5.12)  Z l4 ; + 1 )  =   x( 2 v+ 1> (O)- .v2(O). Aż eby  ruch  był   stabilny,  okreś lone  wyż ej  przyrosty  speł niać  muszą   nastę pują ce  warun ki: P o  wykonaniu  koniecznych  przekształ ceń  warun ki  (5.13)  sprowadzą   się   do  ukł adu równ ań  róż nicowych,  liniowych,  jedn orodn ych ,  wzglę dem  przyrostów  a v,  / ?„, yv,av,  ń v. Rozwią zań  tego  ukł adu  poszukujemy  w  postaci (5.14)  av  =  fc, e v,  (j„  =   k 2 e v ,  y v   =   k 3 e v ,  a v   =  k A e v ,  A y   =  k s s v , gdzie k x ,  ...,  k 5   —  pewne  stał e. Podstawiają c  (5.11)  do  ukł adu  równ ań ,  wyprowadzon ych  z  warun ków  (5.13),  otrzy- mamy  ukł ad  równań  algebraicznych,  liniowych,  wzglę dem  stał ych  k,,...,k s .  Z  kolei AN ALIZA  RU CH U   PEWNEG O  UKŁADU   WIBRO- UDERZENIOWEGO  65 tworzymy  wyznacznik  charakterystyczny  tego  ukł adu równań. Zerowanie się  tego wyznacz- nika jest  warunkiem  koniecznym istnienia rozwią zań  niezerowych  ukł adu. Rozwijają c  wyznacznik  charakterystyczny,  otrzymamy  równanie  algebraiczne  wzglę - dem  s. Aż eby  ukł ad  był  stabilny,  przyrosty  a V)  pv,  yv,  av,  Av  muszą   dą ż yć  do zera przy  liczbie przedział ów  ruchu v  dą ż ą cej  do  nieskoń czonoś ci. Wynika  stą d  warunek (5.15)  | s | < l ; czyli  moduł y  pierwiastków  równania charakterystycznego  muszą   być mniejsze  od jednoś ci. Warunek  konieczny  i  dostateczny  speł nienia zwią zku  (5.15)  wyznaczymy,  posł ugują c się   twierdzeniem  Schura  [1]. Zgodnie  z tym twierdzeniem  moduł y  pierwiastków  równania charakterystycznego  są   mniejsze  od jednoś ci,  gdy  zachodzi nierówność (5.16)  \ W t \ KHHm  orpaH H ^H TejieM . 3KcnepnMeH Ta.ni.H hie  HccjiefloBaHHH   TaiKHOCTLIO  H enpeptiBH oft  peryjiapoBKH   ^acioT h i  yaapoB BO  BpeiYM   pa6oTbI  yCTpOHCTBa. B  HacTOHmeJi  pa6oTe  pacciwoTpHBaeTCJi  MOflent.,  cocTOHmaa: H3 flByx  Macc5  oAHa  H3 KOTopwx  COOT- Macce  norpe>iiKCH H H ,  BtmH CJienbi  raaBH bie  napaM eTpw  CH CTCMŁI,  Bbffle- ycnoBH H   CTpyKTypiioH   ycToft^HBOCTH.  PaccMOTpen  cny^iań  i- iecTaqHOHapHoro flBH >Kem«ij BH 3- BaH H oro  BHe3anHWM   fleiicTBH eiw  IIOCTOH H H OH   BKemH eii  C H JI M .  H eK o ro p we  pesyjibTaTM   HccjieflOBaHHH H3o6pa>KeHM   n a rp ad n iK ax.  BbiBefleHHBie  B  pa6oTe  4)opiwyjiw  MoryT  H aftra  npH MeneiiH e npH  npoeKTH- BH 6poyflapH wx  norpyH