Czułość siatki i weryfikacja eksperymentalna dla losowo zaprojektowanych metamateriałów akustycznych opartych na komorach o dowolnej geometrii zaprojektowanych za pomocą 2D FEM

Uciszanie hałasu za pomocą malutkich labiryntów

Współczesne życie jest głośne: od hal produkcyjnych po biura typu open space, niechciany hałas szkodzi zdrowiu i koncentracji. Inżynierowie zwracają się ku „metamateria łom akustycznym” – precyzyjnie zaprojektowanym strukturom, które tłumią dźwięk w sposób niedostępny dla zwykłej pianki czy włókien szklanych. W niniejszym artykule opisano nową, szybszą metodę projektowania szczególnej klasy tych materiałów, wykorzystujących kanały przypominające labirynty, które pochłaniają dźwięk przy jednoczesnym zachowaniu kompaktowości i niskiej masy.

Budowanie inteligentnych pułapek na dźwięk

Metamateriały akustyczne to powtarzalne elementy wypełnione wąskimi komorami i kanałami, które manipulują falami dźwiękowymi. Wiele współczesnych, wysokowydajnych pochłaniaczy opiera się na rezonatorach – małych kieszeniach i rurkach drgających na określonych częstotliwościach, które przekształcają energię akustyczną w ciepło. Omawiane tutaj projekty to metamateriały „oparte na komorach”, w których dźwięk jest wymuszony do przejścia przez kręte, labiryntowe drogi powietrzne. Gdy fala dźwiękowa przeciska się przez te wąskie przejścia, tarcie i niewielkie zmiany temperatury przy ściankach odbierają jej energię, zmniejszając przenikanie hałasu.

Dlaczego konwencjonalne symulacje natrafiają na barierę

Aby zaprojektować takie skomplikowane struktury, badacze zwykle korzystają z potężnych symulacji komputerowych opartych na metodzie elementów skończonych (FEM). Modele te śledzą, jak porusza się dźwięk i jak tracona jest energia w cienkich „warstwach przyściennych” powietrza przylegających do ścian kanałów. Jednak gdy geometria jest złożona i naprawdę trójwymiarowa, wierne odwzorowanie efektów termowiskowych wymaga ogromnej liczby punktów obliczeniowych, czyli elementów siatki. W praktyce pełny model 3D, który rozdziela te warstwy, może zajmować dni obliczeń dla pojedynczego projektu, co uniemożliwia systematyczną optymalizację wielu kształtów.

Spłaszczanie projektów 3D do map 2D

Autorzy proponują inną strategię: reprezentować komórkę metamateriału 3D jedną przekrojową mapą 2D i symulować tylko ten wycinek. Koncentrują się na strukturach, które można otrzymać przez ekstruzję płaskiego wzoru prostopadle do płaszczyzny, na przykład kanałów labiryntowych. Każdy projekt jest kodowany jako prosty bitmapowy obraz czarno-biały, gdzie jeden piksel odpowiada kwadratowi 2 mm będącemu albo materiałem stałym, albo powietrzem. To zamienia problem projektowy w układanie pikseli na siatce zgodnie z podstawowymi regułami (ciągłe ścieżki powietrza, brak odizolowanych kieszeni, brak jednopikselowych „kolców” materiału), a następnie użycie 2D modelu FEM uwzględniającego straty termowiskowe do przewidzenia, ile dźwięku struktura pochłonie w określonym zakresie częstotliwości.

Testowanie dokładności i skracanie obliczeń

Aby sprawdzić, czy model płaski może zastąpić pełny model 3D, badacze najpierw porównali kilka podejść na prostej strukturze z zaledwie dwoma rezonatorami. Rozważali wzory analityczne (metodę macierzy przejścia), standardowe 3D FEM, ich zredukowany model 2D oraz rzeczywiste pomiary w rurce impedancyjnej. Symulacja 3D z pełną fizyką termowiskową zajęła niemal sześć dni obliczeń i wciąż wykazywała zauważalne przesunięcia częstotliwości. Natomiast 2D model termowiskowy uruchamiał się w ciągu kilku minut i trafiał w mierzony częstotliwościowy szczyt pochłaniania w granicach około ćwierci procenta. Zachęceni tym, autorzy przeszli do bardziej złożonych, losowo generowanych geometrii labiryntowych kodowanych jako mapy 32×32 piksele.

Jak gruba może być siatka, żeby nadal działała?

Ponieważ większość kosztów obliczeniowych wynika z rozdzielania siatki przy ściankach, zespół systematycznie zmieniał dwa czynniki skalujące kontrolujące, jak cienka jest pierwsza warstwa przyścienna i ile takich warstw stosuje się. Dla dwudziestu różnych struktur typu labirynt i po siedemdziesiąt pięć ustawień siatki każdy zmierzyli, jak bardzo przewidywane krzywe pochłaniania dźwięku zmieniają się względem bardzo drobnej siatki „referencyjnej”. Stwierdzili, że nawet gdy siatka w warstwie przyściennej była znacznie grubsza, średni błąd w przewidywanym pochłanianiu pozostawał poniżej 0,5% dla szerokiego zestawu ustawień, podczas gdy liczba niewiadomych w obliczeniu spadała o ponad 70%. W końcu wydrukowali 3D sześć nowych struktur i porównali model 2D z pomiarami w rurce. Model przewidywał częstotliwości rezonansów ze średnim błędem około 2,6%, przy czym większe różnice pojawiały się głównie w wysokości szczytów — prawdopodobnie z powodu chropowatości powierzchni i strat materiałowych w drukowanym plastiku.

Co to oznacza dla przyszłej kontroli hałasu

Dla czytelnika niebędącego specjalistą główne przesłanie jest takie, że autorzy pokazali sposób na zamianę bardzo ciężkiego problemu symulacji dźwięku 3D w znacznie lżejszy problem 2D, bez rezygnacji z praktycznej dokładności dla szerokiej klasy pochłaniaczy labiryntowych. Pracując na pikselowych planach i starannie dostrojonych siatkach, mogą eksplorować znacznie więcej kandydatów do projektów na zwykłych komputerach, torując drogę do automatycznej optymalizacji, a nawet generowania nowych metamateriałów akustycznych wspomaganego przez sztuczną inteligencję. Choć metoda nie obejmuje wszystkich możliwych geometrii i jak dotąd testowana była w ograniczonym paśmie częstotliwości, stanowi potężne skrócenie drogi ku cichszym maszynom, pomieszczeniom i urządzeniom zbudowanym z przemyślanie ułożonych, dźwiękochłonnych labiryntów.

Cytowanie: Książek, P., Chojnacki, B. Mesh sensitivity and experimental verification for randomized arbitrary geometry cavity-based acoustic metamaterials designed with 2D FEM simulations. Sci Rep 16, 6873 (2026). https://doi.org/10.1038/s41598-026-38139-w

Słowa kluczowe: metamateriały akustyczne, pochłanianie dźwięku, modelowanie metodą elementów skończonych, struktury labiryntowe, rurka impedancyjna