Izolacja drgań zmienia układ mechaniczny: wprowadza nowe częstotliwości własne, ruch względny, drogi przenoszenia siły i sposób odpowiedzi na udar. Rozwiązanie skuteczne w jednym paśmie może wzmacniać odpowiedź w pobliżu rezonansu albo dopuścić do uderzenia w ogranicznik. Ocena wymaga zatem jednoczesnego śledzenia ruchu chronionego zespołu, skoku zawieszenia i sił działających w interfejsach.

Ostrzeżenie EOD. Publiczny opis może szczegółowo wyjaśniać fizykę, modelowanie, metrologię i proces kwalifikacji. Nie powinien natomiast ujawniać częstotliwości własnych, rozmieszczenia mocowań, poziomów prób, map odpowiedzi ani słabych punktów konkretnego współczesnego systemu bojowego. Nieznanego lub uszkodzonego sprzętu wojskowego nie wolno poruszać, wzbudzać, mocować do stanowiska, zasilać ani zwalniać jego blokad. Miejsce zabezpiecza się i przekazuje właściwemu zespołowi EOD; element sprężysty może przechowywać energię także wtedy, gdy nic się widocznie nie porusza.

1. Funkcja chroniona i mapa źródło–droga–odbiornik

Wymaganie „zmniejszyć drgania” nie określa jeszcze, co uznaje się za powodzenie. Czujnik może tracić dokładność wskutek obrotu obudowy, złącze może ulegać mikropoślizgowi pod wpływem siły cyklicznej, płytka drukowana może osiągać niebezpieczne odkształcenie na lokalnej postaci drgań, a przewód może zostać uszkodzony przez nadmierny skok między dwoma punktami mocowania. Jednakowa wartość przyspieszenia obudowy nie oznacza jednakowego zagrożenia dla tych funkcji.

Opis funkcji chronionej powinien wskazywać wielkość fizyczną związaną ze stanem granicznym: dopuszczalne przemieszczenie lub obrót, odkształcenie, siłę w złączu, chwilową utratę ciągłości elektrycznej albo błąd pomiaru. Do tego dochodzą pasmo, czas działania, kierunek, temperatura, konfiguracja i obciążenie statyczne. Dopiero takie wymaganie pozwala ocenić, czy właściwą odpowiedzią będzie zmiana źródła, podatne mocowanie, usztywnienie odbiornika, tłumik czy absorber strojony.

Mapa źródło–droga–odbiornik porządkuje fizyczną stronę problemu. Źródłem może być ruch podstawy, niewyważenie maszyny, pulsacja ciśnienia, pole akustyczne albo krótkie zdarzenie udarowe. Droga prowadzi przez zaprojektowane podpory, lecz również przez konstrukcję nośną, przewody, rury, taśmy cieplne, połączenia wyrównawcze, blokady transportowe i przypadkowe styki. Odbiornikiem jest konkretny element wraz z pełnioną funkcją. Punkt, w którym umieszczono akcelerometr, pozostaje jedynie miejscem obserwacji.

Ta mapa zapobiega pozornym sukcesom. Spadek przyspieszenia na pokrywie nie dowodzi ochrony płytki, jeżeli obudowa zaczęła się obracać i silniej pobudza jej lokalny mod. Dobra charakterystyka izolatorów nie opisuje kompletnego zespołu, gdy napięty przewód tworzy równoległą drogę o porównywalnej sztywności. Jeżeli źródło znajduje się wewnątrz izolowanej skrzynki, podatne oddzielenie skrzynki od podstawy może utrudnić przejęcie reakcji zamiast ograniczyć odpowiedź odbiornika.

Wynikiem pierwszego etapu jest więc para: mierzalna funkcja chroniona oraz model dróg, który wyjaśnia, którędy może dotrzeć obciążenie. Każda dalsza transmitancja musi mieć nazwane wejście i wyjście. Bez tego liczba opisana tylko jako „poziom drgań” nie prowadzi do jednoznacznej decyzji.

2. SDOF przy wymuszeniu podstawą

Najprostszy model izolowanego zespołu ma jeden stopień swobody (single degree of freedom, SDOF): masę $m$, liniową sprężynę o sztywności $k$ i tłumik lepki o współczynniku $c$. Niech $x(t)$ oznacza bezwzględne przemieszczenie masy, a $y(t)$ — przemieszczenie podstawy. Sprężyna i tłumik reagują na różnice położenia oraz prędkości, dlatego równanie ruchu ma postać:7,16

$$m\ddot{x}+c(\dot{x}-\dot{y})+k(x-y)=0.$$

Do kontroli skoku wygodniejsza jest współrzędna względna

$$z=x-y.$$

Ponieważ $\ddot{x}=\ddot{z}+\ddot{y}$, po podstawieniu otrzymuje się

$$m\ddot{z}+c\dot{z}+kz=-m\ddot{y}.$$

Prawa strona pokazuje bezpośrednio, że wejściem jest przyspieszenie podstawy. To równanie służy zarówno do wyznaczania odpowiedzi harmonicznej, jak i do całkowania przebiegu czasowego udaru. Z rozwiązania $z(t)$ wynikają odkształcenie izolatora, wykorzystanie prześwitu i chwila ewentualnego kontaktu z ogranicznikiem.

Parametry porządkujące model wynoszą

$$\omega_n=\sqrt{\frac{k}{m}},\qquad f_n=\frac{\omega_n}{2\pi},\qquad \zeta=\frac{c}{2\sqrt{km}},$$

gdzie $\omega_n$ jest kołową częstotliwością własną układu nietłumionego, $f_n$ odpowiadającą jej częstotliwością, a $\zeta$ bezwymiarowym względnym tłumieniem. Dla wymuszenia harmonicznego o częstotliwości kołowej $\omega$ wprowadza się stosunek

$$r=\frac{\omega}{\omega_n}=\frac{f}{f_n}.$$

Częstotliwość własna nie jest własnością izolatora rozpatrywanego bez obciążenia. Ta sama podpora współpracująca z inną masą tworzy inny układ, a zmiana sztywności dynamicznej pod wpływem temperatury lub amplitudy również przesuwa $f_n$. Model SDOF zachowuje wartość dopóty, dopóki masa może być traktowana jak bryła sztywna, ruch odbywa się głównie w jednym kierunku, a charakterystyki $k$ i $c$ pozostają wystarczająco liniowe w analizowanym zakresie.

3. Ruch bezwzględny, ruch względny i siła przekazywana

Dla podstawy poruszającej się harmonicznie z amplitudą $Y$ transmitancję ruchu bezwzględnego definiuje się jako stosunek amplitudy masy $X$ do amplitudy podstawy:

$$T_x=\frac{X}{Y} =\sqrt{\frac{1+(2\zeta r)^2} {(1-r^2)^2+(2\zeta r)^2}}.$$

Ponieważ przyspieszenia harmoniczne $x$ i $y$ są mnożone przez tę samą wartość $\omega^2$, stosunek amplitud przyspieszeń ma w tym szczególnym przypadku tę samą wartość:

$$\frac{|\ddot X|}{|\ddot Y|}=T_x.$$

Nie uprawnia to do utożsamiania ruchu i przyspieszenia w dziedzinie czasu. Równość dotyczy stosunku amplitud dwóch przebiegów sinusoidalnych o jednej częstotliwości. Dla udaru lub procesu szerokopasmowego każde pasmo ma własne $r$, amplitudę i fazę, a szczyt przemieszczenia nie wynika z samego szczytu przyspieszenia.

Transmitancja ruchu względnego jest inną funkcją:

$$T_z=\frac{Z}{Y} =\frac{r^2} {\sqrt{(1-r^2)^2+(2\zeta r)^2}}.$$

Wysoko powyżej rezonansu $T_x$ może być małe, podczas gdy $T_z$ dąży do jedności. Masa pozostaje wtedy prawie nieruchoma w przestrzeni, a podstawa porusza się pod nią; izolacja ruchu bezwzględnego zostaje okupiona skokiem zbliżonym do amplitudy podstawy. Właśnie $Z$, po uwzględnieniu położenia statycznego i tolerancji, porównuje się z prześwitem.

Siła przekazywana przez sprężynę i tłumik do podstawy wynosi

$$F_p=kz+c\dot{z}.$$

W stanie ustalonym można ją znormalizować względem amplitudy $m\omega^2Y$, czyli siły bezwładności potrzebnej do wymuszenia ruchu masy sztywno związanej z podstawą:

$$T_F=\frac{|F_p|}{m\omega^2Y}.$$

Dla tego modelu, harmonicznego wymuszenia podstawą i tak przyjętej normalizacji otrzymuje się $T_F=T_x$. Jest to wynik bilansu sił, a nie dowód zamienności wielkości. $X/Y$ opisuje ruch, natomiast $F_p$ obciąża interfejs i konstrukcję nośną. Zmiana mianownika albo rozpatrywanie wymuszenia siłowego prowadzi do innej funkcji transmitancji.

Moc i energia wymagają kolejnych definicji. Średnia moc rozpraszana w idealnym tłumiku lepkim podczas ruchu harmonicznego wynosi

$$\overline{P}_c=\frac{1}{2}c\omega^2 Z^2.$$

Sprężyna okresowo magazynuje i oddaje energię; w stanie ustalonym nie rozprasza jej średnio w czasie. Amplituda siły, moc średnia, energia rozproszona w cyklu i energia całego zdarzenia udarowego odpowiadają więc na różne pytania. Nie można ich zastępować jedną „transmisją energii”.

Dla zdefiniowanej wyżej transmitancji $T_x$ początek izolacji oznacza $T_x<1$. Z równania wynika warunek

$$r>\sqrt{2}.$$

Warunek ten odnosi się do harmonicznej transmitancji bezwzględnego ruchu masy względem ruchu podstawy. Nie stanowi uniwersalnej granicy dla skoku, mocy, odpowiedzi przejściowej ani dowolnie zdefiniowanej siły.

W pobliżu $r=1$ tłumienie zmniejsza wzmocnienie i skraca wybrzmiewanie. Przy małym $\zeta$ szczyt jest wysoki, a przy $r=1$

$$T_x=\frac{\sqrt{1+4\zeta^2}}{2\zeta}.$$

Daleko powyżej rezonansu pojawia się odmienny kompromis. Dla niezerowego tłumienia licznik jest zdominowany przez $2\zeta r$, więc w przybliżeniu $T_x\sim 2\zeta/r$. W modelu bez tłumienia spadek byłby szybszy, około $1/r^2$. Zwiększenie $\zeta$ ogranicza zatem odpowiedź rezonansową, lecz może podnieść transmitancję ruchu i znormalizowanej siły w paśmie leżącym znacznie powyżej rezonansu. Wyprowadzenie odpowiedzi częstotliwościowej liniowych układów drugiego rzędu i układu wzbudzanego ruchem podstawy przedstawia również podręcznik Hallauera.16

Wymuszenie podstawą i wymuszenie siłowe są różnymi zadaniami

Przy wymuszeniu podstawą wejściem jest $y(t)$, a siła bezwładności $-m\ddot y$ pojawia się po przejściu do współrzędnej względnej. Typowym celem jest izolacja odbiornika: ograniczenie bezwzględnego ruchu, przyspieszenia lub innej funkcji chronionej przy jednoczesnym zachowaniu skoku. Gdy natomiast na masę działa siła harmoniczna $F(t)=F_0\sin\omega t$, a podstawa pozostaje nieruchoma, równanie ma postać

$$m\ddot x+c\dot x+kx=F(t),$$

a amplituda przemieszczenia wynosi

$$X=\frac{F_0}{\sqrt{(k-m\omega^2)^2+(c\omega)^2}}.$$

Siła przekazana nieruchomej podstawie ma amplitudę $F_t=X\sqrt{k^2+(c\omega)^2}$. Po znormalizowaniu przez $F_0$ otrzymuje się

$$\frac{F_t}{F_0}= \sqrt{\frac{1+(2\zeta r)^2}{(1-r^2)^2+(2\zeta r)^2}}.$$

Wyrażenie ma tę samą postać algebraiczną co $T_x$ dla ruchu podstawy, lecz opisuje inną parę wejście–wyjście. W izolacji źródła celem jest ograniczenie reakcji przekazywanej do otoczenia; w izolacji odbiornika — odpowiedzi chronionego zespołu. Równość dwóch funkcji w idealnym SDOF nie przenosi się automatycznie na układ wielopunktowy, podatną podstawę, wymuszenia rozłożone ani nieliniowe drogi równoległe. Dlatego raport musi nazwać zarówno źródło wymuszenia, jak i mianownik transmitancji.

Tłumienie lepkie, histeretyczne i tarciowe nie są wymienne

Współczynnik $c$ jest modelem siły proporcjonalnej do chwilowej prędkości. Jest wygodny, ponieważ zachowuje liniowe równanie ruchu i daje prostą miarę $\zeta$, ale rzeczywisty elastomer, połączenie śrubowe czy splot metalowy nie muszą zachowywać się lepko. Dla liniowego opisu harmonicznego materiału lepkosprężystego stosuje się zespoloną sztywność

$$K^*=K'+jK''=K'(1+j\eta),$$

gdzie $K'$ jest częścią magazynującą, $K''$ — stratną, a $\eta=K''/K'$ współczynnikiem strat. Dla szczytowej amplitudy ruchu względnego $Z$ energia rozproszona w cyklu wynosi $\Delta W=\pi K''Z^2$. Równoważny współczynnik lepki w jednym punkcie pracy można zapisać jako $c_{eq}=K''/\omega$. Nie jest on stałą elementu: zmiana częstotliwości, temperatury, amplitudy, docisku lub liczby cykli może zmienić zarówno $K'$, jak i $K''$.18,20

Tarcie Coulomba daje pętlę zależną od amplitudy, próg ruszenia i możliwy stick–slip; zastąpienie go jednym $c_{eq}$ może odtworzyć energię jednego cyklu, ale niekoniecznie fazę, małoamplitudową odpowiedź ani przejście między przywieraniem i poślizgiem.19 Podobnie współczynnik strat wyznaczony z jednej pętli nie opisuje automatycznie udaru. Równoważność modelu trzeba więc przypisać do konkretnego pasma, amplitudy, temperatury i wielkości wyjściowej, a poza tym zakresem wrócić do zmierzonego prawa siła–przemieszczenie–prędkość.

4. Skok, ugięcie statyczne i ogranicznik

Dla pionowego, liniowego układu podtrzymującego masę $m$ statyczne ugięcie pod ciężarem wynosi

$$\delta_s=\frac{mg}{k}.$$

Po połączeniu tej relacji z $\omega_n=\sqrt{k/m}$ otrzymuje się

$$\omega_n=\sqrt{\frac{g}{\delta_s}},\qquad f_n=\frac{1}{2\pi}\sqrt{\frac{g}{\delta_s}}.$$

Relacja ta obowiązuje dla jednego pionowego stopnia swobody, liniowej sprężyny i tej samej sztywności użytej w równowadze statycznej oraz małych drganiach. Nie można przenosić jej bez korekty na podporę poziomą, układ wstępnie napięty, izolator o charakterystyce progresywnej, mechanizm z przełożeniem geometrycznym ani materiał, którego sztywność dynamiczna różni się od siecznej sztywności statycznej.

Ugięcie statyczne ustala położenie robocze i dzieli dostępny skok na oba kierunki. Do amplitudy $Z$ dochodzą tolerancje wysokości, pełzanie, rozszerzalność cieplna, błąd montażu i możliwe przemieszczenie równowagi pod innym obciążeniem. Prześwit powinien być oceniany osobno w każdym kierunku, ponieważ nominalnie symetryczny ogranicznik może mieć po obciążeniu bardzo niesymetryczną rezerwę.

Gdy $|z|$ osiąga dostępny prześwit $s$, równanie liniowe traci ważność. Najprostszy model kontaktu dodaje siłę równą zeru dla $|z|<s$, a po przekroczeniu granicy wprowadza sztywność kontaktową $k_o$, często znacznie większą od $k$, oraz tłumienie kontaktowe $c_o$. Powstaje odcinkowo nieliniowa droga obciążenia. Zmienia się częstotliwość chwilowa, pojawiają się wyższe składowe częstotliwości, impuls siły i możliwość odbicia.

Ogranicznik nie jest w takim stanie bierną granicą na rysunku. Jego prześwit, podatność, powierzchnia styku, tłumienie i trwałość współtworzą odpowiedź. Do oceny chwili pierwszego kontaktu potrzebny jest przebieg $z(t)$, a do oceny skutków — model oraz pomiar siły lub wiarygodna detekcja kontaktu. Samo widmo liniowej transmitancji nie przewiduje odpowiedzi po zwarciu miękkiej drogi przez znacznie sztywniejszy styk.

Udar porządkuje stosunek czasu impulsu do okresu własnego

W odpowiedzi przejściowej samo maksimum przyspieszenia wejścia nie wystarcza. Istotny jest bezwymiarowy stosunek czasu trwania impulsu $\tau$ do okresu własnego $T_n=1/f_n$, jego kształt, znak i warunki początkowe. Bardzo krótki impuls w porównaniu z $T_n$ przede wszystkim nadaje masie zmianę prędkości; największy ruch względny może wystąpić już po ustaniu wymuszenia. Obciążenie znacznie dłuższe od $T_n$ zbliża odpowiedź do ruchu quasi-statycznego, o ile nie ma gwałtownego zbocza końcowego. Gdy $\tau$ jest porównywalne z $T_n$, faza wymuszenia i ruchu może prowadzić do szczególnie dużej odpowiedzi.

Tłumienie zwykle skraca wybrzmiewanie i ogranicza odpowiedź rezonansową, lecz przy bardzo krótkim impulsie może mieć mały wpływ na pierwszy szczyt. Zmiękczenie zawieszenia może obniżyć przyspieszenie masy, a jednocześnie zwiększyć ruch względny i ryzyko kontaktu. Dlatego udar należy całkować w dziedzinie czasu dla zmierzonego lub uzasadnionego przebiegu, z kontrolą $x(t)$, $z(t)$, reakcji i ewentualnego ogranicznika. Widmo odpowiedzi udarowej (shock response spectrum, SRS) porównuje rodzinę oscylatorów odniesienia, ale nie zachowuje fazy ani nie rozstrzyga samodzielnie, kiedy w konkretnym wyrobie zamknął się luz.8

5. Dydaktyczny przykład obliczeniowy

Rozważmy całkowicie umowny zespół laboratoryjny o masie $m=12 \mathrm{kg}$. Podpory mają utworzyć pionowy model SDOF o częstotliwości własnej $f_n=8 \mathrm{Hz}$ i względnym tłumieniu $\zeta=0{,}12$. Podstawa porusza się harmonicznie z częstotliwością $f=24 \mathrm{Hz}$ i amplitudą przemieszczenia $Y=0{,}40 \mathrm{mm}$. Dane nie opisują żadnego rzeczywistego uzbrojenia.

Najpierw oblicza się

$$\omega_n=2\pi f_n=50{,}3 \mathrm{rad/s},$$

a następnie parametry równoważnego izolatora:

$$k=m\omega_n^2 =12\cdot 50{,}3^2 \approx 30{,}3 \mathrm{kN/m},$$

$$c=2\zeta m\omega_n =2\cdot0{,}12\cdot12\cdot50{,}3 \approx145 \mathrm{N\,s/m}.$$

Stosunek częstotliwości wynosi $r=24/8=3$. Wspólny mianownik transmitancji ma wartość

$$D=\sqrt{(1-r^2)^2+(2\zeta r)^2} =\sqrt{64+0{,}5184} \approx8{,}03.$$

Transmitancja ruchu bezwzględnego jest równa

$$T_x=\frac{\sqrt{1+(2\zeta r)^2}}{D} =\frac{1{,}232}{8{,}03} \approx0{,}153.$$

Amplituda ruchu masy wyniesie zatem

$$X=T_xY\approx0{,}061 \mathrm{mm}.$$

Amplituda przyspieszenia podstawy to

$$A_y=\omega^2Y =(2\pi\cdot24)^2\cdot0{,}00040 \approx9{,}10 \mathrm{m/s^2},$$

a przewidywana amplituda przyspieszenia masy

$$A_x=T_xA_y\approx1{,}40 \mathrm{m/s^2}.$$

Jeżeli wymaganie przewiduje $T_x\leq0{,}20$ i $A_x\leq1{,}8 \mathrm{m/s^2}$, oba kryteria ruchu bezwzględnego są spełnione. Siła przekazywana ma amplitudę około

$$|F_p|=mA_x\approx16{,}8 \mathrm{N},$$

lecz jej dopuszczalność wymaga osobnego limitu interfejsu.

Kontrola skoku prowadzi do innego wyniku:

$$T_z=\frac{r^2}{D} =\frac{9}{8{,}03} \approx1{,}12,$$

$$Z=T_zY\approx0{,}448 \mathrm{mm}.$$

Załóżmy, że najmniejszy zmierzony prześwit w analizowanym kierunku wynosi $0{,}50 \mathrm{mm}$, a wymagany zapas na tolerancje i nieujęte zaburzenia to $0{,}10 \mathrm{mm}$. Dopuszczalna amplituda względna jest wtedy równa $0{,}40 \mathrm{mm}$. Obliczone $0{,}448 \mathrm{mm}$ przekracza tę wartość, chociaż izolacja przyspieszenia wygląda korzystnie.

Decyzja dla podanych danych brzmi: wariant nie może zostać zatwierdzony. Spełnia wymaganie transmitancji ruchu, ale nie zachowuje wymaganego zapasu prześwitu. Różnica $0{,}052 \mathrm{mm}$ do fizycznego ogranicznika również nie stanowi wiarygodnego marginesu, a obliczenie harmoniczne nie obejmuje rozruchu, udaru ani rozrzutu sztywności. Dalsza zmiana wymaga ponownego policzenia obu stron kompromisu. Zwiększenie prześwitu może usunąć ryzyko kontaktu bez zmiany $T_x$. Zmiana sztywności może natomiast zwiększyć albo zmniejszyć skok, zależnie od położenia względem rezonansu; w tym przykładzie jej zwiększenie obniżyłoby $r$, zbliżyło układ do rezonansu i zwiększyło $T_z$, zamiast automatycznie ograniczyć ruch.

6. Studium przypadku: miękkie zawieszenie dobija do ogranicznika

Wymaganie i anomalia.

Umowny zespół badawczy miał zachować $T_x\leq0{,}25$ w paśmie 30–80 Hz. Wymaganie udarowe określono funkcjonalnie: brak kontaktu z ogranicznikiem i co najmniej $0{,}20 \mathrm{mm}$ gwarantowanej rezerwy po uwzględnieniu tolerancji geometrycznej oraz niepewności pomiaru ruchu. Dla nominalnego prześwitu $1{,}00 \mathrm{mm}$, rozszerzonej niepewności geometrii $U_s=0{,}02 \mathrm{mm}$ i rozszerzonej niepewności toru ruchu $U_z=0{,}03 \mathrm{mm}$ dawało to projektowy limit wskazania $|z|_{\max}\leq1{,}00-0{,}02-0{,}20-0{,}03=0{,}75 \mathrm{mm}$. Wartość ta była pochodną geometrii i budżetu niepewności, a nie niezmiennym ograniczeniem funkcji. Próba wibracyjna spełniła kryterium pasmowe. Podczas udaru zarejestrowano jednak krótki pik wysokoczęstotliwościowy, a po oględzinach zauważono ślad w pobliżu ogranicznika.

Konkurencyjne hipotezy i sygnatury.

Rozpatrzono pięć konkurencyjnych hipotez. Kontakt z ogranicznikiem powinien dać osiągnięcie prześwitu przez $z(t)$, równoczesny impuls siły lub sygnał styku oraz szerokopasmową odpowiedź pojawiającą się dopiero od określonego poziomu wejścia. Nieliniowość izolatora bez kontaktu powinna ujawnić zależne od amplitudy przesunięcie rezonansu, zniekształcenie pętli siła–przemieszczenie i różnicę między przemiataniem w górę oraz w dół, bez sygnału ogranicznika. Mostek przez przewód powinien korelować pik z naprężeniem wiązki i zmieniać odpowiedź po kontrolowanej zmianie jej prowadzenia. Rezonans przyrządu mocującego powinien być widoczny w jego własnym kanale oraz w lokalnej funkcji odpowiedzi częstotliwościowej (frequency response function, FRF), bez jednoczesnego osiągnięcia granicy ruchu względnego. Nasycenie czujnika powinno pozostawić spłaszczony przebieg na granicy zakresu, nietypowy czas powrotu albo rozbieżność z równoległym czujnikiem o większym zakresie.

Eksperyment rozstrzygający.

Układ pomiarowy objął przyspieszenie podstawy i chronionej masy, bezkontaktowy pomiar $z(t)$, podkładkę siłową przy ograniczniku, prosty obwód detekcji styku, odkształcenie uchwytu przewodu oraz przyspieszenie przyrządu mocującego. Na masie zamontowano dwa zsynchronizowane akcelerometry o różnych zakresach. Kontrolowane zmiany wykonywano pojedynczo i przy obniżonym poziomie wejścia: zmieniano prowadzenie wiązki bez naruszania jej funkcji, zastępowano lokalne mocowanie czujnika zweryfikowanym mocowaniem oraz przeprowadzano stopniowanie amplitudy udaru. Stanowisko osłonięto przed niekontrolowanym ruchem i możliwym odrzutem elementów.

Wynik i przyczyna.

Pomiar geometryczny przed próbą dał rzeczywisty prześwit $s=0{,}99 \mathrm{mm}$ z rozszerzoną niepewnością $U_s=0{,}02 \mathrm{mm}$, a rozszerzona niepewność toru przemieszczenia wynosiła $U_z=0{,}03 \mathrm{mm}$. Względne przemieszczenie dochodziło do $0{,}98 \mathrm{mm}$. W tym samym przedziale czasowym pojawiały się sygnał styku, impuls na podkładce siłowej i początek pasma wysokoczęstotliwościowego. Kontakt przy wskazaniu zbliżonym do dolnej granicy rzeczywistego prześwitu był zgodny z budżetem niepewności. Kanał przyrządu mocującego nie wykazywał własnego piku poprzedzającego zdarzenie, a zastosowanie zweryfikowanego mocowania czujnika nie przesunęło piku ani nie zmieniło jego fazy. Zmiana prowadzenia przewodu nie zmieniła chwili ani poziomu wystąpienia impulsu, a pomiar jego uchwytu nie pokazał wzrostu siły przed kontaktem. Czujnik o niższym zakresie rzeczywiście się nasycał, lecz kanał o wyższym zakresie potwierdził istnienie odpowiedzi po zderzeniu. Nasycenie zafałszowało wartość szczytu, ale nie spowodowało zjawiska. Stopniowanie wejścia dawało odpowiedź zbliżoną do liniowej aż do osiągnięcia około $0{,}95 \mathrm{mm}$; potem pojawiała się gwałtowna zmiana. Zestaw sygnatur wskazał kontakt z ogranicznikiem jako przyczynę pierwotną.

Korekta.

Z modelu czasowego bez ogranicznika, skorygowanego wynikami niższych poziomów próby, otrzymano obwiednię $|z|_{\max}=1{,}08 \mathrm{mm}$. Zespół odpowiedzialny za wymagania potwierdził, że stanem granicznym jest kontakt i utrata rezerwy, a nie sama wartość $0{,}75 \mathrm{mm}$ wyprowadzona dla starej geometrii. Jedyną wprowadzoną korektą było więc zwiększenie użytecznego prześwitu do nominalnych $1{,}35 \mathrm{mm}$, bez zmiany głównych izolatorów, tłumienia i prowadzenia przewodów. Po wykonaniu zmierzono $s=1{,}35 \mathrm{mm}$ z rozszerzoną niepewnością $U_s=0{,}02 \mathrm{mm}$. Dolna granica prześwitu wynosiła $1{,}33 \mathrm{mm}$; po odjęciu wymaganej rezerwy $0{,}20 \mathrm{mm}$ dopuszczalna górna granica ruchu wynosiła $1{,}13 \mathrm{mm}$. Zmiana nie była odstępstwem od funkcji, lecz kontrolowaną zmianą geometrii, dla której ponownie wyprowadzono kryterium.

Regresja i reguła decyzji.

W regresji udarowej największy zarejestrowany ruch względny wyniósł $1{,}06 \mathrm{mm}$. Po dodaniu rozszerzonej niepewności toru jego konserwatywna górna granica wynosiła $1{,}09 \mathrm{mm}$, czyli o $0{,}04 \mathrm{mm}$ mniej od dopuszczalnych $1{,}13 \mathrm{mm}$. Detektor styku i podkładka siłowa nie wykazały kontaktu, oba akcelerometry pracowały w zakresie, a impuls wysokoczęstotliwościowy zniknął. Powtórzona próba 30–80 Hz potwierdziła $T_x\leq0{,}25$, brak nowego maksimum oraz zgodność fazy i amplitudy z wynikiem sprzed korekty.

Reguła decyzji jest jednoznaczna: rozwiązanie można przyjąć tylko wtedy, gdy we wszystkich ocenianych konfiguracjach zachodzi $|z|_{\max}+U_z\leq s-U_s-\Delta s$, gdzie $U_z$ i $U_s$ są rozszerzonymi niepewnościami ruchu oraz prześwitu, żaden niezależny kanał nie wskazuje kontaktu, czujniki pozostają poza nasyceniem, a pierwotne wymaganie transmitancji nadal jest spełnione. Niespełnienie choć jednego z tych warunków oznacza odrzucenie wariantu albo ponowne otwarcie analizy.

Tak zamknięty model SDOF wyjaśnia dominujący ruch i drogę kontaktową. Gdy odpowiedzi poszczególnych punktów zaczynają się różnić, pojawiają się obroty, nierówny docisk albo kilka ruchomych punktów podstawy, te same wielkości trzeba zapisać już w modelu wielu stopni swobody.

7. Od bryły sztywnej do modelu 6-DOF

Jeżeli najniższe mody zespołu odpowiadają ruchowi niemal nieodkształcalnej obudowy na podatnych podporach, jego stan opisuje wektor $\mathbf q=[\mathbf u^T\;\boldsymbol\theta^T]^T$, zawierający trzy translacje i trzy małe obroty. Równanie ruchu ma postać

$$\mathbf M\ddot{\mathbf q}+\mathbf C\dot{\mathbf q}+\mathbf K\mathbf q=\mathbf f(t).$$

Przy początku układu w środku masy macierz $\mathbf M$ zawiera bloki $m\mathbf I$ oraz tensor bezwładności $\mathbf J$. Poza osiami głównymi jego elementami są również produkty bezwładności, więc przyspieszenia kątowe wokół osi geometrycznych mogą być sprzężone. Przejście do osi głównych diagonalizuje $\mathbf J$, ale nie musi diagonalizować $\mathbf K$ ani $\mathbf C$: sprzężenie powodowane geometrią i anizotropią podpór pozostaje fizyczne.

Charakterystyki podpór składa się dopiero po wyrażeniu ich we wspólnych współrzędnych. Dla podpory $i$, o orientacji $\mathbf R_i$, jest $\mathbf K_i^{(g)}=\mathbf R_i\mathbf K_i^{(l)}\mathbf R_i^T$, analogicznie dla tłumienia. Jej ruch względny wynosi $\boldsymbol\delta_i=\mathbf B_i\mathbf q-\mathbf u_{b,i}$, gdzie $\mathbf B_i=[\mathbf I\;-[\mathbf r_i]_\times]$, a $\mathbf r_i$ prowadzi od środka masy do podpory. Z zasady prac wirtualnych:17

$$\mathbf K=\sum_i\mathbf B_i^T\mathbf K_i^{(g)}\mathbf B_i,\qquad \mathbf C=\sum_i\mathbf B_i^T\mathbf C_i^{(g)}\mathbf B_i.$$

Gdy punkty podstawy poruszają się niezależnie, ogólny wektor $\mathbf f(t)$ trzeba rozpisać. Dla liniowych podpór translacyjnych równanie przyjmuje postać

$$\mathbf M\ddot{\mathbf q}+\mathbf C\dot{\mathbf q}+\mathbf K\mathbf q =\mathbf f_{\mathrm{ext}}+ \sum_i\mathbf B_i^T\left(\mathbf K_i^{(g)}\mathbf u_{b,i} +\mathbf C_i^{(g)}\dot{\mathbf u}_{b,i}\right).$$

Zapis ten wprowadza momenty od reakcji, sprzężenia translacja–obrót i różne charakterystyki wynikające z docisku, napięcia wstępnego oraz tolerancji. Zakłada jednak, że podpora przenosi siły opisane macierzą translacyjną 3×3. Jeżeli ma własną sztywność momentową albo sprzężenie siła–moment, jej opis trzeba rozszerzyć do relacji 6×6 i odpowiednio rozbudować macierz transformacji. Środek sztywności jest jednoznaczny tylko w układach szczególnych, na przykład w płaskim podparciu równoległymi sprężynami skalarnymi. W przestrzennym układzie anizotropowym, zwłaszcza ze sztywnością zależną od częstotliwości, jeden taki punkt może nie istnieć; o kołysaniu rozstrzygają wtedy pozadiagonalne człony macierzy i postacie własne pary $(\mathbf K,\mathbf M)$.

Do podatnej podstawy lub modelu metodą elementów skończonych (MES) przechodzi się, gdy punkty podstawy nie poruszają się zgodnie co do amplitudy i fazy, jej podatność jest porównywalna z podatnością podpór, w paśmie decyzji leży mod obudowy lub ramy albo kilka funkcji odpowiedzi częstotliwościowej (frequency response function, FRF) zależnych od miejsca nie daje się odtworzyć modelem bryłowym. Przesłanką są też istotne odkształcenia między punktami uznanymi wcześniej za sztywne. Podpory łączą wówczas stopnie swobody zespołu z lokalnymi węzłami podstawy, a model 6-DOF pozostaje kontrolą rzędu wielkości i narzędziem identyfikacji brakującego sprzężenia.

8. Zespół na czterech podporach: rozstrzygnięcie przyczyny kołysania

Wymaganie i anomalia. Wszystkie zapisy $a\pm U$ w tym studium oznaczają wynik $a$ oraz przypisaną mu rozszerzoną niepewność $U$ dla przyjętego w budżecie współczynnika rozszerzenia; reguła decyzji używa odpowiedniej górnej albo dolnej granicy, a nie samej wartości centralnej.

Rozważmy dydaktyczny zespół o masie 48 kg na czterech podporach w narożach prostokąta 0,80 × 0,60 m. Środek masy jest przesunięty o 70 mm i leży 210 mm nad płaszczyzną podparcia. Przed próbą ustalono limity: obrót 2,0 mrad, skok naroża 2,5 mm, brak odciążenia podpory i siła w złączu wiązki poniżej 100 N. Podczas przemiatania przy 31 Hz przeciwległe naroża poruszają się niemal w przeciwfazie; obrót osiąga 3,1 mrad, skok 2,9 mm, a siła wiązki 146 N.

Hipotezy i przewidywane sygnatury.

Geometria masy była znanym składnikiem bazowym H0, a nie konkurencyjną przyczyną: powinna wyznaczyć część kołysania i jego oś już przy nominalnych podporach. Rozpatrzono następnie cztery źródła brakującej odpowiedzi. H1, rozrzut właściwości elementu, powinien podążyć za zamienioną podporą. H2, pozycyjnie zależne napięcie wstępne lub osadzenie, powinno pozostać związane z miejscem, lecz zmienić się po kontrolowanym ponownym montażu i wyrównaniu docisków. H3, podatność podstawy, wymaga różnych amplitud lub faz przy czterech mocowaniach. H4, równoległa droga przez wiązkę, przewiduje zgodność fazową jej siły z brakującą częścią bilansu reakcji oraz zmianę modu po kontrolowanej zmianie wyłącznie tej drogi.

Eliminacja H1–H3.

Bilans statyczny potwierdził wpływ położenia środka masy, ale model 6-DOF przewidywał tylko 1,5 mrad. Zmierzone sztywności podpór przy rzeczywistych dociskach różniły się o 7%; ich skrajne przypisania zmieniały obrót najwyżej o 0,2 mrad, a zamiana elementów nie przeniosła asymetrii, więc H1 nie była dominująca. Aby ograniczyć H2, wykonano powtórny montaż według tej samej procedury, zmierzono reakcje statyczne i skorygowano dociski do tolerancji 3%. Oś kołysania pozostała związana z konstrukcją, a obrót zmienił się o 0,1 mrad, czyli mniej niż brakujące około 1,6 mrad. Wynik ten nie falsyfikuje H2: zależna od pozycji podatność styku lub osadzenia może odtworzyć się po identycznej procedurze. Cztery czujniki podstawy były zgodne amplitudowo w granicach 3% i fazowo w granicach 2°. Zamiast uznać te różnice za pomijalne, wprowadzono cztery zmierzone zespolone ruchy $\mathbf u_{b,i}$ do prawej strony równania. Analiza wrażliwości dała dodatkowy obrót $0{,}12\pm0{,}05 \mathrm{mrad}$, a skrajna kombinacja dopuszczona przez niepewność pomiaru — $0{,}19 \mathrm{mrad}$. H3 mogła więc wyjaśnić niewielką część odpowiedzi. H1–H3 pozostały składnikami obwiedni, lecz nie wyjaśniały dominującego niedoboru modelu.

Potwierdzenie H4 i kontrolowana korekta.

Po zmierzeniu zespolonej sztywności wiązki w reprezentatywnym ułożeniu model przewidział 30,6 Hz, właściwą oś kołysania i 3,0 mrad. Przy 31 Hz zespolony bilans reakcji podpór i bezwładności pozostawiał moment $52\pm6 \mathrm{N\,m}$; moment wyznaczony niezależnie z siły wiązki i jej ramienia wynosił $49\pm7 \mathrm{N\,m}$, z różnicą fazy 8°. Eksperyment rozstrzygający nie usuwał wiązki, lecz wprowadzał pętlę odciążającą o znanej, mniejszej sztywności, bez zmian masy, podpór i podstawy. Rezonans przesunął się do 27,8 Hz, obrót spadł do 1,4 mrad, siła w złączu do 41 N, a niewyjaśniony moment do $11\pm6 \mathrm{N\,m}$. Zmiana była zgodna z H4 jako dominującym brakującym sprzężeniem, podczas gdy H0 nadal wyznaczała bazową odpowiedź około 1,5 mrad, a H1–H3 pozostawały w budżecie niepewności.

Regresja i reguła decyzji.

Korektą było trwałe przeprowadzenie wiązki przez taką pętlę i przesunięcie jej podparcia ku płaszczyźnie izolatorów; sztywności izolatorów nie zmieniano. Regresja objęła skrajne położenia środka masy, dociski i zmierzone granice sztywności podpór, zmierzoną obwiednię ruchów podstawy oraz tolerancję ułożenia wiązki. Powtórzono pomiar sześciu ruchów bryły, czterech ruchów podstawy, przemieszczeń względnych naroży, reakcji i siły wiązki. Najgorszy wariant dał obrót $1{,}80\pm0{,}10 \mathrm{mrad}$, maksymalne przemieszczenie względne naroża $2{,}10\pm0{,}10 \mathrm{mm}$ i siłę $73\pm6 \mathrm{N}$. Najmniejsza chwilowa reakcja podpory wynosiła $18\pm6 \mathrm{N}$, zatem jej konserwatywna dolna granica 12 N pozostawała dodatnia. Górne granice pozostałych miar, odpowiednio 1,90 mrad, 2,20 mm i 79 N, były niższe od limitów 2,0 mrad, 2,5 mm i 100 N. Przed regresją ustalono tolerancje korelacji: 5% dla częstotliwości, 10° dla fazy i 15% dla amplitudy odpowiedzi. Największe rozbieżności model–pomiar wyniosły odpowiednio 2,1%, 6° i 9%. Reguła decyzji była koniunkcyjna: każdy limit wraz z niepewnością musiał być spełniony w całej obwiedni, korelacja zmieścić się w podanych tolerancjach, a kontrolowana zmiana drogi — wywołać przewidzianą sygnaturę.

9. Równoległe drogi jako sieć impedancji dynamicznych

Podpory i pozostałe połączenia tworzą sieć dróg dynamicznych. Drogi równoległe łączące te same dwa węzły podlegają temu samemu ruchowi względnemu, lecz w modelu 6-DOF poszczególne punkty mają na ogół różne lokalne ruchy $\boldsymbol\delta_i$. Dla drogi liniowej użyteczna jest zespolona sztywność dynamiczna $K_d(\omega)=F/X$ albo impedancja $\mathcal Z(\omega)=F/V$. Reakcja całkowita jest sumą zespoloną $F_\Sigma=\sum_jF_j$, z zachowaniem fazy; suma modułów ukrywa znoszenie. Przy fazorach RMS średnią moc czynną drogi ocenia się jako $\operatorname{Re}\{F_jV_j^*\}$; czynnik $1/2$ pojawiłby się przy fazorach z amplitudami szczytowymi. Duża reakcja elementu zachowawczego może współistnieć z zerową średnią mocą czynną, ponieważ energia jest oddawana w każdym cyklu; osobno ocenia się więc amplitudę reakcji, moc czynną i energię wymienianą lub rozpraszaną w cyklu.

Wiązka lub rura staje się istotna, gdy jej kierunkowa impedancja w rzeczywistym ułożeniu jest porównywalna z impedancją podpór. Zapas długości może pracować miękko na zginanie, a po wyprostowaniu — sztywno osiowo; temperatura i ciśnienie zmieniają punkt pracy. Taśma pełniąca funkcję cieplną albo połączenia wyrównawczego może wywierać znaczny moment przez długie ramię mimo małej siły. Ogranicznik przed kontaktem ma udział bliski zeru, po zamknięciu luzu tworzy natomiast nieliniową drogę wysokiej impedancji, niewidoczną w małoamplitudowej FRF.

Udział określa bilans zespolonych reakcji lub mocy przy wszystkich interfejsach. Identyfikację przyczynową daje kontrolowana zmiana jednej impedancji — na przykład promienia pętli, miejsca podparcia albo prześwitu — przy niezmienionych innych warunkach. Zgodna z modelem zmiana odpowiedzi potwierdza drogę; brak równowagi między reakcjami a bezwładnością wskazuje brakującą drogę, błąd fazy lub niepełny model.

10. Absorber jako układ dwumasowy

Dla masy głównej $m_1$, połączonej z otoczeniem przez $k_1,c_1$, i absorbera $m_2$ dołączonego przez $k_2,c_2$, równania przy wymuszeniu $F(t)$ mają postać:16

$$\begin{bmatrix}m_1&0\\0&m_2\end{bmatrix}\ddot{\mathbf x}+ \begin{bmatrix}c_1+c_2&-c_2\\-c_2&c_2\end{bmatrix}\dot{\mathbf x}+ \begin{bmatrix}k_1+k_2&-k_2\\-k_2&k_2\end{bmatrix}\mathbf x= \begin{bmatrix}F\\0\end{bmatrix}.$$

Częstotliwość odłączonego absorbera wynosi $\omega_a=\sqrt{k_2/m_2}$, a $\mu=m_2/m_{1,\mathrm{eff}}$ odnosi jego masę do modalnej masy efektywnej kontrolowanego modu. Po sprzęgnięciu powstają dwa mody. W podatności $H_{x_1F}$ między biegunami może pojawić się antyrezonans, lecz inaczej wygląda on w akcelerancji, ruchu absorbera, sile interfejsowej lub odpowiedzi poza obszarem aktywnym modu. „Rozszczepienie piku” ma więc sens dopiero po wskazaniu wejścia, odpowiedzi i punktu pomiaru.

Większe $\mu$ silniej rozdziela częstotliwości i zwykle zmniejsza wymagany ruch względny, ale zwiększa masę oraz reakcję mocowania. Tłumienie $c_2$ wygładza maksima i ogranicza wrażliwość na odstrojenie, zarazem spłyca nominalne minimum. Za małe daje wąski antyrezonans i duże $|x_2-x_1|$, za duże osłabia działanie strojone. Dla małych zmian $\Delta\omega_a/\omega_a\approx\tfrac12(\Delta k_2/k_2-\Delta m_2/m_2)$, lecz odporność ocenia się z pełnej FRF, wraz z dryfem konstrukcji, tłumieniem i limitem skoku. Także ogólne opracowania pasywnej stabilizacji podkreślają, że absorber trzeba analizować jako element sprzężony z konstrukcją, a nie niezależną część o jednej częstotliwości katalogowej.15

11. Absorber po zmianie konfiguracji: przesiewowe ćwiczenie modelowe

Zakres i kryteria.

Wszystkie liczby w tym przykładzie są dydaktyczne i pochodzą z jawnego modelu liniowego z sekcji 10; nie są wynikami prób rzeczywistego urządzenia. Ćwiczenie wykonuje przesiewowe sprawdzenie skończonego zbioru wariantów i nie stanowi dowodu kwalifikacyjnego ani ciągłej analizy całej przestrzeni tolerancji. Przyjęto pasmo decyzji 68–90 Hz i limit akcelerancji $|H_{a_1F}|\leq0{,}70 \mathrm{(m/s^2)/N}$. Dwa lokalne maksima potrzebne do oceny strojenia oraz amplitudę ruchu względnego wyszukuje się w szerszym przedziale 60–100 Hz; stosunek wyższego do niższego maksimum nie może przekroczyć 1,25, a $|x_2-x_1|$ przy harmonicznej sile o amplitudzie 250 N — 3,0 mm. Pierwszy limit chroni odpowiedź panelu, drugi wykrywa utratę zrównoważonego strojenia, a trzeci chroni zakres ruchu absorbera. Kryteria są umowne, lecz przed obliczeniem jednoznacznie wiążą trzy różne wielkości z decyzją.

Punkt bazowy.

Modalną masę efektywną przyjęto jako $m_1=10{,}0 \mathrm{kg}$, tłumienie konstrukcji jako $\zeta_1=0{,}015$, a masę absorbera jako $m_2=0{,}50 \mathrm{kg}$, zatem $\mu=0{,}050$. Przed zmianą konfiguracji mod konstrukcji leżał przy $f_1=73{,}8 \mathrm{Hz}$, częstotliwość odłączonego absorbera wynosiła $f_a=73{,}6 \mathrm{Hz}$, a $\zeta_2=0{,}12$. Odpowiada to $k_1\approx2{,}15 \mathrm{MN/m}$, $c_1\approx139 \mathrm{N\,s/m}$, $k_2\approx107 \mathrm{kN/m}$ i $c_2\approx56 \mathrm{N\,s/m}$. Model daje dwa maksima około 0,60 przy 67,1 Hz oraz 0,50 przy 81,6 Hz, minimum około 0,41 przy 74,2 Hz i największy ruch względny około 2,6 mm. Pierwsze maksimum leży nieznacznie poniżej dolnej granicy pasma decyzji, ale uwzględniono je przy kontroli kształtu FRF i skoku. Konfiguracja początkowa spełnia wszystkie trzy kryteria: w paśmie nie przekracza 0,70, stosunek maksimów wynosi 1,19, a skok pozostaje poniżej 3,0 mm.

Hipoteza odstrojenia.

Scenariusz modelowy zakłada, że po zmianie konstrukcji jej mod przesunął się do $f_1=77{,}8 \mathrm{Hz}$, podczas gdy parametry absorbera pozostały bez zmian. Jest to założenie H1, nie wniosek z fikcyjnego pomiaru. W realnym badaniu należałoby przeciwstawić mu H2 — zmianę tłumienia, H3 — luz lub uszkodzenie mocowania oraz H4 — artefakt pomiarowy. H1 przewiduje nierówne maksima po obu stronach antyrezonansu i przesunięcie modu konstrukcji przy czasowo zablokowanym absorberze. H2 zmieniałaby głównie szerokość pików i zanik swobodny; H3 dawałaby zależność od poziomu, harmoniczne albo histerezę; H4 nie powinien odtwarzać się w niezależnym torze przemieszczenia. Rozstrzygnięcie wymagałoby skalibrowanej siły, zsynchronizowanych odpowiedzi obu mas, fazy, ruchu względnego, stopniowania poziomu, zamiany czujników i przemiatania w obu kierunkach. Ponieważ takich danych tu nie ma, ćwiczenie ocenia wyłącznie następstwa H1, a nie dowodzi jej prawdziwości.

Dla wariantu odstrojenia, bez jakiegokolwiek dopasowywania do fikcyjnych danych, model daje maksima około 0,43 przy 69,5 Hz i 0,65 przy 83,7 Hz, minimum około 0,40 przy 73,7 Hz oraz skok około 2,0 mm. Limity bezwzględne pozostają spełnione, ale stosunek maksimów rośnie do 1,51 i przekracza 1,25. Odrzucenie wynika więc z utraty odpornego, zrównoważonego strojenia, a nie z arbitralnie dobranego przekroczenia akcelerancji.

Korekta modelowa.

Korektę wyznaczono przez przeszukanie $f_a=70$–82 Hz co 0,005 Hz. Krok jest rozdzielczością algorytmu, nie deklarowaną dokładnością fizyczną. Dla każdej wartości obliczono maksimum akcelerancji w paśmie 68–90 Hz na siatce częstotliwościowej co 0,0025 Hz, przy stałych $m_2=0{,}50 \mathrm{kg}$ i $\zeta_2=0{,}12$. Minimum funkcji celu leży przy $f_a\approx76{,}4 \mathrm{Hz}$, czyli $k_2\approx115 \mathrm{kN/m}$ oraz $c_2\approx58 \mathrm{N\,s/m}$. Modelowa korekta jest więc równoczesnym dostrojeniem sztywności i współczynnika tłumienia w celu utrzymania $\zeta_2=0{,}12$, a nie zmianą samej sprężyny. Sztywność rośnie o około 7,7%, a nie o $(77{,}8/73{,}6)^2-1\approx11{,}7\%$: równość częstotliwości odłączonego absorbera i konstrukcji nie jest optimum sprzężonego układu dwumasowego.

Po korekcie oba maksima wynoszą około 0,54 i występują w pobliżu 70,3 oraz 85,2 Hz. Minimum ma wartość około 0,42 przy 76,9 Hz, a amplituda ruchu względnego przy sile 250 N wynosi około 2,2 mm. Stosunek maksimów jest praktycznie równy jedności. Pełną precyzję zachowano w obliczeniu rozwiązującym zespoloną macierz sztywności dynamicznej z sekcji 10; nie przenosi się jej na zapis wyniku inżynierskiego.

Przesiew tolerancji i granica wniosku.

Sprawdzenie nie ogranicza się do punktu nominalnego. Siatka 243 punktów brzegowych i nominalnych obejmuje $m_1=9{,}9;10{,}0;10{,}1 \mathrm{kg}$, $m_2=0{,}495;0{,}500;0{,}505 \mathrm{kg}$, $f_1=77{,}4;77{,}8;78{,}2 \mathrm{Hz}$, $f_a$ równe 99,5%, 100% i 100,5% wartości nominalnej oraz $\zeta_2=0{,}11;0{,}12;0{,}13$. Spośród tych kombinacji największa akcelerancja w paśmie wynosi około 0,60 $\mathrm{(m/s^2)/N}$, największa amplituda ruchu względnego przy 250 N — około 2,5 mm, a największy stosunek dwóch maksimów — 1,15. Po dodaniu umownych rozszerzonych niepewności 0,03 $\mathrm{(m/s^2)/N}$, 0,15 mm oraz osobno wyznaczonej $U_R=0{,}04$ konserwatywne granice wynoszą 0,63, 2,65 mm i 1,19; pozostają poniżej limitów w zbadanym zbiorze. Niepewność ilorazu ocenia się osobno, ponieważ wspólne błędy skali maksimów są skorelowane.

Wniosek jest pozytywny wyłącznie dla 243 obliczonych kombinacji. Zakresów $f_1$, $f_a$ i $\zeta_2$ nie przypisano temperaturom, więc nie jest to regresja temperaturowa, a wyniku nie wolno interpolować na wnętrze przedziałów. Decyzja projektowa pozostaje otwarta do czasu zbadania wnętrza przestrzeni tolerancji, rozstrzygnięcia H1–H4 i doświadczalnego potwierdzenia liniowości, rzeczywistej masy modalnej oraz własności mocowania. Wniosek nie obejmuje innego stosunku mas, tłumienia, punktu siły ani miejsca obserwacji.

12. Wybór elementu jako konsekwencja mechanizmu

Punktem wyjścia jest potrzebna zespolona sztywność dynamiczna $K^*(\omega,T,A,N)=K'+jK''$, określona w paśmie częstotliwości, temperaturze $T$, amplitudzie $A$ i po liczbie cykli $N$. Część zachowawcza ustala częstotliwości i skok, stratna zaś kształtuje rezonans oraz nagrzewanie. Element musi jednocześnie zapewnić nośność, dopuszczalny ruch, kontrolę reakcji i stabilność w środowisku. ISO 10846-2 opisuje bezpośredni pomiar translacyjnej sztywności dynamicznej podatnych podpór przy zadanym napięciu wstępnym; zakres tej normy nie obejmuje automatycznie całego zespołu ani wszystkich kierunków i nieliniowości.10

Elastomer łączy sprężystość z rozpraszaniem, lecz jego moduł i współczynnik strat zależą od temperatury, częstotliwości, amplitudy i wstępnego ścisku. Badania NASA-CR-134939 pokazują wprost potrzebę redukcji danych sztywności i tłumienia do funkcji temperatury, częstotliwości, geometrii i poziomu dyssypacji.18 Pełzanie przesuwa równowagę oraz dociski podpór; w próżni znaczenie mają odgazowanie i zanieczyszczenie sąsiednich powierzchni. Potrzebna jest obwiednia charakterystyk po kondycjonowaniu i starzeniu, nie jedna krzywa pokojowa.

Element metalowo-linowy rozprasza przez mikropoślizgi i tarcie drutów. Badania modeli i pętli histerezy takich izolatorów wykazały istotność tarcia Coulomba oraz zależność odpowiedzi od geometrii i warunków wzbudzenia.19 Element nie wymaga cieczy i może pracować w szerokim zakresie temperatur, ale jego charakterystyka zależy od amplitudy, kierunku, historii i napięcia wstępnego. Wymaga sprawdzenia osiadania i trwałości zmęczeniowej; nie można przypisać mu stałego $c$ bez określenia zakresu ruchu.

Sprężyna metalowa z osobnym tłumikiem rozdziela nośność od dyssypacji i pozwala osobno kształtować $k$ oraz $c$, lecz dodaje interfejsy i tryby uszkodzeń. W tłumiku lepkim zależność siły od prędkości wynika również z lepkości, temperatury i geometrii przepływu; projekt dla środowiska próżniowego wymaga ponadto oceny szczelności i zachowania cieczy.21 Trzeba także sprawdzić stan po utracie tłumienia.

Tarcie wprowadza próg ruszenia, stick–slip, zużycie i zależność od docisku; może ograniczać duże odpowiedzi, ale zniekształca małoamplitudową FRF.19 Warstwa lepkosprężysta zależy od temperatury i częstotliwości oraz wymaga właściwego pola odkształcenia, trwałego połączenia i odporności na pełzanie. Program badań materiału tłumika Hubble’a łączył temperaturę, częstotliwość, przemieszczenie, liczbę cykli i pomiar zespolonej sztywności przed oraz po starzeniu.20 W obu mechanizmach amplituda i historia cykli należą do specyfikacji.

Wybór jest zatem hipotezą, że element dostarczy wymaganą sztywność, stratność i skok bez nieakceptowalnego dryfu, zanieczyszczenia lub nowej drogi reakcji. Hipoteza wyznacza plan pomiarów i prób: charakterystykę elementu pod reprezentatywnym dociskiem, amplitudą i temperaturą; reakcje i ruch kompletnego zespołu; kontrolowane zmiany dróg równoległych; badania po starzeniu i cyklach; wreszcie regresję funkcjonalną zatwierdzonej konfiguracji.

13. Mapa dowodowa pomiaru

Pomiar staje się dowodem dopiero po związaniu go z pytaniem i decyzją. Łańcuch powinien mieć postać: wymaganie lub hipoteza → miara → przewidywana sygnatura → kanały → kryterium → decyzja. Sam wykres, nawet poprawnie obliczony, nie rozstrzyga, czy chroniona funkcja zachowała wymagany margines. Terminologię drgań, udaru i monitorowania stanu porządkuje ISO 2041,5 budżet niepewności opisują rozdziały 3.5–3.6 podręcznika Podstawy pomiarów,12 a próbkowanie, tor akwizycji, dyskretną transformatę Fouriera i skutki doboru okna omawiają rozdziały 1.5, 5.1 i 6.3–6.4 książki Systemy pomiarowe i podstawy analizy sygnałów.11

Wymaganie lub pytanie Miara Przewidywana sygnatura Niezbędne kanały Kryterium Decyzja
Czy ekspozycja losowa przekracza dopuszczalną intensywność w paśmie? PSD odpowiedzi, wkład pasma do wariancji, RMS, czas ekspozycji i miara uszkodzenia właściwa dla funkcji różnica względem widma odniesienia lub zbliżenie do limitu w określonym paśmie, możliwe wąskie piki rezonansowe zsynchronizowane wejście i odpowiedzi w punktach związanych ze stanem granicznym limit PSD/RMS oraz osobne kryterium trwałości lub funkcji spełnione dla określonego profilu i czasu utrzymać profil, uzasadnić ograniczenie wejścia albo zmienić konstrukcję
Czy zostanie wykorzystany prześwit lub wystąpi kontakt? przebieg czasowy przemieszczenia względnego, prędkości i przyspieszenia dojście do wartości prześwitu, zmiana nachylenia odpowiedzi, impuls wysokoczęstotliwościowy po kontakcie ruch podstawy i chronionej masy, pomiar względny oraz niezależny kanał siły lub styku górna granica maksymalnego przemieszczenia względnego mniejsza od dolnej granicy kierunkowego prześwitu o wymagany zapas; brak kontaktu zaakceptować albo zmienić sztywność, tłumienie, ogranicznik lub prowadzenie przewodów
Czy udar pobudza niebezpieczny zakres dynamiczny? określony rodzaj SRS oraz przebieg czasowy maksimum odpowiedzi oscylatorów odniesienia SRS w paśmie istotnym dla chronionej funkcji; pik nie identyfikuje sam przez się modu wyrobu wejście udarowe i odpowiedzi konstrukcji, poprawny zakres bez nasycenia obwiednia właściwego SRS i odpowiedź czasowa mieszczą się w kryteriach przyjąć dowód wyłącznie dla wskazanego tłumienia, rodzaju SRS i konfiguracji
Czy zmieniła się droga transmisji lub parametr modalny? FRF: podatność, mobilność albo akcelerancja; koherencja jako kontrola jakości estymacji przesunięcie rezonansu, zmiana fazy, amplitudy lub antyrezonansu zmierzona siła wejściowa oraz odpowiednio przemieszczenie, prędkość albo przyspieszenie zmiana większa od powtarzalności i budżetu niepewności, przy wiarygodnej estymacji FRF skorelować model, zbadać mocowanie lub odrzucić zmianę konfiguracji
Czy obserwowany ruch odpowiada przewidywanym modom? eksperymentalna identyfikacja modalna wejście–wyjście częstotliwości, tłumienie i postacie własne zgodne lub rozbieżne z modelem zmierzona siła i wielopunktowa odpowiedź amplitudowo-fazowa; metody output-only stanowią odrębną klasę zgodność parametrów istotnych dla przewidywanego użycia modelu zatwierdzić zakres modelu albo skorygować konkretną hipotezę fizyczną
Jak konstrukcja porusza się podczas pracy i która droga jest kandydatem? kształt odkształcenia roboczego (operational deflection shape, ODS) kołysanie, lokalne ugięcie, ruch związany z wiązką lub kontakt zależny od warunków pracy synchroniczne odpowiedzi w odpowiedniej siatce punktów oraz sygnał odniesienia; przy wniosku o reakcjach także siła, odkształcenie lub niezależny bilans obserwowany wzorzec wytypowuje drogę i jest powtarzalny; przyczynę potwierdza bilans albo kontrolowana zmiana tej drogi zaplanować interwencję rozstrzygającą, pomiar reakcji lub identyfikację modalną

Dla sygnału $x(t)$ o średniej $\bar{x}$ PSD dotyczy zwykle fluktuacji $x-\bar{x}$. Przy konsekwentnie przyjętej normalizacji jednostronnej całka PSD w paśmie $[f_1,f_2]$ daje wkład tego pasma do wariancji; pierwiastek całki po całym rozpatrywanym paśmie jest wartością RMS składowej zmiennej.22

$$\sigma^2_{[f_1,f_2]}=\int_{f_1}^{f_2}G_{xx}(f)\,df.$$

Po scałkowaniu całego zdefiniowanego pasma sygnału otrzymuje się całkowitą wariancję $\sigma^2$, a całkowita wartość skuteczna spełnia zależność

$$x_\mathrm{RMS}=\sqrt{\bar{x}^2+\sigma^2}.$$

RMS jest więc pierwiastkiem wariancji tylko po usunięciu średniej albo dla procesu o średniej równej zeru. PSD nie opisuje kolejności zdarzeń, a równe RMS nie oznacza równego maksimum, czasu przebywania przy dużej amplitudzie ani równego zużycia. Fizyczne własności procesu należy odróżnić od ich estymaty: estymata PSD i jej niepewność zależą od przyjętego pasma, rozdzielczości częstotliwościowej, okna, liczby efektywnie niezależnych uśrednień i czasu rejestracji. Nie są to wyłącznie ustawienia wykresu — wpływają na obciążenie estymatora, jego wariancję oraz zdolność rozdzielania bliskich składowych.23 Interpretacja statystyczna wymaga procesu dostatecznie stacjonarnego w analizowanym przedziale; rozruch, dryf temperatury lub przejście do kontaktu należy wydzielić i analizować czasowo.

SRS również nie jest jedną uniwersalną krzywą. Trzeba podać, czy przedstawia odpowiedź przyspieszeniową, przemieszczenie względne, pseudoprędkość czy inną wielkość, a także tłumienie oscylatorów, warunki początkowe, sposób traktowania dodatnich i ujemnych maksimów oraz algorytm obliczenia.8 SRS traci fazę i pełną historię czasową, dlatego nie wystarcza do oceny luzu, liczby kontaktów ani trwałego przesunięcia.

FRF musi być nazwana przez parę wejście–odpowiedź: podatność to przemieszczenie na jednostkę siły, mobilność — prędkość na jednostkę siły, a akcelerancja — przyspieszenie na jednostkę siły. Ich wartości mają sens przy znanej konwencji amplitudy, fazy i jednostek; eksperymentalne wyznaczanie FRF przez wymuszenie udarowe wymaga kontroli siły oraz odpowiedzi.9 FRF zakłada w przybliżeniu liniowy i stacjonarny układ w czasie pomiaru. Zależność wyniku od amplitudy, kierunku przemiatania lub historii wskazuje, że jedna liniowa FRF nie opisuje całego zakresu pracy.

Koherencja zwyczajna dla pary sygnałów jest zależną od częstotliwości wielkością $\gamma^2_{xy}(f)=|G_{xy}(f)|^2/[G_{xx}(f)G_{yy}(f)]$, wyznaczaną z estymat autospektrów i widma wzajemnego; jej wartość i rozrzut zależą zatem także od segmentacji, okien i uśredniania.24 W eksperymencie wejście–wyjście jest wskaźnikiem zgodności danych z liniową relacją między konkretnym zmierzonym wejściem i odpowiedzią, nie uniwersalnym miernikiem „jakości pomiaru”. Niska wartość może wynikać z szumu, niezmierzonego lub wielokrotnego wejścia, niewystarczającego wzbudzenia, przecieku widmowego albo nieliniowości; dlatego jej interpretacja musi korzystać z konfiguracji próby i innych kanałów.25 Wysoka koherencja nie dowodzi reprezentatywności toru pomiarowego, poprawnego montażu czujnika, kompletności zestawu wejść ani braku systematycznego błędu. Identyfikacja modalna wyodrębnia parametry modelu modalnego, natomiast ODS przedstawia ruch roboczy będący mieszaniną modów, odpowiedzi wymuszonej i wpływu wielu wejść. ODS może wskazać miejsce problemu, lecz sama nie jest postacią własną ani dowodem drogi przenoszenia siły.

14. Próby jako eksperyment

Program próby powinien rozdzielać hipotezy, a nie tylko wykonywać zadany profil. Sinus ujawnia rezonanse, fazę i zależność od amplitudy; różnica między przemiataniem w górę i w dół może wskazywać luz, tarcie lub zmianę temperatury. Próba losowa kontroluje zadane PSD w punktach sterowania, ale energia rzeczywiście wprowadzona do wyrobu zależy od impedancji wyrobu, przyrządu mocującego i stanowiska. Udar sprawdza odpowiedź przejściową, maksymalne przemieszczenie względne i kontakt. Wymuszenie akustyczne działa na rozłożone powierzchnie, dlatego nie jest zamienne z ruchem podstawy nawet wtedy, gdy pojedynczy kanał przyspieszenia osiąga podobny poziom.

Reprezentatywność dotyczy przede wszystkim dróg wprowadzania obciążenia: punktów i geometrii mocowania, podparcia podstawy, ciśnienia akustycznego, przewodów, masy, środka masy, temperatury oraz stanu blokad. Zakres informacji wymienianych przy stosowaniu podatnych systemów mocowania porządkuje ISO 2017-1; norma ta nie zastępuje kompletnej metodyki projektowania ani kwalifikacji.6 Charakterystyka samej podpory może wymagać pomiaru sztywności dynamicznej przy reprezentatywnym obciążeniu wstępnym i sposobie zamocowania.10

Przyrząd mocujący, wyrób i ruchoma część stanowiska tworzą wspólny układ. Mody przyrządu mocującego, lokalna elastyczność interfejsu albo niejednorodny ruch punktów mocowania mogą wytworzyć pozorną odpowiedź wyrobu. Sterowanie wielopunktowe musi zatem określać sposób łączenia kanałów, tolerancje, pasma oraz reakcję na rozbieżność między punktami. Raport powinien zachować rzeczywiście zrealizowany profil, nie tylko profil zadany.

Force limiting i response limiting są odrębnymi kryteriami ochronnymi: pierwsze steruje na podstawie zmierzonej siły interfejsowej, drugie na podstawie odpowiedzi wybranego punktu lub funkcji. Osiągnięcie limitu może wywołać zmniejszenie wejścia w określonym paśmie, czyli notching, ale także wstrzymanie albo przerwanie próby. Notching może być ponadto zaplanowany przed próbą na podstawie uzasadnionego modelu. Pojęcia są więc powiązane, lecz nie tożsame; NASA-RP-1403 definiuje próbę force-limited przez pomiar i kontrolę siły, wyraźnie odróżniając ją od samego analitycznego wyznaczenia wcięcia profilu.14 Każdy limit trzeba powiązać z modelem środowiska użytkowego, czujnikiem, niepewnością i chronionym stanem granicznym.

Korelację planuje się przed próbą: przewiduje się częstotliwości, fazę, skok, siły i reakcje, a następnie określa dopuszczalne odchylenia. Kryteria przerwania powinny obejmować przekroczenie odpowiedzi funkcjonalnej, ruchu względnego, siły, temperatury, limitu stanowiska, utratę kanału krytycznego i sygnaturę rozwijającego się uszkodzenia. Po przerwaniu zabezpiecza się dane i konfigurację przed demontażem.

Stanowisko wymaga osłon i wtórnego zatrzymania odpowiedniego do możliwych następstw awarii. Należy uwzględnić energię sprężystą zgromadzoną w podporach, napinaczach i ogranicznikach, wyrwanie elementu lub czujnika oraz niekontrolowany ruch wyrobu, przyrządu mocującego albo części wzbudnika. Bezpieczny stan musi być określony także dla zaniku zasilania, utraty sterowania i przerwania mocowania.

15. Hierarchia modeli i wiarygodność

Obliczenie ręczne sprawdza rząd wielkości częstotliwości, ugięcia, siły i skoku. Model skupiony dodaje kilka mas, równoległe drogi i ograniczniki. Model 6-DOF służy do oceny translacji, obrotu i reakcji podpór. MES jest uzasadniony, gdy o decyzji rozstrzygają elastyczne mody konstrukcji, lokalne odkształcenia albo podatność podstawy. Nieliniowy model kontaktu jest potrzebny po pojawieniu się luzu, tarcia, histerezy, odrywania lub ogranicznika. Wyższy poziom szczegółowości nie naprawia błędnej masy, warunków brzegowych ani charakterystyki połączenia.

Walidacja dotyczy przewidywanego użycia modelu. Zgodność częstotliwości niskich modów nie waliduje automatycznie sił kontaktu ani odpowiedzi udarowej. Dla każdej decyzji określa się wielkości wyjściowe, zakres amplitudy, temperaturę, konfigurację i akceptowalną niepewność. Wiarygodność obejmuje weryfikację obliczeń, walidację względem właściwych danych, ocenę niepewności oraz granice zastosowania.3,4

Niepewność trzeba propagować do kryterium, nie tylko do parametrów wejściowych. Rozrzuty masy, sztywności, tłumienia, docisku i temperatury mogą tworzyć inne najgorsze kombinacje dla przemieszczenia względnego, siły i obrotu. Budżet powinien rozdzielać składniki ocenione statystycznie i z innych danych, a wynik wskazywać przyjęty współczynnik rozszerzenia.12 Zarządzanie konfiguracją wiąże wynik z wersją geometrii, właściwościami materiałów, partią podpór, modelem przyrządu mocującego, oprogramowaniem, ustawieniami solvera i zestawem danych korelacyjnych. Aktualizacja parametru po próbie wymaga udokumentowanej hipotezy fizycznej i sprawdzenia w kanałach niewykorzystanych do strojenia.

16. Kwalifikacja, akceptacja, diagnostyka i cykl życia

Na potrzeby tego artykułu kwalifikacja oznacza wykazanie zdolności dojrzałego projektu do zachowania funkcji z wymaganym marginesem, a akceptacja — ocenę konkretnego egzemplarza ukierunkowaną na wykrycie wad wykonania bez nieuzasadnionego zużywania jego trwałości. Próba rozwojowa służy uczeniu się i korelacji, natomiast diagnostyka rozstrzyga hipotezę o zmianie stanu. Są to definicje robocze; dokument programu może stosować terminologię szczegółową dla własnej branży. Różne cele oznaczają różne poziomy, czasy, konfiguracje i reguły oceny, dlatego jedna ekspozycja nie uzyskuje automatycznie wszystkich tych znaczeń.

NASA-STD-7001C odnosi się do określonej klasy ładunków kosmicznych i nie powinien być uogólniany na dowolne urządzenie, pojazd ani system rakietowy.1 NASA-STD-7002B należy opisywać jako wydanie z 2018 roku z Change 1/reaffirmation z 2023 roku; jego logikę wymagań prób można stosować tylko po sprawdzeniu zakresu i adekwatności do danego programu.2

Diagnostyczna miara ma wynikać z hipotezy. Analiza rzędów jest właściwa, gdy podejrzewana składowa pozostaje związana z prędkością obrotową i dostępny jest wiarygodny sygnał fazy lub obrotów. Analiza obwiedni odpowiada na hipotezę o okresowych impulsach modulujących pasmo rezonansowe, na przykład wskutek cyklicznego kontaktu. Rozdziały 3.1, 4.1, 4.4 oraz 4.10–4.12 książki Utrzymanie ruchu w przemyśle pokazują dobór wielkości, widma, fazy i obwiedni do różnych symptomów oraz znaczenie porównywalnego stanu pracy.13 Żadna z tych miar nie jest uniwersalnym wskaźnikiem stanu: wynik zależy od pasma, prędkości, obciążenia, toru pomiarowego i konfiguracji.

Dowód cyklu życia musi identyfikować partię elementu, jego pozycję i orientację, historię obciążeń, temperaturę, czas magazynowania, demontaże oraz naprawy. Zamiana podpór miejscami może zniszczyć ciągłość trendu, ponieważ każda pozycja ma inną reakcję i temperaturę. Po naprawie nie wystarcza nazwa wymienionej części: trzeba odtworzyć napięcie wstępne, podkładki, momenty połączeń, prowadzenie przewodów i stan powierzchni. Kryterium ponownego dopuszczenia powinno łączyć kontrolę funkcji, inspekcję, powtarzalny pomiar dynamiczny i ocenę wykorzystanej trwałości.

17. Reguła decyzji i granica jawności

Końcowa decyzja nie może opierać się na najkorzystniejszym pojedynczym kanale. Łączy odpowiedź chronionej funkcji, największy skok, siłę transmisji, obrót, reakcje wszystkich podpór, trwałość oraz niepewność pomiaru i modelu. Wynik jest pozytywny dopiero wtedy, gdy wszystkie istotne stany graniczne pozostają spełnione w tej samej dopuszczonej konfiguracji i obwiedni środowiska. Spadek przyspieszenia przy jednoczesnym wykorzystaniu prześwitu, przeciążeniu jednej podpory albo skróceniu trwałości nie stanowi poprawy systemowej.

Odchylenie od przewidywania wymaga rozstrzygnięcia, czy pochodzi z wyrobu, przyrządu mocującego, toru pomiarowego, niepełnego modelu czy zmiany konfiguracji. Margines musi być porównany z niepewnością; nominalne spełnienie limitu bez zdolności rozróżnienia wyniku od przekroczenia nie zamyka dowodu. Decyzja powinna wskazywać zakres ważności, nierozstrzygnięte ryzyka, właściciela akceptacji i warunki ponownej oceny po zmianie.

Publikowanie wyników wymaga odrębnego zatwierdzenia zgodnego z zasadami ochrony informacji obowiązującymi właściciela systemu. Wrażliwa charakterystyka może zostać ujawniona przez zestawienie pozornie jawnych danych: częstotliwości własnych, tłumienia, poziomów prób, rozmieszczenia czujników, reakcji podpór, strojenia absorbera i progów przerwania. Ocenie podlega więc cały pakiet wraz z metadanymi, wykresami i pochodnymi modelu, a nie każdy fragment osobno. Materiał publiczny może wyjaśniać metodę i fizykę, lecz dane pozwalające odtworzyć zachowanie konkretnego współczesnego systemu powinny pozostać w kontrolowanym obiegu.

18. Synteza

Izolator, tłumik ani absorber nie mają samodzielnej „skuteczności”. Tworzą ją dopiero wraz z masą, geometrią, podstawą, przewodami, ogranicznikami, temperaturą i historią obciążeń. Z tego powodu dowód należy do określonej konfiguracji: katalogowa część dostarcza charakterystyki wejściowej, lecz nie przesądza o skoku, obrocie, reakcjach, trwałości ani zachowaniu funkcji po integracji.

Hierarchię modelu wyznacza pytanie. SDOF wystarcza do kontroli dominującego ruchu i podstawowego kompromisu transmitancja–skok; model 6-DOF jest potrzebny przy obrocie i nierównych reakcjach; podatna podstawa oraz lokalne mody prowadzą do MES; kontakt, tarcie i histereza wymagają prawa nieliniowego oraz przebiegu czasowego. Wyższy model jest uzasadniony dopiero wtedy, gdy prostszy nie odtwarza wielkości sterującej decyzją.

Hierarchię miar wyznacza stan graniczny. Ruch bezwzględny opisuje ochronę odbiornika, ruch względny — wykorzystanie skoku, siła interfejsowa — obciążenie drogi, a moc i energia — dyssypację. FRF porządkuje odpowiedź liniową, SRS rodzinę odpowiedzi na udar, ODS ruch roboczy, lecz żadna z tych reprezentacji nie zastępuje brakującej fazy, historii czasowej ani pomiaru chronionej funkcji.

Przejście od obliczenia do dowodu wymaga kolejno: jawnej hipotezy i przewidywanej sygnatury, reprezentatywnej konfiguracji, skalibrowanych kanałów, budżetu niepewności, testu rozstrzygającego oraz regresji po zmianie. Kryterium musi jednocześnie obejmować funkcję, skok, reakcje i trwałość. Dopiero ich łączne spełnienie w określonej obwiedni pozwala zamknąć decyzję; korzystny pojedynczy wykres pozostaje obserwacją, nie kwalifikacją.

Źródła