Od przyspieszenia do prędkości drgań bez błędów

W skrócie

• Dla sinusoidy amplituda prędkości wynosi A/(2πf). Ten skrót nie działa dla dowolnego impulsu przez podstawienie jednej „dominującej” częstotliwości.
• Stały offset 0,005 m/s² po 60 s tworzy sztuczne 0,30 m/s, jeśli zostanie bezpośrednio scałkowany.
• Filtr dolnych częstotliwości nie jest kosmetyką. Może usunąć dryf razem z prawdziwą odpowiedzią konstrukcji.
• PPV wymaga jawnej definicji kierunku, pasma, filtracji i maksimum. Zmiana etykiety z m/s² na mm/s niczego nie liczy.
• Obecna analiza dynamiczna Inclify pracuje na przebiegach przyspieszenia i nie ma ogólnego integratora do prędkości ani gotowego PPV. Pokazuje falę, FFT oraz RMS/MAX w tercjach.

Na ekranie łatwo zrobić prędkość. Wybierasz kanał przyspieszenia, dodajesz blok integracji i dostajesz mm/s. Wykres jest gładki, a liczba ma właściwą jednostkę. To jeszcze nie znaczy, że opisuje ruch konstrukcji.

Każdy niewielki błąd stały w przyspieszeniu rośnie po całkowaniu liniowo z czasem. Powolna zmiana temperatury, przechylenie czujnika, fragment grawitacji, offset przetwornika albo niesymetryczne nasycenie mogą zamienić się w potężny dryf prędkości. Filtr usuwa część problemu, ale robi to kosztem informacji. Granica między błędem i prawdziwym ruchem nie znajduje się sama.

To dlatego w zapytaniu ofertowym pytanie „czy system liczy prędkość?” jest za słabe. Trzeba zapytać, z czego, jak, w jakim paśmie, z jaką niepewnością i czy można odtworzyć wynik.

Trzy wielkości opisują ten sam ruch, ale ważą częstotliwość inaczej

Przemieszczenie x(t), prędkość v(t) i przyspieszenie a(t) łączą pochodne:

v(t) = dx(t)/dt

a(t) = dv(t)/dt = d²x(t)/dt²

W drugą stronę prędkość jest całką przyspieszenia:

v(t) = v(t₀) + ∫[t₀,t] a(τ) dτ

Pojawia się warunek początkowy v(t₀). Dla drgań często przyjmuje się, że średnia prędkość w odpowiednio wybranym oknie wynosi zero, ale to założenie musi pasować do zjawiska. Nie wolno go stosować bezmyślnie do ruchu, który rzeczywiście zawiera wolną składową lub trwałe przemieszczenie.

Dla czystej sinusoidy:

a(t) = A sin(2πft)

amplituda prędkości wynosi:

V = A/(2πf)

amplituda przemieszczenia:

X = A/(2πf)²

Im niższa częstotliwość przy tym samym przyspieszeniu, tym większa prędkość i przemieszczenie. Akcelerometr naturalnie podkreśla wyższe częstotliwości, a całkowanie wzmacnia znaczenie dołu pasma. Ta różnica wyjaśnia, dlaczego ta sama wartość przyspieszenia może prowadzić do innej oceny po przeliczeniu.

Metryka wynika z dokumentu oceny i celu. Zanim zaprojektujesz integrację, rozdziel pytania, na które odpowiadają PPV, RMS, FFT i pasma tercjowe. W analizie maszyn ISO 20816-1 dopuszcza różne wielkości oraz metody, pod warunkiem że są zdefiniowane, a ich ograniczenia opisane. Dla konstrukcji ISO 4866 wiąże pomiar z częstotliwością, czasem trwania i amplitudą odpowiedzi. Nie przenoś procedury maszynowej do budynku tylko dlatego, że obie pokazują mm/s.

Przykład liczbowy: ta sama amplituda przyspieszenia, inna prędkość

Przyjmijmy amplitudę przyspieszenia A = 0,20 m/s².

Dla sinusa 10 Hz:

V = 0,20/(2π × 10) = 0,00318 m/s ≈ 3,18 mm/s

Dla sinusa 2 Hz:

V = 0,20/(2π × 2) = 0,0159 m/s ≈ 15,9 mm/s

Przyspieszenie jest takie samo, a amplituda prędkości przy 2 Hz jest pięć razy większa. Wniosek ma sens tylko dla czystych sinusów. Krótki udar zawiera wiele częstotliwości i faz. Nie można wziąć jego szczytu przyspieszenia, podzielić przez 2π oraz jedną częstotliwość z FFT, a potem nazwać wyniku PPV.

Dla sygnału wieloczęstotliwościowego integrujesz cały skalibrowany przebieg albo każdą składową w domenie częstotliwości z poprawną obsługą fazy, zera i brzegów rekordu. Następnie liczysz maksimum według wymaganej definicji.

Skąd bierze się dryf

Załóżmy stały offset przyspieszenia b = 0,005 m/s². To niewielka liczba wobec krótkiego impulsu. Po prostym całkowaniu przez 60 s daje:

Δv = b × t = 0,005 × 60 = 0,30 m/s = 300 mm/s

Powstała wartość nie opisuje drgań. Jest matematyczną konsekwencją błędu stałego. W realnym torze offset może pochodzić z elektroniki, kalibracji, zmiany temperatury lub składowej grawitacji.

Akcelerometr mierzący DC rejestruje także grawitację. Jeżeli oś obróci się o mały kąt, zmieni się rzut g na tę oś. Przybliżenie dla małego kąta pokazuje skalę problemu: obrót o 0,1° daje zmianę rzędu g × sin(0,1°), czyli około 0,017 m/s². To ponad trzy razy więcej niż offset z naszego przykładu. Dokładny wpływ zależy od typu czujnika, orientacji i sposobu usuwania grawitacji.

Źródła dryfu obejmują także:

• asymetrię dodatnich i ujemnych próbek;
• błąd warunku początkowego;
• przerwę lub duplikat w czasie;
• niezerową średnią wynikającą z obciętego rekordu;
• nasycenie jednego znaku impulsu;
• szum bardzo niskiej częstotliwości;
• skok offsetu po przeciążeniu;
• ruch podstawy albo przewodu.

Usunięcie średniej z całego rekordu pomaga tylko w części przypadków. Jeżeli offset zmienia się w czasie albo rekord obejmuje niesymetryczny impuls, jedno odejmowanie nie naprawia przebiegu. Wielomianowy detrend może z kolei usunąć prawdziwą wolną odpowiedź. Każdy wybór jest elementem modelu pomiarowego.

Filtr górnoprzepustowy rozwiązuje problem i tworzy nowy

Dokumentacja Dewesoft dla integracji dynamicznych przebiegów wskazuje, że offset wzmacnia dryf, dlatego przed lub w ramach integracji stosuje się filtr górnoprzepustowy. Częstotliwość odcięcia decyduje o kompromisie. Wyższe odcięcie stabilizuje wynik szybciej, ale usuwa większą część rzeczywistego niskiego pasma. Niższe zachowuje więcej sygnału, lecz wymaga dłuższego czasu ustalenia i czujnika zdolnego mierzyć ten zakres.

Filtr ma typ, rząd, częstotliwość, charakterystykę amplitudową i fazową. Sam wpis „HP 1 Hz” jest niepełny. Przy częstotliwości bliskiej odcięciu amplituda jest już zmieniona, a filtr przyczynowy przesuwa fazę. Filtracja w przód i wstecz może usunąć przesunięcie fazowe offline, ale korzysta z danych przyszłych i ma inne zachowanie na brzegach. Nie wolno udawać, że jest tym samym co alarm działający online.

W domenie częstotliwości integracja odpowiada dzieleniu przez j2πf:

V(f) = A(f)/(j2πf)

Dla f = 0 wyrażenie jest nieokreślone. Dla bardzo małych f nawet mały szum zostaje silnie wzmocniony. Trzeba jawnie ustalić, co dzieje się ze składową stałą i dolną granicą pasma.

Zerowanie prążka DC nie wystarcza do rozwiązania przecieku z wolnego trendu. Krótki rekord i nieciągłość na jego końcach rozlewają energię na sąsiednie częstotliwości. Funkcja okna ogranicza przeciek, ale także zmienia amplitudę. Parametry muszą być zapisane razem z wynikiem.

Procedura obliczenia, którą da się odebrać

1. Zdefiniuj wynik

Ustal, czy potrzebujesz przebiegu prędkości, wartości RMS, maksimum jednej osi, największej składowej, sumy wektorowej czy PPV według konkretnej metody. Zapisz pasmo i punkt pomiaru. Bez tego nie ma sposobu oceny filtra.

2. Potwierdź jednostkę wejścia

Surowa liczba z przetwornika może oznaczać wolty, g, m/s² albo jednostkę urządzenia. Zastosuj współczynnik kalibracji, czułość i orientację przed integracją. Jeśli wejście jest w g, przelicz je na m/s², a dopiero potem dobieraj skalę wyniku w mm/s. Karta przejścia od wartości surowej do jednostki inżynierskiej pomaga zapisać współczynniki, wersję i odpowiedzialność za ten etap.

3. Sprawdź czas i kompletność

Całkowanie zakłada znane odstępy między próbkami. Luka, duplikat lub jitter zmienia pole pod krzywą. Nie interpoluj brakującego piku bez oznaczenia. Rekord niekompletny nie powinien automatycznie generować końcowego PPV. Dobór czasu, triggera i bufora trzeba odebrać według zasad opisanych w poradniku o próbkowaniu i pre-triggerze.

4. Oceń nasycenie oraz offset

Sprawdź spłaszczone piki, zmianę średniej, zachowanie przed zdarzeniem i powrót do tła. Pre-trigger daje materiał do oceny offsetu oraz ustalenia filtrów. Bez niego algorytm może wymusić błędny warunek początkowy.

5. Ustal filtrację

Zapisz filtr antyaliasingowy toru, filtr górnoprzepustowy integracji, ewentualny dolnoprzepustowy, ich typ, rząd i częstotliwości. Uzasadnij dolne pasmo na podstawie zjawiska oraz dokumentu oceny, nie na podstawie tego, gdzie wykres przestaje dryfować.

6. Całkuj numerycznie

Metoda trapezów, integrator cyfrowy lub operacja w domenie częstotliwości mogą dawać zbliżone wyniki w paśmie, ale różnią się na brzegach i dla składowych wolnych. Nazwij implementację, wersję oraz precyzję. Nie traktuj jej jak neutralnego szczegółu programistycznego.

7. Zweryfikuj sygnałami odniesienia

Użyj sinusów w kilku częstotliwościach, sygnału wielotonowego i impulsu. Porównaj amplitudę oraz fazę z niezależnym torem albo obliczeniem. Dodaj kontrolowany offset, aby sprawdzić, czy procedura nie produkuje sztucznej prędkości.

8. Zachowaj dane potrzebne do powtórzenia

Archiwizuj przebieg przyspieszenia, czas, kalibrację, orientację, filtry, algorytm, wynik i status kompletności. ISO 21289 wskazuje parametry, które należy specyfikować dla pozyskiwania, analizy, raportowania i archiwizacji danych drganiowych. To dobry szkielet wymagań, nawet gdy szczegółową metrykę definiuje inny dokument.

Filtr przed czy po całkowaniu

Kolejność operacji ma znaczenie dla sygnału skończonego i implementacji. Idealne liniowe operacje w nieskończonym, poprawnie reprezentowanym sygnale można opisywać w prosty sposób częstotliwościowy. Rzeczywisty rekord ma początek, koniec, offset, szum i nasycenie.

Filtrowanie przyspieszenia przed integracją ogranicza dolne pasmo, zanim dzielenie przez częstotliwość wzmocni szum. Integrator z wbudowanym filtrem realizuje oba działania jednym torem. Filtrowanie gotowej prędkości może dać inny stan początkowy i przejście brzegowe. Nie wybieraj kolejności po tym, który wykres wygląda spokojniej.

Podobna ostrożność dotyczy tercji. RMS przyspieszenia w paśmie, zintegrowany przez prosty współczynnik dla częstotliwości środkowej, nie jest dokładnie tym samym co RMS przebiegu prędkości po filtracji. Pasmo zawiera zakres częstotliwości, a amplituda nie musi leżeć w środku. Jeśli kryterium wymaga prędkości pasmowej, licz ją zgodnie z jego metodą.

Trzy strategie zakupowe

Bezpośredni przetwornik nie usuwa wszystkich problemów. Ma własne dolne i górne pasmo, tłumienie, czułość poprzeczną i wymagania montażowe. Akcelerometr nie jest automatycznie lepszy, lecz daje elastyczność, jeśli zachowujesz pełny sygnał. Urządzenie brzegowe może być najlepsze w terenie, o ile jego obliczenie jest przejrzyste i kwalifikowane.

Pytania do dostawcy przed zakupem

• Jaka wielkość fizyczna jest mierzona bezpośrednio?
• Czy pełny przebieg wejściowy jest zachowywany i eksportowalny?
• Jaki czujnik, pasmo, zakres, szum i odpowiedź fazową ma tor?
• Jak usuwane są grawitacja, offset i stan przejściowy po przeciążeniu?
• Jaki filtr działa przed próbkowaniem, a jaki przed integracją?
• Jaki jest typ, rząd i odcięcie każdego filtra?
• Jak liczony jest warunek początkowy prędkości?
• Czy algorytm działa online, offline, przyczynowo czy dwukierunkowo?
• Jak system traktuje luki i niekompletny rekord?
• Jak definiuje PPV dla wielu osi i jakiego używa pasma?
• Jakie sygnały oraz tolerancje obejmuje FAT i SAT?
• Czy po zmianie algorytmu można przeliczyć historię bez nadpisania dowodu?
• Czy wynik zawiera wersję profilu i konfigurację?

Odpowiedź „zgodnie z normą” bez numeru, wydania i parametrów nie zalicza żadnego z tych pytań. Wymagane sygnały, tolerancje i dowody należy wpisać do FAT i SAT systemu monitoringu drgań.

Jak wspiera to Inclify i gdzie jest granica

W obecnej implementacji dynamicznej Inclify pełny przebieg jest traktowany jako przyspieszenie. Platforma pokazuje go w czasie, wyznacza FFT oraz 21 pasm tercjowych od 1 do 100 Hz. Profil z filtrami Butterwortha 4. rzędu liczy RMS i MAX przyspieszenia w pasmach. Kalibrację ustawia się osobno dla każdej osi, a profil użyty dla zdarzenia pozostaje zapisany z wynikiem.

Inclify nie ma ogólnego integratora przyspieszenia do prędkości i nie wystawia gotowego PPV. Nie należy więc opisywać wykresu przyspieszenia jako prędkości ani porównywać tercji przyspieszenia z progiem wyrażonym w mm/s. Jeśli projekt wymaga PPV, obliczenie musi powstać w kwalifikowanym torze urządzenia albo poza obecną analizą platformy, z jawnymi filtrami i procedurą walidacji.

Zachowanie pełnego rekordu, czasu, współczynników osi, kompletności i profilu analizy ułatwia niezależne przeliczenie. Konkretny sygnał dynamiczny można wyeksportować do XLSX. To materiał do audytu, nie automatyczne potwierdzenie metryki prędkości.

Alarm dynamiczny Inclify działa na dostępnych wynikach RMS/MAX w tercjach względem obwiedni referencyjnej. Nie zamienia ich w PPV. Na tym kończy się obecny zakres analizy drgań w platformie; wymagania na prędkość trzeba dopisać osobno do projektu integracji.

Ograniczenia, których nie wolno wygładzić

Nie odzyskasz prawdziwej prędkości z przyspieszenia po nasyceniu czujnika. Brakujący fragment także może zawierać maksimum. Filtr ani AI nie odtworzą jednoznacznie informacji, której tor nie zapisał.

Nie ma też jednego „dobrego” odcięcia. Dla maszyny 10 Hz może być rozsądne, a dla budynku usunąć istotną część odpowiedzi. Dokumentacja Dewesoft podaje typowe zakresy jako przykłady aplikacyjne, nie uniwersalne ustawienia do norm budowlanych.

Prędkość obliczona dwoma programami może różnić się mimo tego samego pliku, jeżeli inaczej usuwają średnią, filtrują, obsługują brzegi i skalują FFT. Z tego powodu wynik bez konfiguracji nie nadaje się do porównania dostawców.

Źródła i dalsza lektura

• ISO 4866:2010, Vibration of fixed structures – zasady pomiaru i przetwarzania drgań konstrukcji w odniesieniu do częstotliwości, czasu oraz amplitudy.
• ISO 21289:2008, parameters for acquisition of vibration data – parametry do specyfikacji pomiaru, analizy, raportowania i archiwizacji przebiegów drganiowych.
• ISO 20816-1:2016, measurement and evaluation of machine vibration – przykład ram, w których wielkość pomiarowa i ograniczenia metody muszą być zdefiniowane.
• Dewesoft X Manual, Integration and Double Integration – wpływ offsetu, filtracji górnoprzepustowej, rzędu filtra i skali jednostek.
• PCB Piezotronics, Accelerometer Selection Based on Applications – dobór czujnika oraz wpływ małego offsetu na błąd po integracji.
• NIST Technical Note 334, sampled data and switching logic – próbkowanie, aliasing i ograniczenia rekonstrukcji sygnału.

Decyzja przed zakupem

Poproś dostawcę o jeden plik przyspieszenia oraz wynik prędkości razem z filtrami, kodem wersji i obliczeniem referencyjnym. Potem zmień dolne odcięcie i sprawdź, jak mocno przesuwa się PPV. Jeśli nie można tego odtworzyć, kupujesz liczbę, nie metodę pomiarową. Porozmawiaj z zespołem Inclify, jeśli chcesz ocenić, które dane z istniejącego toru można zachować w platformie i czego brakuje do kwalifikowanego wyniku prędkości.