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 —cos
0, 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,2 0 0,5 1,0 - f.5 2,0 S R ys. 7. Obszary istn ien ia rzeczywistych wartoś ci sin
)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 vgo 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
\