Rezgésdiagnosztika 1. Bevezetés A rezgésdiagnosztika a műszaki diagnosztika egy meghatározott területe. A gépek állapotvizsgálatánál talán a legelterjedtebb vizsgálati módszer a rezgésmérés. Ebben a jegyzetben először a rezgésmérés elméleti alapjait tekintjük át, majd bemutatjuk a rezgésmérés eszközeit. A diagnosztika a gépek, berendezések járművek megbontás nélküli ellenőrzését jelenti. A műszaki állapotvizsgálatot csak abban az esetben szabad alkalmazni, ha a mérési sorozat eredményeként a vizsgált berendezésről, annak műszaki állapotváltozásáról tudunk értékelhető adatokat beszerezni. A jegyzetben esetenként keverednek a rezgéstan és lengéstan fogalmak, de mi minden esetben ugyanazt értjük a két fogalom alatt. Az anyag elsajátítását, az egyes fejezetek végén található ellenőrző kérdések segítik.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
1
2. Rezgéstani alapfogalmak 2.1. Harmonikus rezgőmozgás Rezgésnek nevezünk általában minden olyan fizikai jelenséget, amely az időben periódikusan ismétlődik.
x
A
x
ω
ϕ t
A
R=A T 1.ábra
A rezgések felosztása többféle szempont szerint történhet: a) A rezgést hordozó rugalmas közeg alapján: a rezgések gázokban, folyadékokban, szilárd testekben jöhetnek létre. b) A rezgőmozgást végző anyag alakja szempontjából: •
Húr rezgések (fonál, húr, drót)
•
Sík rezgések (hártya, héj, lemez)
•
Kontinum rezgések (szilárd testek, víz, levegő) különböztethetők meg.
c) Periódikusan változó mennyiségek szerint: Ezek általában fizikai vagy mechanikai mennyiségek. Pl. áram, feszültség, teljesítmény, fordulatszám, nyomaték, nyomás, térfogatáram, út, sebesség. d) Frekvencia szerinti felosztás egyes rezgések esetén: •
Elektromágneses hullámok (pl. rádió, infravörös, látható fény, ultraibolya, stb.
•
Hangrezgések (infrahang, hallgató hangok, ultrahang, hiperhang.)
e) Matematikai megfontolások, a rezgéseket leíró függvények jellege szerint. ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
2
A rezgések legegyszerűbb fajtája a harmonikus rezgés. Igazolható, hogy a függőleges síkban megfelelő fordulatszámú, egyenletes körmozgást végző pontszerű test vetülete és egy rugóra függesztett tömegpont rezgése egyforma, tehát a vetületi mozgás is harmonikus rezgőmozgás. A tömegpont kitérése az időnek szinusz függvénye: x = A ⋅ sin ϕ = A ⋅ sin ωt Itt A az amplitúdó, a rezgő testnek az egyensúlyi helyzetétől mért legnagyobb kitérése. T a rezgésidő (vagy periódus) egy teljes rezgés megtételéhez szükséges idő. T=
2π (s) ω
Az 1 s alatt végzett rezgések száma f f =
1 ω = T 2π
[H z ]
A harmonikus rezgő mozgás kitérését csak akkor adja meg az x = A sin ωt függvény, ha a t időt az egyensúlyi helyzettől számítjuk. Ha az időt tetszőleges helyzetből kiindulva mérjük, akkor t idő múlva a kitérés: x = A sin(ωt + α 1 ) A t=0 időpontban a kitérés: x0 = A . sin α1 E szerint α1 , az un. fázisállandó, megadja a rezgő test helyzetét a t = 0 időpontban.
x
x
A
x0
α1
ωt α1
2.ábra
A tömegpont sebessége:
v=
dx = A ⋅ ω cos (ωt + α 1 ) dt
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
3
A gyorsulás: a=
dv d 2 x = 2 = − Aω 2 ⋅ sin (ωt + α 1 ) dt dt
a=
d 2x = −ω 2 ⋅ x 2 dt
irható továbbá
A harmonikus rezgőmozgás gyorsulása a kitéréssel arányos. Az utóbbi egyenletből a harmonikus rezgőmozgás differenciálegyenlete is felírható: d 2x +ω2 ⋅x =0 dt 2 A d . e – t csupán kinematikai megfontolásból írtuk fel.
2.2. Harmonikus rezgések összetevése a) Két egyirányú, egyenlő frekvenciájú rezgés eredője: x1 = A1 ⋅ sin (ωt + α 1 )
x 2 = A2 ⋅ sin (ωt + α 2 ) A1 és A2 az amplitúdók, ω a körfrekvencia α1 és α2 pedig a rezgés kezdő fázisai. A két rezgés eredő kitérése: x = x1 + x2 = A1 ⋅ sin (ωt + α1 ) + A2 ⋅ sin (ωt + α 2 ) Azonban a sin (α + β ) = sin α ⋅ cos β + cos α ⋅ sin β összefüggés alapján: x = A1 sin ωt ⋅ cos α 1 + A1 ⋅ cos ωt ⋅ sin α 1 + A2 ⋅ sin ωt ⋅ cos 2 + A2 ⋅ cos ωt ⋅ sin α 2 Rendezve: x = ( A1 ⋅ cos α 1 + A2 ⋅ cos α 2 ) sin ωt + ( A1 ⋅ sin α 1 + A2 ⋅ sin α 2 ) cos ωt Legyen: A1 ⋅ cos α 1 + A2 ⋅ cos α 2 = A ⋅ cos α A1 ⋅ sin α1 + A2 ⋅ sin α 2 = A ⋅ sin α akkor x = A cos α ⋅ sin ωt + A ⋅ sin α ⋅ cos ωt = A sin (ωt + α ) Az eredő amplitúdó: A=
A12 + 2 A1 A2 cos(α 1 − α 2 ) + A22
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
4
Az eredő fázisszöge: tgα =
A1 ⋅ sin α1 + A2 ⋅ sin α 2 A1 ⋅ cosα1 + A2 ⋅ cos α 2
A rezgés amplitúdója a legnagyobb, ha
(n = 0, 1, 2,....)
α = α 1 − α 2 = 2π ⋅ n Ebben az esetben ugyanis: cos (α 1 − α 2 ) = 1 A = A1 + A2
x
x
x2
A A2
x1
A1 t
T 3.ábra Ha a fáziskülönbség: x = x1 − x 2 = (2 n + 1)π A = A1 − A2
(n = 0, 1, 2,.......)
b) Két egyirányú, különböző frekvenciájú rezgés eredője. Ez esetben az összetevő rezgések kitérései: x1 = A1 ⋅ sin (ω1 t + α1 )
x 2 = A2 ⋅ sin (ω 2t + α 2 ) Legyen: A1 = A2 = A α1 = α 2 = 0 ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
5
Az eredő rezgés: x = x1 + x 2 = A(sin ω1t + sin ω 2 t ) Azonban a
sin α + sin β = 2 cos
α −β α+β ⋅ sin 2 2
összefüggés felhasználásával az eredő rezgés kitérése: ω − ω2 ω + ω2 x = 2 A cos 1 t ⋅ sin 1 ⋅t 2 2 Látható, hogy az eredő nem harmonikus rezgés. Frekvenciája az összetevő rezgések frekvenciájának számtani közepe; amplitúdója azonban nem állandó, hanem az idő periodikus függvénye: A, (t ) = 2 A ⋅ cos
ω1 − ω 2 ⋅t 2
Ha azonban a frekvenciák közel egyenlők (ω1 ≈ ω 2 ) akkor az az eredő jó közelítésben olyan színuszrezgésnek vehető, melynek körfrekvenciája
(ω1 + ω 2 ) / 2
és amplitúdójának értéke 0 és 2A között viszonylag lassan
változik.
x 2A A t
T 4.ábra A váltakozás (lebegés) frekvenciáját a következőképpen határozzuk meg. Mivel cos (2n + 1)
π =0 2
(n = 0, 1, 2,.....)
Tehát az amplitúdó annyi idő alatt változik zérus értékről ismét zérus értékre, amennyi idő alatt az α szög π értékkel változik meg. ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
6
ω1 − ω 2 ω − ω2 ⋅ t1 = π t2 − 1 2 2
Ebből a lebegés periódusa: T = t2 − t1 =
2π ω1 − ω 2
A lebegés frekvenciája: f =
1 T
c) Két egymásra merőleges, egyenlő frekvenciájú rezgés eredője. Vizsgáljuk meg az általános esetet, amikor a rezgések közötti fáziskülönbség α ≠0.
Ebben az esetben a rezgések egyenlete: x = A1 ⋅ sin ωt
y = A2 ⋅ sin (ωt + α )
Az első egyenletből: 2
x x = sin ωt ; cos ωt = 1 − sin 2 ωt = 1 − 2 A1 A1 A másodikból: y = sin (ωt + α ) = sin ωt ⋅ cos α + cos ωt ⋅ sin ϕ A2 Helyettesítsük ide az első egyenletből kifejezett szögfüggvény értékeket 2
y x x = ⋅ cos α + 1 − 2 sin α A2 A1 A1 Rendezve: x2 y2 2 xy + − ⋅ cos ϕ = sin 2 ϕ 2 2 A1 A2 A1 ⋅ A2 y x x2 − ⋅ cos ϕ = 1 − 2 sin ϕ A2 A1 A1 y2 x2 2 xy + 2 − cos ϕ = sin 2 ϕ 2 A2 A1 A1 ⋅ A2 Ha α =
π 4
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
7
y2 x2 2 xy + 2 − 2 A2 A1 A1 ⋅ A2
1 2
2=
A kapott eredmény egy ellipszis egyenlete, amelynek középpontja a koordináta rendszer középpontjában van és tengelye a vízszintessel 450 –os szöget zár be. (5. ábra)
y
45° x
5.ábra Két egymásra merőleges és különböző frekvenciájú rezgés eredője az összetevők amplitúdójától, a frekvenciák viszonyától és a fáziskülönbségtől függ. Ezek az un. Lissajous-görbék. (6. ábra)
6.ábra ω1 1 = ω2 2
; α =0
ω1 2 = ω2 3
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
; α =0
8
3. A rezgőtest dinamikája A mechanikában a rezgésekkel a lengéstan foglalkozik. A valós szerkezetek rezgéstani (lengéstani) vizsgálatai összetett feladatot jelentenek. A legegyszerűbb lengéstani modell egy rugóból és egy tömegpontból áll. (7. ábra)
y
x
c
c
m
µ=0 m a.
b. 7.ábra
Ha a tömegpontot nyugalmi helyzetből kimozdítjuk, úgy a tömegpontra a kitéréssel ellentétes értelmű rugóerő hat. Ha a rugó visszatérítő erején túl más erő nem befolyásolja a rezgést, akkor örökké tartó mozgás jön létre. Az ilyen lengés a szabad lengés. A rugóerő számítása: F = −D ⋅ x = −
x c
D - a merevség, vagy direkciós állandó c - az un. rugó állandó c=
1 D
mm / N
A tárgyalás során az un. lineáris karakterisztikájú esetet tárgyaljuk.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
9
kemény
F
F
lágy
α x
x
lineáris 8.ábra A 8. ábrán különböző karakterisztikákat tüntettünk fel. Az impulzus tétel felhasználásával d2x dt 2 x F = −D ⋅ x = − c d 2x x + =0 dt 2 m ⋅ c 1 α2 = m⋅c F = m⋅a = m
A behelyettesítések elvégzésével a harmonikus rezgőmozgás differenciál egyenletét kapjuk. d 2x +α2 ⋅ x = 0 dt 2 A differenciál egyenlet általános megoldása: x = A ⋅ sin (α ⋅ t + α1 ) Hasonló megoldásra jutottunk tisztán kinematikai meggondolásokból is.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
10
3.1. A harmonikus rezgések fajtái: A rugózó elem eltérő tulajdonsága miatt különböző lengések lehetnek. Így beszélhetünk: egyenesvonalú lengésről torziós lengésről hajlító lengésről
3.1.1. Egyenesvonalú lengés. Ezt a lengést longitudinális lengésnek is nevezzük. A tömegre ható visszatérítő erő a rúd rugalmasságából ered és rugóállandója: c=
x l = F AE
A lengés körfrekvenciája: 1 = m⋅ c
∅
α=
A⋅ E m⋅l
d
m
x
l
9.ábra
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
11
3.1.2.Torziós lengés. Ezen lengésnél a lengő tömeg változó értelmű körmozgást végez. A rugózást a tengely valósítja meg. A rugóállandó ebben az esetben az egységnyi nyomaték által okozott szögelfordulás lesz. l I p ⋅G
∅
c0 =
d
l J m D
Mt
10.ábra
A körfrekvencia:
α=
1 J ⋅ c0
=
Jp ⋅ G l ⋅ c0
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
12
3.1.3. Hajlítólengés. Ha a lengő tömeg a rudat hajlításra veszi igénybe. A rugóállandó ebben az esetben az egységnyi erő által okozott lehajlás.
∅
c=
l3 3I ⋅ E
d
m f
l 11.ábra A lengés körfrekvenciája: α
1 31E = ml 3 m⋅ c
3.2. Ellenőrző kérdések 1. Ismertesse a lengési modellt. 2. Mi a rugókarakterisztika? 3. Hogy írhatjuk fel a harmonikus lengés differenciál egyenletét? 4. Ismertesse a harmonikus lengés foronómiai görbéit. 5. Hogyan határozza meg az azonos frekvenciájú két harmonikus lengés eredő amplitúdóját és fázis szögét? 6. Mikor jön létre a lebegés? 7. Mekkora a lebegési idő? 8. Ismertesse a harmonikus rezgés fajtáit. 9. Írja fel a különböző harmonikus lengésekre az α és T értékeit.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
13
4. A harmonikus lengés csillapítása 4.1. Nedves csillapítás A valóságban a rezgések amplitúdója egyre jobban csökken, ezt a jelenséget csillapodásnak nevezzük. A csillapodásnak az oka a belső súrlódás, a közegellenállás és egyéb veszteségek. A csillapító erő a kitéréssel ellentétes és a sebességgel arányos (ez un. nedves csillapítás). F = −k ⋅
dx dt
A csillapított lengés modelljét a 12. ábra szerint válasszuk.
c
m
k
12.ábra A csillapított lengés alapegyenlete: m
dx x d2x dx = −D ⋅ x − k =− −k 2 dt c dt dt
Bevezetve a következő jelöléseket: D 1 = =α2 m c⋅ m
;
k =β 2m
β - a lengés csillapítási tényezője
A csillapodó rezgőmozgás differenciálegyenlete: d 2x dx + 2β ⋅ + α 2 ⋅ x = 0 2 dt dt
A differenciál egyenlet. egy partikuláris megoldása: x = e λt λt ekkor: x = λ ⋅ e
;
x = λ 2 e λt
Behelyettesítve a differenciál egyenletbe:
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
14
(
)
e λ t λ 2 + 2β ⋅ λ + α 2 = 0 λ1 , 2 = −β ± β 2 − α 2 A gyakorlat szempontjából fontos eset, ha: β 〈α
így
α 2 − β 2〉0
γ = α2 − β2 Ezek után a partikuláris megoldás: x1 = e − β t ⋅ cos γt x 2 = e − βt ⋅ sin γt
Az általános megoldás tehát: x = e − βt (C1 ⋅ cos γt + C 2 ⋅ sin γt ) A premfeltételek: x (0 ) = 0
x (0 ) = v0 Ha:
x(0) = 0
akkor
C1 = 0
x = − β ⋅ e − βt (C 2 ⋅ sin γt ) + e − βt ⋅ (C 2 ⋅ cos γt ⋅ γ ) x (0 ) = v0 x (t ) =
akkor
C2 =
v0 γ
v 0 − βt ⋅ e ⋅ sin γt γ
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
15
x
A
T/2
T/2
T
a1 x1
a3
a5
x3
T
a6
a4
t
a2 A T= 2π γ
T
-Ae- βt
13.ábra
T=
2π 2π = γ α −β2
áll
a1 C ⋅ e − βt t = 2 − β (t +T ) = e β T a2 C 2 e 1 A logaritmus dekramentum: υ = ln
a1 = βT a2
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
16
4.2. Száraz csillapítás Ha a csillapítóerő iránya a mozgással ellentétes, nagysága állandó, akkor súrlódásos csillapításról beszélünk.
c
m
0 x
A
µ a2
a1
B
14.ábra Ha a tömegpont A.-tól – 0-n keresztül B-ig halad a súrlódó erő a – x irányba mutat. Amíg a tömegpont B-től A-ig halad a súrlódó erő +x irányú. A-tól – B-ig: −
x d 2x − Fs = m 2 = 0 c dt
B-től A-ig: x d 2x + Fs − m 2 = 0 c dt Fs = µ ⋅ m ⋅ g
−
Vezessük be: Fs =
x0 c
x0 = c⋅µ ⋅m⋅ g A differenciál egyenletbe irható: −
Vezessük be:
x ± x0 d 2x −m 2 = 0 c dt
z = x ± x0 d2z +z=0 dt 2 z +α2z = 0
m
A differenciál egyenlet megoldása:
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
17
z = a ⋅ sin (α ⋅ t + α 1 )
x = a ⋅ sin (αt + α 1 ) ± x 0
x
a1
2π T 0= α
T0
T0/2
a1 a3
a4
x0
α α
t x0 a2
a2
T0/4 T0/4 T0/4 T0/4 T0/4 T0/4 T0/4 T0/4 T0/4 T0/4
15.ábra
Az amplitúdók a munkatétellel határozhatók meg. a22 a12 − = − Fs (a1 + a 2 ) 2e 2c (a 2 − a1 ) ⋅ (a 2 + a1 ) = −2 Fs ⋅ c(a1 + a2 ) a 2 = a1 − 2 ⋅ x 0
Példa 1. A vízszintes érdes síkon μ= 0,05,
G = m . g = 50 N súlyú testet c = 1,5 mm/N
rugóállandójú rugóhoz kapcsolunk. Hány teljes lengés után fog megállni, ha az első kitérés a1 = 30 cm?
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
18
Megoldás: Fs = μ . m . g = 0,05 . 50 = 2,5 N x0 = Fs . c = 3,75 mm a1 = 30 mm
; a2 = a1 – 2x0 = 22,5 mm
a3 = a2 – 2x0 = 15 ; a4 = 7,5 a5 = 7,5 – 2. 2,75 = 0 T=
2π = 2π m ⋅ c = 2π ⋅ 5 ⋅ 2,5 ⋅ 10 −3 = 0,545 sec α
A mozgás: ?T = 2T = 1,08 sec
x
a1 a3
+x0 -x0
t a2
a4
T
T
16. ábra
A tömeg az x0 szélességű sávon belül leáll, ugyanis az
x visszatérítő erő már kisebb, c
mint a súrlódó erő.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
19
5. A harmonikus lengés gerjesztése (rugón keresztül történő gerjesztés) O1 O2
A’ O2 A1
A O1
η
B1
c
η
m O3 x
ω
B
x
r
O3
ωt 17.ábra
A deformáció miatt x ≠ η η = r ⋅ sin (ωt ) A rugóerő:
F=
η−x c
η−x − m⋅ x = 0 c r ⋅ sin (ωt ) − x − m⋅ x = 0 c r ⋅ sin (ωt ) x − −x=0 m⋅c m⋅c 1 α2 = m⋅c 2 & & x + α ⋅ x = r ⋅ α 2 ⋅ sin (ωt ) A differenciál egyenlet megoldását keressük: x = K ⋅ sin (ωt )
& x&= − K ⋅ ω 2 sin (ωt ) − K ⋅ ω 2 ⋅ sin (ωt ) + α 2 ⋅ sin (ωt ) = rα 2 ⋅ sin (ωt ) − K ⋅ω 2 + α 2 ⋅ K = r ⋅α 2 K=r
α2 α 2 −ω2
Az un. rezonancia függvény: ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
20
R=
K α2 = 2 = r α −ω2
1 ω 1− α
2
A differenciál egyenlet megoldása tehát
x =r⋅
1 ω 1− α
R=
2
⋅ sin (ωt − α 1 )
K2 1 = r 1-ξ2
R
1 0
ω ξ= α
1 18.ábra
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
21
5.1 Hajlításra igénybevett tengely kritikus fordulatszáma: FcF = m ⋅ ( x + e) ⋅ ω 2 x 2 − + m( x + e)ω = 0 c x 2 2 xω + e ⋅ ω − =0 m ⋅c xω 2 + eω 2 − x ⋅ α 2 = 0
(
ω
ϕ
)
x ω 2 − α 2 + eω 2 = 0 x =e
ω2 1 = e⋅ 2 2 2 α −ω α −1 ω
e o
Fr
S
Fc
X ecos ϕ=e
19.ábra
Rezonancia esetén: ω =α 2π ⋅ n n ≈ 60 9,55 1 9,55 = n kr = 9,55 ⋅ m ⋅c m⋅c n=
c –a tömeg középpontjába helyezett 1 N erő által okozott x elmozdulással egyenlő.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
22
6. Gerjesztett és csillapított lengések Vizsgáljunk egy színuszfüggvény szerint változó erővel a rugón keresztül gerjesztett nedves csillapítású lengőrendszert. F0 ⋅ sin ωt =
r sin ωt c
k
c m
20.ábra
Írjuk fel a kinetikai alapegyenletet:
− m ⋅ x − kx −
x + F0 ⋅ sin ωt = 0 c
Rendezve: mx + kx +
x = F0 sin ωt c
Bevezetve: k 1 ; α2 = 2m m⋅c 2 2 x + 2β x + α ⋅ x = rα ⋅ sin ωt β=
A differenciál egyenlet megoldása a homogén és partikuláris megoldások összegeként nyerhető: x h = C 2 ⋅ e − βt ⋅ sin γt
x p = K ⋅ sin (ωt − ε ) x = xh + x p
A K és az ε meghatározható a Gümbel-féle forgó vektorokkal. (21. ábra)
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
23
r F0= c
ω
m⋅K⋅ω
2
y
ε k⋅K⋅ω
m⋅K⋅ω
r F 0= c ε
K c
2
k⋅K⋅ω x K =mα2⋅K c
K
21.ábra
A visszatérítő erőrektor a K amplitúdóval kifejezve:
K = K ⋅ m ⋅α 2 c
A tömeg tehetetlenségi erejének vektora: K . m ω2 A csillapító erő vektora: K.K.ω A gerjesztő erő vektora a kitérés K vektorával ε fázisszöget zár be és az
r sin (ωt) c
gerjesztő erőnek maximális értékével egyenlő, r c
Tehát:
Most vizsgáljuk ezen erők egyensúlyának feltételeit: ∑ Xi = 0 ∑ Yi = 0
r kK ⋅ ω − sin ε = 0 c r K ⋅ mω 2 + cos ε − K ⋅ mα 2 = 0 c
; ;
E két egyenletből megfelelő rendezés után az ismeretlen ε és
K
értékei
meghatározhatók. r sin ε = K ⋅ k ⋅ ω c r cos ε = K ⋅ m α 2 − ω 2 c
(
)
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
24
A két egyenlet osztásával: tgε = χ=
k ⋅ω ξ =χ 2 2 m α −ω 1−ξ 2
(
)
1 mα
ω α
ξ=
;
Emeljük a két egyenletet négyzetre és adjuk őket össze: 2
[
2
(
2 r r 2 2 sin ε + cos ε = (k ⋅ K ⋅ ω ) + K ⋅ m α − ω c c
[
2
(
r 2 2 2 2 2 2 = K k ω + m α −ω c
)]
2
)] 2
és ebből K=
r c
(
k 2ω 2 + m 2 α 2 − ω 2
)
2
=
(
r ⋅ m ⋅α 2
)
m2 α 2 − ω 2 + k 2ω 2
K – a gerjesztett és nedves csillapítású lengés amplitúdója. K2 – a csillapítás nélküli gerjesztett lengés amplitúdója. K < K2 Ugyanis K=
rα 2
(α
2
−ω
)
2 2
kω + m
2
〈 K=
α2 α 2 −ω2
Ha ω = α - akkor rezonancia jön létre. Ekkor a K értéke maximális, de nem végtelen nagy. K rez = K max =
m ⋅ r ⋅α 2 k ⋅ω
A rezonancia függvény az előzőekből ismert (18. ábra): R=
K = r
1
(1 − ξ )
2 2
+ χ 2 ⋅ξ 2
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
25
κ=0
K R= r
6,0 5,0
κ=0,2
4,0 3,0 κ=0,5
2,0
κ=1,0
1,0
κ=2,0 0
0,5
1,5
1,0
2,0
2,5
3,0
ξ= ω α
ε
κ=1,0 κ=2,0
π 2
B A
0
0,5
κ=0,2 κ=0,5 κ=1,0 κ=2,0
κ=0,0
C
π
κ=0,5 κ=0,2 κ=0,0 1,5
1,0
2,0
2,5
3,0
ω ξ= α
22.ábra A rezonancia ábra alatt feltüntettük az ε = f (ξ ) fázisfüggvényt is.
Az ábrából is kitűnik, hogy a rezonanciának (ξ = 1) megfelelő fázis szög: ε=
π 2
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
26
6.1. Ellenőrző kérdések 1. Ismertesse a főbb csillapítási fajtákat. 2. Ismertesse a száraz csillapítást. 3. Írja fel a száraz csillapítás differenciál egyenletét. 4. Mekkora a száraz csillapítás lengésideje és csillapítás tényezője? 5. Írja fel a nedves csillapítás mozgás egyenletét. 6. Írja fel a nedves csillapítás körfrekvenciáját, lengésidejét. 7. Mi a logaritmus dekrementum? 8. Ismertesse a rezonancia függvényt. 9. A rugón keresztüli gerjesztésnél mutassa be a kinetikai egyensúly vektoros ábráit. 10. Ismertesse a hajlított tengely kritikus fordulatszámát.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
27
7. Többtömegű lengőrendszer Ha nem egyetlen tömeg képez egy vagy több rugóval lengőrendszert, akkor többtömegű lengőrendszerről beszélünk. Az egyes tömegek szabad mozgásukat illetően többféleképpen lehetnek beépítve, így beszélhetünk szabad és kötött rendszerekről. A szabad rendszer minden tömege a hozzákapcsolt rugó által korlátozottan bár, de mégis szabadon mozoghat. (23. ábra) m1
m2
c1
m3
c2 23.ábra
A kötött rendszer egy vagy több rugója egyik végének merev kapcsolata miatt szabad mozgásában gátolva van. (24. ábra) m1
c1
m2
c2
m0 24.ábra 7.1. Többtömegű szabad rendszer A következőkben három tömegű szabad rendszert vizsgálunk és az így nyert eredményeket általánosítjuk. (25. ábra) m1 m2 c1 F1 x1
F1
m3
c2 F2
F2
x2
x3
25.ábra A c1 rugó deformációja (x1 – x2 ), a rugó erő pedig: ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
28
F1 =
x1 − x 2 c1
A c2 rugó deformációja (x2 – x3), a rugóerő pedig: F2 =
x 2 − x3 c2
Írjuk fel mindhárom tömegre külön-külön a kinetikai egyensúlyi egyenleteket: − m1 ⋅ x1 −
x1 − x2 =0 c1
x − x3 x1 − x 2 =0 − m2 − x 2 − 2 − c1 c2 x − x3 − m3 ⋅ x 3 + 2 =0 c2 Az átalakítások után a következő egyenletet kapjuk ,melyet lengési főegyenletnek is nevezünk. x1 (m1 + m2 + m3 ) + x ( IV ) (m2 ⋅ m3 ⋅ c 2 + m1 ⋅ m2 ⋅ c1 + m1 ⋅ m 3 ⋅ c1 + m1 ⋅ m3 ⋅ c 2 ) +
+ x1(VI ) (c1 ⋅ c 2 ⋅ m1 ⋅ m2 ⋅ m3 ) = 0
x1 (m1 + m2 + m3 ) + x IV [m1 ⋅ m2 ⋅ c1 + m1 ⋅ m3 (c1 + c 2 ) + m2 ⋅ m3 ⋅] +
+ x VI ⋅ m1 ⋅ m2 ⋅ m3 ⋅ c1 ⋅ c 2 = 0 A megoldást keressük: x = A ⋅ sin (αt ) A behelyettesítések után az un. karakterisztikus egyenletet kapjuk: m1 + m 2 + m 3 − α 2 [m1 ⋅ m 2 ⋅ c1 + m1 ⋅ m3 (c1 + c2 ) + m 2 m 3 ⋅ c 2 ] + + α 4 ⋅ m1 ⋅ m 2 ⋅ m3 c1 ⋅ c 2 = 0
Az egyenletben az ismeretlen α legmagasabb kitevője 2 (n - 1), ahol n a lengő tömegek száma. A gyökök száma ugyancsak (2n – 1). A karakterisztikus egyenlet
n
lengőtömeg esetében a fenti példa mintájára
mindenkor könnyen felírhatjuk. Az egyenlet felépítése ugyanis a következő: Az első tag a lengő tömegek összegéből áll, és pozitív előjelű. A második tag negatív előjelű és az α2-tel szorozzuk a tömegek és a közöttük lévő rugóállandó szorzatának összegét. A harmadik tag ismét pozitív előjelű és α4 szorozva a tömegekkel és a köztük lévő rugóállandókkal. ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
29
A gyökök közül a legkisebb α a rendszer alap lengésének körfrekvenciája. Ha a megfigyelés kezdetén már volt a tömegeknek kilengése, akkor x = K ⋅ sin (α ⋅ t + α1 ) A tömegek kilengésének amplitúdói: K1, K2 és K3 a kitérésekre a következő függvényeket kapjuk: x1 = K1 ⋅ sin (αt + α 1 )
x 2 = K 2 ⋅ sin (αt + α 1 ) x 3 = K 3 ⋅ sin (αt + α1 ) Behelyettesítve a kinetikai egyenletekbe, majd rendezve az egyenleteket, valamint az egyenletek összeadása után kapjuk: − α 2 (m1 ⋅ K 1 + m2 + m3 ⋅ K 3 ) = 0 Tetszőleges több tömeg esetére is fennáll: ∑ m, ⋅ K , = 0
Ez a súlyponttétel azt mutatja, hogy szabad lengés esetén a lengőrendszer súlypontja helyben marad.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
30
7.2. A többtömegű kötött rendszer A kötött rendszer a szabad rendszerből úgy származtatható, hogy a megfogás helyén lévő tömeget végtelen nagynak tekintjük (26. ábra) m0
m1 c1
m2 c2
m3 c3
26.ábra Az ennek figyelembevételével felírt karakterisztikus egyenletből a körfrekvenciák meghatározhatók. A számítások igazolják, hogy a karakterisztikus egyenlet α-nak kettőnél magasabb kitevőjű hatványai az alaplengés numerikus értékét alig befolyásolják. Így, ha az α2 –nél magasabb hatványú tagokat elhanyagoljuk, akkor 5 – 10 %-nál nagyobb hibát az α1 = αmin meghatározásában nem követünk el. A karakterisztikus egyenlet az elhanyagolásokkal:
(m0 + m1 + m2 + m3 ) − α 2 [m0 ⋅ m1 ⋅ c1 + m0 ⋅ m2 (c1 + c 2 ) + + m0 ⋅ m 3 (c1 + c 2 + c3 ) + m1 ⋅ m2 c 2 + m1 ⋅ m3 (c 2 + c3 ) + m2 ⋅ m3 ⋅ c 3 ] = 0 Az egyenletet m0 –val osztva, kapjuk: 1 − α 2 [m1 ⋅ c1 + m 2 (c1 + c 2 ) + m3 (c1 + c 2 + c3 )] = 0 Az m0 = ∞ , hiszen így vettük figyelembe a rendszer kötöttségét. Dunkerley formulája A kötött rendszer legkisebb körfrekvenciájának meghatározására alkalmas. A 27. ábra szerint a háromtömegű kötött rendszert három egytömegű lengőrendszerre bontjuk.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
31
m1 c1
αI m1 c1
m2 c2
m3
m2
c3
c1
+
c2
αII m3 c1
+
c2
+
c3
αIII 27.ábra
1 1 1 1 = 2 + 2 + 2 = m1 ⋅ c1 + m2 (c1 + c2 ) + m3 (c1 + c 2 + c3 ) 2 α 1 α I α II α III
Érdemes megfigyelni, hogy a karakterisztikus egyenletből is ezt az eredményt kapjuk.
7.2.1.Többtömegű rendszer hajlító és torziós lengése. A többtömegű rendszerben a hajlító és torziós lengést vizsgáljuk. A gyakorlatban azonban mindkettő egyszeri fellépése is megtörténhet. Az egyszerűsítés érdekében azonban csak a hajlító, illetve a torziós lengés viszonyait tárgyaljuk.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
32
7.2.2. Hajlító lengés A 28. ábrán feltüntetett modell két tömege hajlító lengést végez. Az m1 tömeg m1 . g súlyereje a saját függőlegesében y1 lehajlást okoz, tehát a rugóállandó: c1 =
y1 m1 ⋅ g
1 m1
2 m2 y1
1N y1
1N y2 28.ábra Az m2 . g súlyerő a 2-es függőlegesben y2 lehajlást okoz, tehát a rugóállandó e helyen: c2 =
y2 m2 ⋅ g
Írjuk fel a két tömegre a Dunkerley formulát:
1 1 1 = 2 + 2 = m1 ⋅ c1 + m2 ⋅ c 2 2 α I α I α II és ebből a rendszer legkisebb körfrekvenciája: α 1 = α min =
1 m1 ⋅ c1 + m2 ⋅ c 2
7.2.3. Torziós lengés
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
33
A 29. ábrán feltüntetett kéttömegű torziós lengési modell megoldásához használjuk fel, a szabad rendszert. J1
J2 ∅d
l
29.ábra
Példa 2. Határozza meg a vázolt tengely legkisebb saját frekvenciáit hajlító és torziós lengések esetére. m1 = 6 kg m2 = 9 kg J 1 = 0,03 kgm
2
J 2 = 0,045 kgm 2 d1=∅200
d2=∅300
d=60
300
600
300
1N 1200 30.ábra
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
34
A tengely statikai adatai: d 4π = 63,6 cm 4 64 E = 210 GPa ; I=
IE ⋅ y1 =
;
I p = 2 I = 127,2 cm 4
G = 81 GPa
(
)
(
93 0,9 4 ⋅ 0,1125 2 + 0, 225 2 + 0, 2252 + 4 ⋅ 0,1125 2 6 6
)
0,2025 = 0,0015 mm / N 63,6 ⋅ 10 −8 ⋅ 210 ⋅ 10 9 1 = 210,9 1 / s α 1 = α min = 15 ⋅ 0,0015 ⋅ 10 − 3 c1 = y1 =
c0 =
0,6 = 0,0045 ⋅ 10 − 3 1 / Nm 127,2 ⋅ 810
α 1 = α min =
0,075 = 3513,6 1 / s 0,03 ⋅ 0,045 ⋅ 0,0045 ⋅ 10 − 3
Egyenesvonalú lengésre nyert karakterisztikus egyenletét, mely szerint:
(m1 + m2 ) − α 2 (m1 ⋅ m2 ⋅ c1 ) = 0 A felírt egyenletben a tömeg helyett a tehetetlenségi nyomatékot, a rugóállandó ebben az esetben az egységnyi nyomatékhoz tartozó szögelfordulást kell helyettesíteni. c0 =
l 32 ⋅ l = 4 I p ⋅G d π ⋅G
a legkisebb körfrekvencia: α 1 = α min =
J1 + J 2 J1 ⋅ J 2 ⋅ c 0
7.3. Ellenőrző kérdések 1. Hogyan osztályozzuk a mechanikai lengéseket (rezgéseket)? 2. Írja fel egy háromtömegű szabad rendszer karakterisztikus egyenletét. 3. Írja fel a kötött rendszer karakterisztikus egyenletét. 4. Ismertesse a Dunkerley formulát. 5. Ismertesse a többtömegű szabad rendszer hajlító tengését. 6. Ismertesse a többtömegű szabad rendszer torziós lengését.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
35
8. Rezgések csoportosítása A rezgőmozgások felosztását a 31. ábra szerinti csoportosításban vizsgáljuk. A stacionárius (időben állandósult) rezgések és a nem stacionárius rezgések közötti különbség elsősorban abban van, hogy a stacionárius rezgés esetén az átlagos jellemzők (pl. négyzetes középérték…stb.) állandó. A nem stacionárius rezgések esetén, különböző időpillanatokban mérve a jeleket, átlagos jellemzőik különbözőek lesznek.
REZGÉSE
Stacioner
Determisztikus
Periódikus
Harmonikus
Nem stacioner
Sztochasztikus
Folyamato vagy stacionáris sztochasztikus
Tranziens
Kvázi dikus
Nem harmonikus
31. ábra Rezgések felosztása A stacionárius determinisztikus rezgések esetén a jelek pillanatnyi értéke előre megadható, míg a stacionárius sztochasztikus rezgések esetén csak a statisztikai jellemzők (átlagos érték, négyzetes középérték….stb.) adhatók meg előre. Periodikus rezgésről akkor beszélünk, ha egy test mozgásciklusa T időnként ismétlődik. Jellegzetes példája a 31. ábrán látható test mozgása, mely harmonikus rezgőmozgás. A gyakorlatban azonban a rezgések csak rendkívül ritkán tiszta harmonikus mozgások. Egy tipikusan nem harmonikus, periodikus rezgés például egy olyan szinuszos összetevőkből eredő rezgés,mely a belső égésű motor hengerfedelén gyorsulás-idő függvény értékeként mérhető ERFP-DD2002-HU-B-01 PROJECT 4. MODUL Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
1
(32. ábra). Kváziperiodikus rezgés jellemzője, hogy minden esetben összetett rezgésként adódik, de az alaprezgések egymásnak nem harmonikusai, így az eredő rezgés T periódus szerint nem ismétlődik. Ilyen rezgés adódik, ha két egyébként harmonikus alaprezgést összegzünk, de ω1/ω2 hányados nem racionális szám. A gyakorlatban rendszerint olyan véletlenszerű (sztochasztikus) rezgések lépnek fel, ahol a szabálytalan mozgásciklusok sohasem ismétlik meg tökéletesen önmagukat (33. ábra). Ilyen rezgés egzakt leíráshoz elméletileg hosszú megfigyelési időre lenne szükség, ami felesleges, mivel más, például a statisztikai jellemzők állandósága miatt már viszonylag rövid idő alatt is megfelelő pontossággal mérhető a rezgés erőssége. a
t T 32.ábra Nem harmonikus, periodikus rezgés x
t
33.ábra ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
2
A nem stacionárius sztochasztikus rezgések jellemzésére az egyszerű átlagolás nem alkalmas, mivel ezek statisztikai pereméterei is változnak. Ez gyakorlatilag azt jelenti, hogy ugyanazt a mérést bizonyos idő után (1-2 óra) megismételjük, más átlagos értéket kapunk. Tipikus példa rá egy útkereszteződésnél mérhető rezgésszint, mely időjárástól, évszaktól, forgalomtól függően változik. Ezért a méréseket meghatározott időszakonként meg kell ismételni és a mérési eredményeket egy újabb átlagolásnak kell alávetni (34. ábra). A tranziens rezgések (35. ábra) szintén gyakran előfordulnak a mindennapi életben. Ezek olyan ütésszerű, rövid időtartalmú rezgések, melyek jellemzésére, időben bonyolult lefutásuk miatt, a környezetükre kifejtett hatásuk szolgál. Ezt időbeli lefutásuk mellett elsősorban energiatartalmuk határozza meg.
(4)
x
(3)
x
(2)
x
idõ-átlagolás csoport-átlagolás (1)
x
t1
t2
t3
idõ-átlagolás; többszöri átlagolás 34.ábra Nem stacionárius sztochasztikus rezgések
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
3
x
t
35.ábra Tranziens rezgés
x (t)
A
x (f)
x=A⋅sinωt t
f0
frekvencia
f
f0= 1 T
T 36.ábra
Frekvenciaspektrum x(f) függvénye harmónikus rezgőmozgás esetén.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
4
8.1. Rezgések jellemzése A rezgések leírása időbeli alakjuk mellett egyéb jellemzőik is használatosak. Ilyen jellemző a rezgés frekvencia tartomány szerinti ábrázolása. Ez az ábrázolás sokkal egyszerűbb és mégsem tartalmaz információveszteséget. Ez különösen bonyolultabb rezgések elemzésénél hasznos, ahol több összetett függvény együttes értelmezése már áttekinthetetlen, a frekvenciaspektrum szerinti ábrázolás viszont jól kiemeli az egyes komponenseket. A frekvenciaspektrum x(f) függvénye matematikailag az x(t) időfüggvény Fouriertranszformáltjaként számítható. Nézzük meg ezt egy harmonikus rezgőmozgás esetén (36. ábra). Ennek matematikai alakja:
x( f ) = F {x(t )} =
∞
∫ x(t )e
− iωt
dt
−∞
= A, ha f = f 0 x( f ) = F {A ⋅ sin ω } = 0; ha f ≠ f 0
Természetesen ez nemcsak kitérés esetére alkalmazható, a sebesség és a gyorsulás is ábrázolható frekvenciatartomány szerint A harmonikus rezgőmozgás a kitérés, a sebesség, vagy a gyorsulás legnagyobb értékével és a körfrekvenciával, a rezgés számmal vagy egy rezgés idejével egyértelműen megadható, ezért álszerű a csúcsértékek mérése. (A, v0 .a0) Ha viszont a rezgés nem tiszta szinuszos, akkor nem jellemezhető a csúcsértékkel, mert ez semmit sem mond arról, hogy milyen a mozgás jellege a két csúcsérték között. Ilyenkor valamivel többet mond a számtani, illetve a négyzetes középérték (az effektívérték). A számtani középértéket a gyakorlatban kevésbé használják. Mind két értéket későbbiekben tárgyaljuk. Rezgésnél alapvető követelmény, hogy a magas szintértékek kijelzése mellett, az alacsony szintek közötti finom különbség is kimutatható
legyen.
Ezért
logaritmikus
skálát
szokás
alkalmazni
az
ún.
decibelszámítás szerint: dB = 20 ⋅ lg
a mért a ref .
ahol: ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
36
amért a rezgés jellemző mért értéke; aref
a referencia vagy vonatkoztatási szint.
Természetesen itt is tetszőleges jellemzőt használhatunk. A
0
dB érték nem
rezgéstelen állapotot jelent, hanem azt, hogy a mért szint és a referenciaszint megegyezik. A decibelskála legnagyobb előnye a széles tartomány (nagy dinamika) átfogása úgy, hogy az alacsony szintértékeknél is megmarad a megfelelő felbontás. A decibelskáláról leolvasott szintek abszolút értékei csak a referenciaértékek ismeretében határozhatók meg. A referenciaértékekre léteznek szabványos ajánlatok (2.1. táblázat), de ezt a megfelelő dinamika érdekében nem mindig tartják be. 2. táblázat
Magyarország Franciaország Németország Oroszország
A referenciamennyiségek használatos értéke néhány országban vr m/s 5.10 -8 1.10 -8 5.10 -8 4,18.10 -8
ar m/s2 1.10-6 1.10-5 1.10-6 3,102.10 -4
A csillapítás kifejezésére is szokás használni a decibelszámítást. Ebben az esetben azonban nem egy meghatározott referenciaszinthez viszonyítunk, hanem a bemenő és a kimenő rezgéseket vizsgáljuk, pl.:
dB = 20 ⋅ lg
a bemenő a kimenő
Ezek után nézzünk néhány példát a különféle rezgések frekvenciatartomány szerinti ábrázolására. A nem harmonikus periodikus rezgés képe látható a 2.15. ábrán. A kváziperiodikus rezgés esetén a frekvencianövekmény nem állandó, mert a rezgések nem egymás harmonikusai (38.ábra), A sztochasztikus rezgések esetén a rezgés frekvencia spektruma folytonos, mivel mindenféle rezgésösszetevő előfordul benne (39. ábra) A rezgések elemzése esetén a csúcsértékek (amplitúdó maximum) mellett célszerű átlagos szintek definiálása is. Ez elsősorban - a nem tiszta szinuszos rezgések leírásakor pontosabb értelmezést jelent. Átlagos rezgésérték (számtani középérték): ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
37
T
x átl = ahol: xátl
1 x(t )dt T ∫0
számtani középérték;
x(t) a rezgésre jellemző függvény
x (t)
a (f)
t
f 0= 1 T
f
f 2= 1 T2
T2 T1 37.ábra
x
f 38.ábra Kváziperiodikus rezgés frekvenciaspektruma
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
38
x (t)
x (f)
t
t
39.ábra
Négyzetes középérték (effektív középérték): x RMS =
1T 2 x (t ) dt T ∫0
ahol: xRMS Root Mean Square = RMS, négyzetes középértéket jelent. x(t)
a rezgésre jellemző függvény.
A függvények jellemzésére használjuk még az alak- és a sisaktényezőt is. Az alaktényező (form): Ff =
x RMS x átl
ahol: xRMS négyzetes középérték; xátl
számtani középérték.
A sisaktényező (crest): Fc =
x csúcs x RMS
ahol: xcsúcs az amplitúdó maximális értéke; xRS
négyzetes középérték.
Természetesen itt is értelmezhető a sebesség és a gyorsulás értékekre az átlagolás, az alak- és a sisaktényező meghatározása. Határozzuk meg ezeket az értékeket egy harmonikus rezgés esetén (40. ábra)
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
39
x (t)
r
xátl
xRMS t
40.ábra Átlagos rezgésértékek értelmezése harmonikus rezgés esetén
A kitérés függvény:
x(t ) = A ⋅ sin (ωt + α )
ahol: A = xcsúcs amplitúdó maximum; α=0 Számtani középérték: π/2 1 2 x átl = 4 ⋅ xcsúcs ∫ sin ωtdt = ⋅ xcsúcs 2π π 0 Négyzetes középérték: xRMS
2 2⋅ = x π csúcs
π /2
∫ sin ωtdt = 2
0
xcsúcs 2
Az alaktényező:
x RMS π = x átl 2⋅ 2 ahol: xRMS = négyzetes középérték xátl = számtani középérték Ff =
A sisaktényező: Fc =
xcsúcs = 2 x RMS
ahol: xcsúcs az amplitudó max. értéke xRMS
négyzetes középérték
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
40
9. A rezgésmérés eszközei A rezgésméréshez a 41. ábrán látható eszközök szükségesek. Az érzékelők általában rendelkeznek egy un. szeizmikus tömeggel, melynek rezgéseit maga az érzékelő alakítja elektromos jellé. Kijelző
Érzékelő
Elemző/Analiz.
Erősítő
Regisztráló
41.ábra
Természetesen vannak mechanikus elven működő műszerek is, de ezek pontatlanok. Célszerűbb a gépek mechanikai rezgéseit villamos jellé átalakítani és a villamos jelet mérni és elemezni. 9.1. Relatív rezgésmérők A relatív rezgésmérés lényege, hogy a mérőeszközt egy – a tárgy rezgésétől független – un. fix ponthoz rögzítjük. Ilyenkor azt vizsgáljuk, hogy a lengő mozgást végző tárgy hogy mozdul el a fix ponthoz képest. ( 42. ábra).
lengõ mozgást végzõ test m c nyugalomban lévõ test
42.ábra 9.2. Abszolút rezgésmérők Ha fix pontunk nincs, abszolút rezgésmérőt választunk. Ilyen mérések esetén a mérendő tárgyhoz egy egyszabadságfokú lengőrendszert csatlakoztatunk, melynek tömege a vizsgált tárgy lengőmozgása a következtében gerjesztett lengést végez. (43. ábra)
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
37
a ház mozgása
y
π
c
u
2π
3π
4π
ωt
4π
ωt
a tömeg mozgása
m π q
k
2π
3π
a felrajzolt diagram
π ϕ
2π
3π
fázisszög
4π q(t)=u(t)-x(t)
ωt
43.ábra
A mérendő lengés: x = A ⋅ sin ωt
Tekintettel arra, hogy a ház is elmozdul a műszer tömege ahhoz képest csak relatív elmozdulást végez, melynek bizonyos ? fáziskésése van. A felrajzolt diagram tehát a q (t ) = Q ⋅ sin (ωt − ϕ ) ϕ - a műszer fázisszög hibája. Ezen felírásnál feltételezzük, hogy a műszernek a mérés kezdetekor jelentkező saját lengése már lecsillapodott. A műszer által észlelt lengés amplitúdó (Q) és a mérendő lengés amplitúdójának (A) viszonya, az út – nagyítás függvény Nu =
Q = A 1 − λ2
(
)
λ2 2
+ 4 ⋅ C 2 ⋅ λ2
c – a rendszer csillapítására jellemző un. Lehr-féle szám C=
k 2m ⋅ α
λ=
f ω = fm α
; az un. frekvencia viszony.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
38
A műszer frekvenciája: fm =
1 1 ⋅ 2π m⋅c
Az út – nagyítás függvényből (44. ábra) látható, hogy útérzékelés csak nagy λ értékek esetében lehetséges. Nu D=0
2,0
0,4 0,5
1,5
1,0 0,6 0,7 1,0
0,5
0
1
2
3
4
f λ= α
mérési tartomány 44.ábra Ha a mérés pontatlansága nem lehet nagyobb mint ± 5 %, akkor a csillapatás mértékét 0,6….0,7 között kell választani és a mérendő frekvencia a műszer saját frekvenciájának 1,6szorosánál kisebb nem lehet. Mivel a csillapítást nehéz a kívánt határok közé beszabályozni, ezért a műszer mérési tartományát úgy választják meg, hogy λ nagyobb legyen mint 2,5 – 3. Az útérzékelő műszereket tehát nagy λ ill. kis α vagyis nagy m tömeg jellemzi. Ipari viszonylatban legfeljebb nagyon nagy vizsgált frekvenciák esetében adódnak elfogadható műszertömegek. Ezért az ipar csaknem kizárólag gyorsulásérzékelőket alkalmaz, melyek kisméretű tömegeket igényelnek. A mérendő lengésgyorsulás amplitúdója: a0 = A ⋅ ω 2 ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
39
A gyorsulás – nagyítás függvénye: N gy =
Q⋅α 2 = A ⋅ω 2 1 − λ2
(
)
2
1 + 4 ⋅ c 2 ⋅ λ2
Ez a függvény grafikusan (45. ábrán) látható. Nu 2,0
D=0
0,4 0,5
1,5
1,0 0,6 0,7 1,0
0,5
0
1
2
3
4
f λ= α
mérési tartomány 44.ábra Ha a gyorsulásérzékelő műszer pontatlanságát ± 5 %-on belül akarjuk tartani, akkor a csillapítási tényező c = 0,6 – 0,7 határok között kívánatos tartani. A gyorsulás mérése a saját frekvencia 0,7 –szereséig végezhető.
9.3. A gyorsulásérzékelők gyakorlati kivitele A gyakorlatban az alábbi típusokat különböztetik meg: - piezoelektromos, amely tovább osztható un. kompressziós és nyíráson alapuló típusra. - nyúlásmérő bélyeges, - szénoszlopos, - elektrodinamikus Részletesebb jellemzőik a következők: ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
40
a) Piezoelektromos A mérés fizikai alapja a következő, a piezokristályokra ható erő hatására a kristályokra felvitt fegyverzetek között villamos töltés keletkezik, melyet ellenálláson levezetve feszültség keletkezik. A gyakorlatban két fajta típus terjedt el az un. kompressziós és a nyíráson alapuló típus (46. ábra)
Elõterhelõ rugó Tömeg Piezoelektromos elem
Elektromos kimenet
Alap Rögzítõ furat Rezgés hatóerõ
46.ábra A piezoelektromos megoldás a gyorsulások mérésére igen alkalmas, rendkívül kis tömege miatt saját frekvenciája nagy, ezért rendkívül kis méretei miatt szinte minden ipari igényt kielégít. A gyorsulás érzékelő tömege legfeljebb egytizede annak a tömegnek, melyet mérünk.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
41
10% Viszonylagos érzékenység dB
30 20 10 Töltés vagy feszültségérzékenység 0
Használható frekvenciatartomány ∼0,33
10 0,0001
0,001
0,01
1
0,1
f/f0
10
47.ábra b) Nyúlásmérő bélyeges A berendezés elvi vázlata 47. ábrán látható. Lényege egy házban laprugókra felerősített m tömeg. Ha a ház a rezgő berendezésre történt felhelyezés, vagy a felerősítés után a berendezéssel együtt rezeg, az m tömeg a tehetetlenség miatt a gyorsulás irányával ellentétesen, de a gyorsulás nagyságával arányosan kitér a laprugók deformációját okozva. A laprugók legjobban igénybevett részén felragasztott 1, 2, 3 és 4 jelű nyúlásmérő bélyegek a lehajlással arányos villamos jelet szolgáltatnak. (48. ábra) csillapító 1 2 3
1 4
2 UKí
3
m
4 Ut 48.ábra c) Szénoszlopos Ismeretes, hogy az egymással érintkező szénlapok vezetőképessége összenyomásuk esetén az átmeneti ellenállás csökkenése miatt nő. A 49. ábrán látható megoldásnál a tömeget lemezrugók tartják, amelyeknek másik oldala csavarokkal előfeszített kb ? 6 x 0,2 mm méretű szénERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
42
pasztillákra támaszkodik. Csillapítás érdekében a ház és a tömeg közötti részben olaj van. A műszert csak úgy, mint a nyúlásmérőbélyeges mérés esetén Wheatstone – hídba kapcsolják.
m y
49.ábra
50.ábra
Szénoszlopos gyorsulás érzékelő
Elektrodinamikus rezgésérzékelő
d) Elektrodinamikus Az elektrodinamikus rezgésérzékelő fő eleme egy mágneses térben mozgó tekercs (50.ábra), amelyben a mozgás következtében feszültség indukálódik. Az indukált feszültség B lengésindukáció D tekercs átmérő és n menetszám esetén: U i = −n ⋅ B ⋅ D ⋅ π ⋅ v Az elektrodinamikus rezgésmérők kimeneti feszültsége csakis a rezgés v sebességétől függ. 9.4. Ellenőrző kérdések: 1. Hogyan osztályozhatjuk a rezgéseket? 2. Hogyan jellemezhetők a rezgések? 3. Ismertesse a rezgésmérés elvét. 4. Ismertesse az abszolút rezgésmérés elvét. 5. Milyen frekvencia tartományban használhatók a gyorsulásérzékelők. 6. Ismertesse a gyorsulásérzékelők gyakorlati kivitelét.
ERFP-DD2002-HU-B-01 PROJECT 4. MODUL
Ipari hátterű alternáló képzés előkészítése a Gépészmérnöki Szakon
PDF created with pdfFactory Pro trial version www.pdffactory.com
43