Redukcja modeli nieliniowych dużych konstrukcji przez dwuetapowe poddzielanie

Dlaczego zmniejszanie wielkich cyfrowych budowli ma znaczenie

Inżynierowie często symulują, jak duże konstrukcje, takie jak zakłady przemysłowe, mosty czy kadłuby samolotów, kołyszą się i drgają pod wpływem wiatru, trzęsień ziemi czy maszyn. Takie cyfrowe testy mogą obejmować setki tysięcy punktów ruchomych i trwać godziny lub dni na potężnych komputerach. W artykule przedstawiono sposób zmniejszenia takich rozległych modeli do znacznie mniejszych, które nadal zachowują się jak oryginał, nawet gdy konstrukcja ma kłopotliwe, silnie nieliniowe złącza i realistyczne, złożone formy tłumienia.

Rozbijanie olbrzymiej konstrukcji na mniejsze części

Punktem wyjścia jest obserwacja, że duże konstrukcje zwykle składają się z powtarzających się elementów: podobnych ram, kondygnacji czy paneli. Zamiast traktować cały budynek naraz, metoda dzieli go na podstruktury. Każdą podstrukturę analizuje się osobno, a następnie łączy przez siły na wspólnych granicach. Ta filozofia, znana jako poddzielanie (substructuring), od dawna stosowana jest w prostszych, liniowych układach, w których odpowiedź jest proporcjonalna do przyłożonych obciążeń. Wkład tej pracy polega na sposobie radzenia sobie z bardziej realistycznym zachowaniem, w którym niektóre złącza działają nieliniowo, a energia tracona na skutek tłumienia nie podąża za uproszczonymi, podręcznikowymi wzorcami.

Uchwycenie złożonego ruchu za pomocą prostych wzorców

Aby zmniejszyć rozmiar każdej podstruktury bez utraty istotnej fizyki, autor korzysta z pojęcia nieliniowych trybów własnych (nonlinear normal modes). W istocie tryb to charakterystyczny sposób, w jaki struktura preferuje drgać. Dla układów liniowych te tryby są proste i dobrze ustrukturyzowane. Gdy ruch staje się duży albo złącza zachowują się jak sztywne sprężyny reagujące w sposób sześcienny zamiast prostolinijny, te wzorce ulegają zagięciu i deformacji. Artykuł przedstawia matematyczny przepis, który reprezentuje każdy nieliniowy tryb jako gładką, zakrzywioną powierzchnię w przestrzeni wszystkich możliwych ruchów. Ruch każdego punktu podstruktury jest wyrażany jako wielomian tylko w kilku kluczowych przemieszczeniach i prędkościach zlokalizowanych na interfejsach, gdzie spotykają się podstruktury. To zamienia ogromny zbiór zmiennych w bardzo zwartą reprezentację, która wciąż odzwierciedla nieliniowy charakter złączy.

Utrzymanie równowagi statycznej i realistycznego tłumienia

Metoda dzieli odpowiedź każdej podstruktury na część dynamiczną, w której żyją nieliniowe tryby, oraz część statyczną, która radzi sobie z powolnymi odkształceniami spowodowanymi siłami na interfejsach. Dla części statycznej podejście zapożycza pomysły z istniejącej ramy zwanej dwuetapową metodą Craig–Bampton (dual Craig–Bampton). Tam kompatybilność między podstrukturami wymuszana jest przez siły na interfejsie, zamiast bezpośredniego sklejania przemieszczeń brzegowych. To prowadzi do mniejszych macierzy i większej elastyczności w łączeniu elementów. Ważnym udoskonaleniem tej pracy jest uwzględnienie ogólnych form tłumienia bezpośrednio w równaniach, zamiast zakładać, że tłumienie jest po prostu proporcjonalne do masy albo sztywności. W rezultacie zredukowany model może wiernie odtwarzać konstrukcje wyposażone w dodatkowe tłumiki lub materiały, które rozpraszają energię w sposób niejednorodny.

Testowanie pomysłu na cyfrowym budynku przemysłowym

Aby pokazać, że metoda jest praktyczna, autor stosuje ją do szczegółowego modelu stalowego budynku przemysłowego. Ramy budynku zawierają złącza modelowane jako sprężyny skrętne, których opór rośnie z sześcianem kąta obrotu — silna forma nieliniowości. Budynek jest wymuszany bocznie siłą sinusoidalną dostrojoną w pobliże jednej z jego naturalnych częstotliwości drgań. Najpierw pełny model MES rozwiązywany jest standardowym algorytmem całkowania w czasie, zużywając kilka setek sekund czasu obliczeń i setki megabajtów pamięci. Następnie budynek dzieli się na powtarzające się podstruktury ram i część pozostałą. Dla ram zachowywane są tylko cztery nieliniowe tryby, skoncentrowane na ruchu poziomym i skręcaniu najistotniejszych węzłów. Rozwiązanie tego zredukowanego układu daje przebiegi przemieszczeń, które niemal idealnie pokrywają się z wynikami pełnego modelu, skracając czas obliczeń o około dwie trzecie i znacznie redukując zużycie pamięci.

Dlaczego mniejsza liczba trybów nadal daje wiarygodne odpowiedzi

Badanie analizuje też, jak dokładność zależy od liczby i doboru nieliniowych trybów. Przy użyciu tylko jednego trybu błąd w przewidywanym ruchu jest większy. Dodanie drugiego trybu, który bezpośrednio obejmuje złącze o zachowaniu sześciennym, prowadzi do gwałtownego spadku błędu, co podkreśla znaczenie uwzględnienia stopni swobody, gdzie nieliniowość jest najsilniejsza. Przy trzech i czterech trybach błąd nadal maleje do poziomów uważanych za bardzo małe w inżynierskim projektowaniu, przy zachowaniu zwartego modelu. Drugi zestaw symulacji dodaje zewnętrzne tłumiki, które tworzą silnie nieproporcjonalny wzorzec tłumienia. Nawet w tym bardziej wymagającym scenariuszu zredukowany model ściśle śledzi pełne rozwiązanie, oferując przy tym znaczące oszczędności czasu i pamięci.

Co to oznacza dla przyszłych cyfrowych konstrukcji

Mówiąc prosto, artykuł pokazuje, jak przekształcić nieporęczny cyfrowy model budynku w zwinny reprezentant, który reaguje niemal identycznie na wstrząsy, nawet gdy jego złącza zachowują się w skomplikowany, nieliniowy sposób, a straty energii są nieregularne. Łącząc poddzielanie, nieliniowe wzorce drgań oraz formułę uwzględniającą tłumienie, metoda otwiera drogę do szybkich, a jednocześnie wiarygodnych symulacji bardzo dużych konstrukcji. To może pomóc inżynierom przeprowadzać znacznie więcej scenariuszy „co-jeśli”, optymalizować projekty i badać nowe materiały i urządzenia, bez ograniczeń wynikających z nadmiernych kosztów obliczeniowych.

Cytowanie: Flores, P.A. Nonlinear model reduction for large-scale structures via dual substructuring. Sci Rep 16, 9286 (2026). https://doi.org/10.1038/s41598-026-38015-7

Słowa kluczowe: dynamiczna analiza struktur, redukcja modeli, drgania nieliniowe, analiza metodą elementów skończonych, poddzielanie