Analiza odpowiedzi dynamicznych mostu na przejeżdżający pojazd testowy jest powszechną metodą identyfikacji modalnych kształtów drgań konstrukcji. Jednak bezpośredni pomiar przyspieszeń punktów kontaktu kół z mostem (CP) jest w praktyce niezwykle trudny, a często wręcz niemożliwy, ponieważ punkty te poruszają się wraz z pojazdem. W celu rozwiązania tego problemu stosuje się podejście odwrotne: na podstawie pomiarów przyspieszeń nadwozia pojazdu, za pomocą równań równowagi, można obliczyć przyspieszenia punktów kontaktu CP. Metoda ta została potwierdzona w testach polowych, gdzie dwuszynowy pojazd umożliwił skuteczne wyznaczenie przyspieszeń CP na podstawie odpowiedzi układu pojazd–most.

Pomiar odpowiedzi kontaktowych pojazdu pozwala na wyodrębnienie składników modalnych konstrukcji mostu, jednakże uzyskane w ten sposób kształty modalne są zniekształcone przez efekt tłumienia mostu. Kluczowym problemem jest, że dotychczasowe metody, w tym transformacja Hilberta (HT) czy transformacja falkowa (WT), nie eliminują wpływu tłumienia na wyznaczone kształty modalne, co prowadzi do błędów w interpretacji wyników.

W rozwiązaniu tego problemu proponuje się zastosowanie rekursywnej formuły pozwalającej na eliminację efektu tłumienia bez konieczności wcześniejszej znajomości współczynnika tłumienia mostu. Idea opiera się na wykorzystaniu korelacji przestrzennej między dwoma punktami kontaktu — przednim i tylnym kołem pojazdu. Poruszające się koła pojazdu tworzą opóźnienie czasowe między sygnałami odpowiedzi mostu w tych punktach, co umożliwia matematyczną separację wpływu tłumienia.

W praktyce, wyrażenie opisujące amplitudę chwilową składowej modalnej w punktach kontaktu jest funkcją sinusoidalną przesuniętą czasowo oraz zawierającą składnik wykładniczego tłumienia. Poprzez operację ilorazu amplitud zarejestrowanych przy obu kołach pojazdu oraz odpowiedniego uwzględnienia przesunięcia czasowego, można zrekonstruować kształt modalny mostu wolny od zniekształceń wynikających z tłumienia.

Podkreślić należy, że metoda ta nie wymaga znajomości ani oszacowania wartości współczynnika tłumienia mostu — jest to jej ogromna zaleta, znacząco upraszczająca i uwiarygodniająca analizę modalną. Zastosowanie takiego podejścia pozwala na uzyskanie czystszych, bardziej wiarygodnych wyników modalnych, co jest kluczowe w ocenie stanu technicznego mostów i podejmowaniu decyzji konserwacyjnych.

Ważnym aspektem jest również prawidłowa kalibracja parametrów pojazdu testowego. Ich dokładność determinuje wiarygodność wyznaczonych odpowiedzi kontaktowych, a co za tym idzie — skuteczność całej metody. Testy kalibracyjne na równych nawierzchniach lub badania na stole wibracyjnym są niezbędne do potwierdzenia właściwych wartości parametrów pojazdu.

Należy pamiętać, że choć metoda eliminuje efekt tłumienia z modalnych kształtów mostu, to nadal nie zastępuje ona kompleksowej diagnostyki, która powinna uwzględniać także inne aspekty dynamiczne i nieliniowe zachowania konstrukcji. Rozumienie złożoności zjawisk dynamicznych, prawidłowe dobranie modeli teoretycznych oraz ich weryfikacja eksperymentalna stanowią fundament rzetelnej oceny stanu mostów.

Ponadto, istotne jest, aby badacz lub inżynier posiadał świadomość ograniczeń metody, zwłaszcza dotyczących założeń dotyczących liniowości układu oraz równomierności parametrów pojazdu i mostu. W praktyce, zmienne warunki eksploatacji, niedoskonałości pomiarowe oraz efekty nieliniowe mogą wpływać na wyniki i powinny być uwzględnione przy interpretacji danych.

Jak można dokładnie zidentyfikować częstotliwości i kształty drgań mostu z odpowiedzi kontaktowych pojazdu?

Identyfikacja własnych częstotliwości drgań mostu, szczególnie tych o charakterze pionowym oraz skrętnym–ugięciowym, jest kluczowa dla zrozumienia dynamicznego zachowania konstrukcji cienkościennych dźwigarów. W rzeczywistej praktyce inżynierskiej, wyodrębnienie tych częstotliwości ze spektrum odpowiedzi dynamicznej pojazdu napotykającego na most jest wysoce problematyczne, głównie z powodu dominacji częstotliwości własnych samego pojazdu – zwłaszcza kołysania i ruchów pionowych. Jednak poprzez odpowiednie przekształcenie odpowiedzi kontaktowych w miejscu styku kół z nawierzchnią mostu, możliwe jest znaczne osłabienie wpływu częstotliwości pojazdu, co pozwala uwidocznić częstotliwości mostu, szczególnie w wyższych modach.

Odpowiedzi kontaktowe, mimo że nie są bezpośrednio mierzalne w warunkach terenowych, mogą zostać obliczone zarówno metodą elementów skończonych (MES), jak i przy pomocy modeli analitycznych. Porównanie trzech zestawów danych – odpowiedzi uzyskanych analitycznie, odpowiedzi MES oraz tych wyznaczonych metodą odwrotną z odpowiedzi pojazdu – wykazuje niemal idealną zgodność w zakresie częstotliwości skrętnych i pionowych dla obu kół. Tak wysoki stopień zbieżności potwierdza rzetelność zarówno teorii, jak i procedur odwrotnego wyznaczania odpowiedzi kontaktowych.

W przypadku odpowiedzi kontaktowych widoczna jest całkowita eliminacja częstotliwości własnych pojazdu, co skutkuje wyeksponowaniem własnych częstotliwości mostu. Szczególnie interesującym rezultatem jest uwidocznienie częstotliwości wyższych modów, które zwykle są trudne do wyodrębnienia ze spektrum drgań całkowitych. Dodatkowo, różnice w amplitudach pomiędzy odpowiedziami lewego i prawego koła wskazują na różną czułość kół na drgania skrętne mostu – koło położone dalej od osi mostu wykazuje większą czułość na ruchy skrętne, co wynika bezpośrednio z jego oddalenia od osi symetrii przekroju.

Dotychczas separacja pionowych i skrętno–ugięciowych częstotliwości możliwa była jedynie poprzez analizę kształtów modów. Zastosowanie równań transformacyjnych pozwala jednak na precyzyjne oddzielenie tych dwóch typów częstotliwości bez konieczności znajomości form modów. Spektra uzyskane przy użyciu wyprowadzonych równań cechują się wysoką zgodnością z wynikami MES, a oddzielenie typów częstotliwości staje się wyraźne: częstotliwości skrętno–ugięciowe pojawiają się jedynie w spektrum skrętnym, a pionowe – jedynie w spektrum pionowym. Taka rozdzielczość widmowa umożliwia jednoznaczną identyfikację częstotliwości bez zakłóceń wynikających z ich nakładania się.

W kolejnym kroku, użycie transformacji falkowej (WT) do przekształcenia czasowo–częstotliwościowego odpowiedzi dynamicznych pozwala na dokładną rekonstrukcję kształtów modów mostu. Analiza falkowa przeprowadzona z częstotliwością centralną 3 Hz oraz zakresem częstotliwości od 0,5 do 50 Hz umożliwia wykrycie punktów maksymalnych współczynników falkowych, które odpowiadają modom drgań. Linie czasowo–częstotliwościowe oznaczone wartościami maksymalnymi pozwalają zidentyfikować częstotliwości każdego z modów. Po odpowiednim znormalizowaniu danych możliwa jest rekonstrukcja przestrzennego rozkładu kształtów modów – zarówno pionowych, jak i skrętnych.

Analiza wykazuje, że zaproponowana metoda pozwala odtworzyć pierwszy i drugi pionowy kształt modu oraz pierwszy skrętny z wysoką zgodnością w stosunku do wyników analitycznych. Zdolność metody do odwzorowania rzeczywistego rozkładu drgań otwiera drogę do praktycznego zastosowania tych narzędzi w diagnostyce i monitoringu dynamicznego infrastruktury mostowej. W szczególności umożliwia ona nieinwazyjne określenie kluczowych parametrów dynamicznych bez konieczności stosowania gęstej sieci czujników na samym obiekcie.

Warto podkreślić, że uzyskane rezultaty nie tylko wskazują na skuteczność podejścia analitycznego i numerycznego, ale również podważają dotychczasowe założenie o konieczności znajomości kształtów modów dla identyfikacji częstotliwości. Umożliwia to w praktyce szybsze, tańsze i bardziej dostępne procedury diagnostyczne dla istniejących konstrukcji mostowych, które narażone są na zmienne i często trudne do przewidzenia warunki obciążenia dynamicznego.

Dla czytelnika istotne jest zrozumienie, że analiza odpowiedzi dynamicznych na poziomie kontaktowym pozwala na efektywną eliminację wpływu ruchu pojazdu, co z kolei umożliwia precyzyjne uchwycenie właściwości dynamicznych samego mostu. Znajomość lokalizacji kół względem osi mostu wpływa na amplitudę wykrytych modów skrętnych, co może być wykorzystane przy optymalizacji rozmieszczenia czujników w badaniach terenowych. Dodatkowo, zastosowanie transformacji falkowej umożliwia dynamiczną analizę czasu i częstotliwości, co czyni tę metodę szczególnie przydatną w kontekście identyfikacji zmian w czasie, takich jak degradacja konstrukcji czy lokalne uszkodzenia.

Jak wyodrębnić i zidentyfikować właściwości modalne mostu przy użyciu technik analizy kontaktu?

Równanie (15.13) po zastosowaniu reguły de L'Hôpitala oraz wzoru odwrotnego obliczania wprowadza formułę, którą można zapisać w następujący sposób:

kc,ij(t)=m(w,ij¨d2y¨w,ijcdy¨)+y¨s,ijdy¨v,ijw,ijs,ij(y¨ktw,ij2y¨w,ij)k_{c,ij}(t) = m \left( \ddot{w,ij} \frac{d^2ÿ_{w,ij}}{c dÿ} \right) + ÿ_{s,ij} \frac{dÿ_{v,ij}}{w,ij} - s,ij (ÿ k t^2_{w,ij} - ÿw,ij)

gdzie i=F,Ri = F,R, a j=l,rj = l, r. Należy zauważyć, że równanie to zostało wyprowadzone przy założeniu zerowego początkowego przyspieszenia kontaktu u¨c,ij(t)üc,ij(t). Z doświadczeń autorów wynika, że przyjęcie warunku początkowego o niezerowym przyspieszeniu kontaktu ma minimalny wpływ na modalne właściwości mostu, ponieważ te są wyprowadzane za pomocą funkcji GijG_{ij} w wewnętrznej funkcji całkowej (równanie 15.13). Dodatkowo, niezerowa stała całkowa jest szybko wygaszana przez tłumienie pojazdu. Zastosowanie warunku początkowego równego zeru zostało sprawdzone na podstawie wyników pomiarów terenowych, które uzyskano z mostu i porównano je z wynikami obliczeń odwrotnych (Xu et al., 2021).

W kolejnych badaniach pojazd testowy będzie ręcznie wykonany, a każdy jego komponent, np. zestaw kół opon, zostanie przetestowany statycznie i dynamicznie przed montażem. Sztywność pojazdu, rozdzielona równomiernie na cztery koła, będzie mogła zostać określona po uzyskaniu danych na temat wagi i częstotliwości pojazdu próbnego, np. na podstawie testu na płaskiej nawierzchni drogi (Yang et al., 2020f) lub testu na stole wibracyjnym. Ogólnie rzecz biorąc, błędy produkcyjne pojazdu testowego traktowane są jako błąd środowiskowy, który w poprzednich zastosowaniach odpowiedzi kontaktowej pojazdu okazał się akceptowalny (Xu et al., 2022).

Kolejnym krokiem w identyfikacji właściwości modalnych mostu jest wyodrębnienie częstotliwości pionowych i skrętnych mostu. Zgodnie z kinetyczną hipotezą sztywnych przekrojów, odpowiedź skrętną mostu można obliczyć z odpowiedzi kontaktowych, stosując przestrzenną zależność między lewymi i prawymi kołami pojazdu skanującego. Zgodnie z równaniem (15.16):

u¨c,il(t)u¨c,ir(t)=θ¨c,i,i=F,R\ddot{u}_{c,il}(t) - \ddot{u}_{c,ir}(t) = \ddot{\theta}_{c,i}, \quad i = F,R

gdzie bl+brbl + br to szerokości lewego i prawego koła pojazdu. W klasycznej teorii belki, belka jest uznawana za smukłą, jeśli jej szerokość i wysokość są małe w porównaniu z długością. Równanie to jest ogólnie ważne dla takich belek, co ma szczególne zastosowanie w mostach, które mają odpowiednie przegrody (Hambly, 1990). Takie założenie smukłości lub sztywności przekrojów działa również w przypadku wielu cienkościennych belek stosowanych w inżynierii mostowej.

W przypadku, gdy wszystkie cztery koła pojazdu działają na most jednocześnie, a dla n-tej częstotliwości jazdy nπV/Ln\pi V/L jest znacznie mniejsze niż n-ta częstotliwość mostu ωD,b\omega_{D,b}, równanie (15.16) może zostać rozszerzone do postaci:

θc,p(t)=i=F,R;j=l,r(Eˉbz,n,ij,1cos(ωbz,n(tti))+Eˉbz,n,ij,2sin(ωbz,n(tti)))\theta_{c,p}(t) = \sum_{i=F,R; j=l,r} \left( Ē_{bz,n,ij,1} \cos(\omega_{bz,n}(t - t_i)) + Ē_{bz,n,ij,2} \sin(\omega_{bz,n}(t - t_i)) \right) \cdots

Wartości współczynników, takich jak Eˉbz,n,ij,1Ē_{bz,n,ij,1}, Eˉbz,n,ij,2Ē_{bz,n,ij,2} i inne, są wyrażone za pomocą parametrów systemu i stanowią ważną część dalszej analizy dynamicznej mostu.

Oddzielne wyodrębnienie odpowiedzi pionowej i skrętnej mostu można uzyskać dzięki równaniom (15.16) i (15.19), które pozwalają na identyfikację częstotliwości pionowych oraz skrętno-zginających mostu. Te odpowiedzi mogą być analizowane za pomocą transformacji Fouriera w celu uzyskania wyraźnych sygnałów drgań mostu.

Również, identyfikacja współczynnika tłumienia mostu stanowi ważną część tego procesu. Zasadniczo, tłumienie mostu może powodować wygaszenie amplitud odpowiedzi kontaktowych w miarę przejazdu pojazdu. Aby wyodrębnić chwilowe amplitudy odpowiedzi komponentów, stosowane są różne techniki przetwarzania sygnałów, takie jak transformacja Hough'a (HT) oraz transformacja falkowa (WT). Do analizy drgań mostu, ze względu na swoją wysoką efektywność obliczeniową oraz dobrą lokalizację czasowo-częstotliwościową, wykorzystywana jest również transformacja Fouriera okna (transformacja Gabor'a).

Procedura identyfikacji współczynnika tłumienia mostu za pomocą transformacji Gabor'a opiera się na następującym równaniu:

G(τG,nG)=s(t),Γ(t)=+s(t)Γ(tτGTG)exp[i2π(nGFG)t]dtG(\tau_G, n_G) = \langle s(t), \Gamma(t) \rangle = \int_{ -\infty}^{+\infty} s(t)\Gamma(t-\tau_G T_G) \exp[-i2\pi(n_G F_G)t] dt

gdzie s(t)s(t) to zarejestrowany sygnał, a Γ(t)\Gamma(t) jest funkcją jądra. Parametry τG\tau_G oraz nGn_G kontrolują modulację funkcji jądra w dziedzinach czasu i częstotliwości.

Wszystkie te podejścia umożliwiają precyzyjne wyodrębnienie i identyfikację właściwości modalnych mostu, takich jak częstotliwości pionowe, skrętne i współczynniki tłumienia, co stanowi kluczowy element analizy dynamiki mostu.

Jak wzmacniacze poprawiają identyfikację częstotliwości mostów za pomocą pojazdu testowego?

Metoda wykorzystująca pojazd testowy z wbudowanym wzmacniaczem do identyfikacji modalnych parametrów mostów zyskuje na znaczeniu w monitoringu stanu konstrukcji. Tradycyjne podejście z czujnikami montowanymi bezpośrednio na moście, choć dokładne, obarczone jest znacznymi kosztami instalacji i utrzymania oraz brakiem elastyczności sprzętowej. Wobec ogromnej liczby mostów wymagających stałego nadzoru, pojazdy testowe wyposażone w odpowiednie wzmacniacze stają się ekonomiczną i efektywną alternatywą.

Podstawą działania tej metody jest dynamiczne oddziaływanie wzmacniacza, zamontowanego na pojeździe testowym, które pozwala na wyraźniejsze wyodrębnienie częstotliwości drgań mostu z sygnałów zakłóconych przez szumy środowiskowe oraz nierówności nawierzchni. Szczególnie istotne jest dobranie częstotliwości wzmacniacza – gdy stosunek częstotliwości wzmacniacza do częstotliwości mostu (βs,bn = ωs/ωb,n) zbliża się do jedności, następuje wzmocnienie n-tego trybu drgań mostu. Jednak osiągnięcie dokładnego warunku rezonansu (βs,bn = 1) może prowadzić do zakłóceń pomiaru, dlatego ważne jest wyważenie częstotliwości wzmacniacza.

Miejsce instalacji wzmacniacza na moście również wpływa na skuteczność pomiarów. Wzmocnienie częstotliwości drgań zanika w punktach, gdzie nxs/L = 1, 2, ..., natomiast największe efekty obserwuje się dla wartości pośrednich, takich jak nxs/L = 0,5, 1,5 itd. Równocześnie, zastosowanie wzmacniacza na moście powoduje, że energia „wtryskiwana” przez pojazd w system jest proporcjonalnie mniejsza, co pomaga uniknąć dominacji sygnału pojazdu nad sygnałem mostu.

Aby zminimalizować efekt maskowania częstotliwości mostu przez własne częstotliwości pojazdu, rekomendowane jest korzystanie z odpowiedzi kontaktowej pomiędzy pojazdem a mostem. Połączenie danych z odpowiedzi kontaktowej i sygnału wzmacniacza pozwala znacznie zwiększyć widoczność modalnych częstotliwości mostu, w tym trybów wyższych.

Ważnym aspektem jest również sposób tłumienia własnych częstotliwości pojazdu, które mogą zniekształcać sygnał pomiarowy. Wzmacniacz zaprojektowany jako tłumik masowy (TMD) dla pojazdu skutecznie redukuje te zakłócenia, pozwalając na wyraźniejsze wykrycie częstotliwości mostu nawet w warunkach trudnej nawierzchni.

Metody te nie tylko umożliwiają bardziej precyzyjne pozyskanie danych o stanie mostu, ale także stanowią odpowiedź na rosnące zapotrzebowanie na ekonomiczne i mobilne systemy monitoringu infrastruktury. Zastosowanie takich systemów pozwala na systematyczne, ciągłe i mniej inwazyjne monitorowanie mostów, co jest kluczowe dla zapewnienia ich bezpieczeństwa oraz efektywnego zarządzania konserwacją.

Ważne jest, aby czytelnik rozumiał, że pomiar modalnych parametrów mostu za pomocą pojazdu testowego wyposażonego we wzmacniacz to złożony proces, w którym zarówno właściwy dobór częstotliwości wzmacniacza, jak i jego pozycjonowanie mają fundamentalne znaczenie. Ponadto, konieczne jest uwzględnienie interakcji pojazd–most oraz efektów środowiskowych i konstrukcyjnych, które mogą wpłynąć na jakość pozyskiwanych danych. Skuteczne wykorzystanie metody wymaga także zaawansowanej analizy sygnałów i odpowiednich algorytmów separacji modalnych częstotliwości od szumów oraz częstotliwości własnych pojazdu.

Jak stosunek masy pojazdu do mostu wpływa na identyfikację parametrów mostu?

Identyfikacja parametrów modalnych mostów może zostać znacząco zakłócona przez stosunek masy pojazdu do masy mostu. W tym celu, do analizy zastosowano pięć różnych stosunków mas pojazdu do mostu (Rm), przy czym masę pojazdu oraz sztywność zmieniano, zachowując jednak identyczne częstotliwości pojazdu, jak przedstawiono w tabeli 9.2. Wartości stosunku masy Rm wynosiły odpowiednio: 0,1%, 0,5%, 1%, 5% i 10%. Ważne jest, by w tym przypadku zachować stałość częstotliwości pojazdu, ponieważ wpływ stosunku mas pojazdu do mostu stanowi istotny element, który należy dokładnie zbadać.

Wykresy FFT (szybka transformacja Fouriera) dla przyspieszeń kontaktu przednich kół w pięciu rozważanych przypadkach przedstawiono na rysunku 9.13(a), a normalizowane spektra na rysunku 9.13(b). Z tych wykresów wynika, że amplituda częstotliwości mostu rośnie wraz ze wzrostem stosunku masy pojazdu do mostu. Co więcej, dla stosunków masy Rm = 5% i 10% częstotliwość pojazdu (fv) również zaczyna być obecna w odpowiedziach kontaktowych, co jest wynikiem efektu sprzężenia układu pojazd-most.

Dla pięciu rozważanych przypadków, pierwsze częstotliwości mostu (fbD,1) zostały zidentyfikowane i przedstawione w tabeli 9.5. Jak widać, dla mniejszych stosunków masy (Rm = 0,1%, 0,5% i 1%), pierwsza częstotliwość mostu praktycznie nie ulega zmianie. Natomiast dla większych wartości stosunku (Rm = 5% i 10%), częstotliwość ta zaczyna spadać. Warto zauważyć, że dla przypadków z mniejszym stosunkiem masy, wyniki identyfikacji wykazują dużą precyzję, a błędy są stosunkowo małe (poniżej 2%). Natomiast dla większych stosunków masy pojazdu do mostu, precyzja identyfikacji wyraźnie się zmniejsza, a błędy mogą wzrosnąć do 8,5%.

Wnioski płynące z tych wyników są jednoznaczne: stosunek masy pojazdu do mostu ma decydujący wpływ na dokładność identyfikacji parametrów mostu, w tym współczynnika tłumienia. Okazuje się, że dla wartości Rm poniżej 1% wyniki identyfikacji są niezależne od tego parametru. Jednak w przypadku wyższych stosunków masy, wpływ pojazdu na identyfikację staje się istotny. Ważnym aspektem jest fakt, że pojazdy o mniejszej masie zmniejszają wpływ na wyniki, co może być korzystne dla dokładności pomiaru, jednak równocześnie zmniejsza się wzmocnienie drgań mostu.

Projektując pojazd do testów, należy zatem zwrócić szczególną uwagę na masę pojazdu, a także na możliwość łatwego dostosowywania tej masy do różnych potrzeb eksperymentalnych. Na przykład, poprzez dodawanie metalowych płyt nad osią pojazdu, można uzyskać różne warianty masy pojazdu, co pozwala na lepszą kontrolę nad parametrami badania.

Należy jednak pamiętać, że stosunek masy pojazdu do mostu nie jest jedynym czynnikiem, który wpływa na dokładność identyfikacji. Kluczowe znaczenie mają także inne parametry, takie jak geometria mostu, rodzaj materiałów, z których jest zbudowany, czy też stan techniczny samego mostu. Ponadto, w przypadku mostów o kilku przęsłach, wpływ stosunku masy na wyniki identyfikacji może być różny w zależności od liczby przęseł oraz rozmieszczenia wsparć. Analizując mosty wielopręsłowe, należy uwzględnić dodatkowe czynniki związane z reakcjami kontaktowymi i dynamiką układu pojazd-most, które mogą prowadzić do zniekształceń w pomiarach i wynikach identyfikacji.