## Elastic wave control in reticulated plates using Schwarz primitive cells

Aida Hejazi Nooghabi, 1 Henrik R. Thomsen, 2 Bao Zhao, 1 and Andrea Colombi 3

1 Department of Civil, Environmental and Geomatic Engineering, ETH Zurich, Zurich, Switzerland 2 Institute for Geophysics, ETH Zurich, Zurich, Switzerland 3 Zurich University of Applied Sciences (ZHAW), Winterthur, Switzerland

## I. ABSTRACT

In this work the Schwarz primitive unit cell is used as the building block of different types of metastructures for steering and focusing elastic vibrations. The emergence of a Bragg-type bandgap when constructing a 2D plate from such unit cells is experimentally validated. It is demonstrated that both increasing the mass and the porosity of the Schwarz primitive lead to a decrease in the frequency of the out-of-plane propagating wave targeted in this study. By arranging these modified Schwarz primitive unit cells in constant and graded layouts, 2D plates with an embedded metabarrier and a metalens are numerically designed. The metabarrier protects an interior area of the plate from the propagating waves on a wide frequency band ( ∼ 1.4 to ∼ 3.4 kHz). Equally, the refractive index profile necessary for gradient index lenses is obtained via a progressive variation of the added mass or, alternatively, the porosity of the unit cell over a rectangular area. For the first time, bending of the out-of-plane mode towards the focusing point is practically validated in a challenging mesoscale experiment requiring the assembly of different 3D printed sections of the plate. The increased porosity design is advantageous not only in terms of overall lightweight, but also towards additive manufacturing as it requires less material.

keywords: Schwarz primitive cell, elastic wave propagation, bandgap, metabarrier, Luneburg lens

## II. INTRODUCTION

Metamaterials - with their engineered structures - have been heavily exploited in the past couple of decades both for static and dynamic applications. They have proven useful for a variety of applications such as sound absorption [1], elastic wave guiding [2], energy harvesting [3, 4], mechanical neural smart devices [5], and seismic and vibration shielding [6-8], to name a few. One of the most important features of such structures is the so-called bandgap [9], defined as the frequency range where the propagation of waves is prohibited [9]. The well-known mechanisms behind this feature are either Bragg scattering [10], local resonance [11, 12] or the coalescence of the two [13-16]. While Bragg scattering emerges due to the periodicity in the arrangement of the unit cells [9] constituting the metamaterial structure, the local resonance mechanism relies on the resonance of the unit cell and is independent upon its arrangement. Local resonances could appear in both periodic and non-periodic distributions of such cells [17]. The opening frequency of the bandgap due to Bragg scattering in phononic crystals is determined by the relationship between the wavelength of the propagating wave and the spacing between unit cells. Contrary, in locally resonant metamaterials, the opening frequency is only dependent on the resonant frequency of a single resonator. Additionally, the frequency range of the bandgap in phononic crystals is often wider than the one in resonant metamaterials. These complementary characteristics have led to studies where a coalescence of the Bragg scattering and local resonance has been implemented to obtain tunable wide bandgaps [16].

In acoustics, metamaterials have been widely used for noise absorption and wave guiding, whereas, in elasticity a broader range of applications such as guiding, focusing, and energy harvesting have been explored. More specifically, elastic wave focusing based on GRadient INdex (GRIN) lenses has been used in a variety of structures and gradient profiles. For instance, Maxwell and Luneburg lenses have been numerically tested based on lattice structures in 2D by leveraging both Bragg and local resonance mechanisms in order to focus the out-of-plane waves in the kHz regime [16]. A GRIN metasurface for microwave imaging has been numerically and experimentally tested in 1D and 2D at the GHz regime [18] and at kHz regime, focusing of surface acoustic waves has been experiementally achieved based on resonant metasurfaces [19]. A wider variety of GRIN lenses including not only Maxwell and Luneburg, but also Eaton and 90 ◦ rotating lenses have been built based on the variation in the thickness of the plates for the flexural waves [20].

In this work, we investigate how the Schwarz primitive cell [21] which belongs to the family of Triply Periodic Minimal Surfaces (TPMS) can be leveraged in the design of meta structures. The Schwarz primitive is symmetric in 3D and is known to show significant strength to compression and fracture toughness, therefore making it a candidate with high potentials for scaffolds in biomimetic designs and TPMS-based scaffolds [22-24]. TPMS structures can offer enhanced mechanical properties compared to truss-like structures. The continuous surface geometry of TPMS structures and gradual change in curvature allows for an efficient load distribution and stress transfer across the surface. This minimizes stress concentrations, improving stress distribution and eventually makes the frame less prone to fatigue damages [24, 25]. The Schwarz primitive cell has previously been investigated for its buckling characteristics [26], numerically and experimentally validated for topological elastic wave guiding [27], and has also been explored for vibration mitigation purposes in the realm of acoustic black holes [28]. An in-depth investigation on the relationship between the volume ratio and bandgap behavior [29] as well as the acoustic properties of Schwarz primitive cell [30] and mechanical vibration bandgaps [31] have been performed. Despite such different types of investigations on the dynamic properties of the Schwarz primitive cell in the context of acoustic and elastic metamaterials [29-31], their potential in the realm of focusing and mitigation of elastic waves has not been thoroughly studied.

Interestingly, biological arrangements with geometries similar to TPMS have been observed [32]. For instance, a slice of a sea urchin skeletal plate displays similarities to the Schwarz primitive geometry [32, 33]. Also, the structure of amphiphilic membranes that separate oil and water, also known as the 'plumber's nightmare', form arrangements that resemble the Schwarz primitive surface [34]. Despite different types of investigations on the dynamic properties of TPMS (more precisely the Schwarz primitive cell) in the context of elastic metamaterials [27, 28, 35], their potential in the context of focusing and mitigation has not been thoroughly studied.

In this paper, we first study the dynamic behaviour of a Schwarz primitive cell by extracting the dispersion relation of a unit cell made of versatile plastic. We validate the presence of a bandgap both numerically and experimentally. We next build plate-like metastructure by locally varying the geometry in the medium to manipulate elastic waves for vibration mitigation (metabarriers) and focusing (GRIN lenses) purposes. The particle vibration of interest is dominated by out-ofplane motion, which is analogous to the A 0 mode of Lamb waves propagating in thin plates [36]. Through large scale, parallel finite element time domain simulations we numerically validate the performance of such devices based on adding mass and increasing porosity in the unit cell. We find that the increased porosity not only adds to the metastructure's lightweight characteristics, but also makes it very suitable for 3D printing (less material) and hands-on experimentation. The concept of increasing porosity as a design variable for a GRIN lens is additionally verified through laboratory experiments conducted on a 3D printed prototype. Our findings on the functionality of such a metastructure pave the way for further applications of bio-mimetic structures in relation to

## III. DISPERSION RELATION OF THE SCHWARZ PRIMITIVE CELL

## A. Unit cell design

The geometry of the TPMS structures can be defined by trigonometric functions and in the case of the Schwarz primitive cell, the 3D surface is obtained by solving the following implicit equation:

$$c o s ( X ) + c o s ( Y ) + c o s ( Z ) = C$$

where X = 2 πx/L x , Y = 2 πy/L y , Z = 2 πz/L z , L x , L y and L z are the size of the unit cell in the x , y , and z directions and C is a constant [21, 24]. In this study, the size of the unit cell along the x , y and z directions is 3 cm and the value of C is zero. The resulting geometry is generated in MATLAB (See Figure. 1a). By adding an offset, the generated surface is thickened by 0.5 mm and the obtained geometry is then exported in STL format. The selected dimensions and thickness of the unit cell ensure the wave propagation properties required for designing metastructures, making it optimal for both experimental and numerical analyses. For the following time domain simulations, the full geometry is processed in Coreform Cubit, where the solid hexahedral mesh elements are generated. Once the geometry is built, thickened, healed, and meshed, the dispersion relation of the cell can be numerically extracted. In the following, we show the dynamic behaviour of the Schwarz primitive unit cell through dispersion relations obtained both numerically and experimentally.

## B. Numerical validation of the dispersion relation

Through simulations, we obtain the dispersion relation i.e., the relationship between the frequency f and real part of the wave number k based on two different methods, namely as eigenfrequency analysis [37] and 2D FFT method [38]. The eigenfrequency analysis is performed using the software COMSOL Multiphysics 6.1, where the dispersion relation of the unit cell discretized into solid tetrahedral mesh elements is solved through the finite element method. This analysis is based on the Bloch theorem, where Bloch periodic boundary conditions are applied to the unit cell and the dispersion relation ω = ω k ( ) is solved by varying the wave vector k along direction of the wave propagation [39]. We choose a conventional cubic unit cell and focus only on the wave propagation along ΓX . As the out-of-plane mode, propagating along the z -direction, is of interest, we retrieve and show the propagating modes only along this direction, however, in COMSOL all the modes (including in-plane and out-of plane) are generated (See Figure. 1b).

In order to correctly identify which modes are excited by an out-of-plane source, we next extract the dispersion relation along the z -direction through time domain simulations. The finite reticulated plate is obtained by tessellating 30 unit cells along both the y and z directions. Solid hexahedral mesh elements are generated through the Python APIs of the Coreform Cubit software (See Figure. 1a). In this set of computations, the material properties were set to: Young's modulus E = 1.0 GPa, Poisson's ratio ν = 0.39 and mass density ρ = 930 kg/m . The model shown here consists of 3 approx 10 6 hex elements. Time domain simulations are all performed in SPECFEM3D, an opensource code that implements a highly parallelised spectral element method for solving elastodynamic equations [40]. Since we are interested in the out-of-plane mode, only the polarization along the

FIG. 1: (a) The Schwarz Primitive unit cell with a length of 3 cm along the x , y and z directions and the configuration of the 2D plate used in time domain simulations built by tessellating the unit cell in two directions. (b) Dispersion relation of the unit cell along the ΓX direction obtained via Bloch analysis in COMSOL (black curves) and the outcome of 2D Fourier transform applied to time domain data (colored background plotted in dB). (c) Imaginary part of the wavenumber obtained from the time domain data. (d) Transmissibility retrieved from time domain data. (e) Mode shapes at the opening and closing frequencies of the bandgap.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00590v1-figure-86783228-f8de-4a9a-8787-177b8be91956.jpg)

x -axis is excited at points marked in Figure.1a. The dispersion relation for the given mode can be obtained by extracting the displacement field along the propagation direction and applying a 2D Fourier transform to the data (See Figure. 1b). The dark orange features show the propagating out-of-plane modes. Both the Bloch and time-domain analyses reveal a bandgap between 2.021 and 3.437 kHz for the A 0 mode. The results of both methods displayed in Fig. 1b suggest a Bragg-type mechanism for the observed bandgap of the out-of-plane modes. The latter can be validated via recovering the imaginary part of the wavenumber using the IWC method proposed by [41]. A behaviour typical of the Bragg bandgap is revealed (See Figure. 1c). The transmissibilty shown in Figure. 1d also indicates the presence of the bandgap already identified in Figure. 1b. Mode shapes at the opening and closing frequencies of the bandgap (indicated by a circle and triangle in Figure. 1b) are shown in Figure. 1e. While the opening mode of the bandgap is a purely out-of-plane mode, the closing mode includes a mixture of the out-of-plane and a torsional mode along the z direction.

## C. Experimental validation of the bandgap

The versatility of additive manufacturing, also for thin walled structures like this one, allows to validate the bandgap of the Schwarz primitive cell experimentally. In order to determine the size of the plate (i.e., the number of cells along the two directions) that allows the wavefield to be identified but also satisfies printability constraints in terms of size, we must perform a numerical simulation on a finite structure in the frequency domain using COMSOL Multiphysics 6.1. A plate consisting of 10 by 6 unit cells along the y and z directions, respectively mounted on top of a thin support plate of 4 mm thick necessary to attach the structure to an electrodynamic shaker, provided promising results by considering the elastic parameters defined earlier and an isotropic loss factor of η = 0 002. The structure is excited at discrete frequencies between 1 and 5 kHz with . a constant displacement ( u ) along the x direction through the supporting plate ( u =( u x ,0,0)). The transmissibility is then computed in dB:

$$T r a n s m i s s i b i l i t y = 2 0 * \log \left ( \frac { u _ { t o p } ( f ) } { u _ { b o t t o m } ( f ) } \right ),$$

with u bottom and u top defined as the out-of-plane component of displacement extracted at specific points in the model. The numerical study shows that close to the borders, the wavefield experiences boundary effects from the free sides and the bottom plate. Figure. 2b illustrates an example of such existing modes at a frequency of 2900 Hz, which lies within the bandgap. Although no propagation in the structure is expected at this frequency, the displacement at the borders of the modelled prototype does not completely vanish. To mitigate the effects of this boundary noise, u top is selected after the fourth row and not at the last row. Additionally, both the u bottom and u top are averaged over several points (See Figure. 2a).

The numerical transfer function is shown in Figure. 2c. The pink region specifies the bandgap range obtained through eigenfrequency analysis on an infinite model calculated previously. We observe a very good match between the two indicating that 10 by 6 unit cells guarantee the expected functionality. Therefore, we proceed with a prototype of identical size 3D printed via Selective Laser Sintering (SLS). Similar to the numerical case, a thin plate of 30.5 by 3 cm and a thickness of 4 mm was added to the bottom of the design to act as an excitation base. Due to limits in the maximum printable size along the z -direction imposed by the manufacturer, the prototype was assembled out of two smaller samples each consisting of 5 by 6 unit cells (along the y and z directions, respectively). The structure was then attached to a shaker using super glue. The A 0 mode is excited in the structure by sending a sweep signal with the duration of 1.2 seconds carrying frequencies between 1 and 5 kHz as shown in the bottom panels of Figure. 2e. As seen from the experimental setup in Figure. 2e, we measure the out-of-plane velocity of the excited propagating mode using a one-component Scanner Laser Doppler Vibrometer (SLDV) Polytec PSV500. The experimental transmissibility of the prototype structure is determined the same way as in the numerical tests. In this case, we take the ratio of the velocity between the measurements at the top and the bottom of the structure. The results are shown in Figure. 2c.

While both experimental and numerical data on the finite model reveal a bandgap between 1.990 and 3.355 kHz, the infinite model predicts a range between 2.021 and 3.427 kHz. The comparison implies a relative error of 1.5% in the lower frequency and of 2.3% in the upper frequency of the bandgap. We related this small difference to the limited number of unit cells and the presence of material damping when compared to the infinite model. The imaginary part of the wave-number retrieved experimentally (as shown in Figure. 2d) validates the Bragg type bandgap mechanism. It should be mentioned that numerical calculations were initially performed based on the nominal

FIG. 2: (a) Finite plate for frequency domain simulations built by replicating the unit cell in 2D. The green and red points mark the positions where data were extracted to calculate the transmissibility. (b) Mode shape of the finite structure at frequency 2900 Hz. Although the frequency lies within the bandgap, some modes appear in the boundary cells. (c) Transmissibility obtained from the numerical (black) and experimental (blue) data. The pink shaded regions marks the range of the expected bandgap, obtained from the infinite model. (d) Imaginary part of the wave-number based on experimental data. (e) Experimental set-up, along with the signal used as the source wavelet in the experiment.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00590v1-figure-f7b8b724-49e6-4a8f-9b02-4a1996223028.jpg)

values of the elastic parameters that the manufacturer provided in the data sheet for versatile plastic. The nominal value for Young's modulus was E = 1.7 GPa, but after running the experiments, the experimental data revealed a bandgap which was shifted toward lower frequencies compared to the one expected on the basis of the infinite model and nominal values. Numerical testing determined that an E = 1 GPa results in a better match between the experimental and numerical results. Throughout the paper, the numerical results are the ones obtained from the updated value for the Young's modulus. The discrepancy between the nominal value of the Young's modulus and the one of the 3D printed structure has also been reported in other studies [42, 43].

Once the dynamic properties of the Schwarz primitive cell are characterised, we may proceed with designing metastructures for propagation mitigation (metabarriers) and focusing (metalenses) purposes at a given frequency. Both devices are obtained by varying the mass distribution on the unit cells (either by adding or removing the mass) and in a target area on the plate. These variations yield a variation of the dispersion properties.

## IV. METABARRIERS FOR ELASTIC WAVE MITIGATION

The first step of the metabarrier design is to investigate the influence of increasing the mass of the unit cell on the lowest out-of-plane mode. A visual inspection of the bandgap opening/closing mode suggests to locate the mass where the inertia will be maximised, i.e., in the center of the unit cell face. The comparison of dispersion relations and bandgap range between the standard unit cell and the one where a mass density of 8000 kg/m 3 is added to the unit cell (see the insets in Figure. 3a) is shown in Figure. 3b. While the standard unit cell has a bandgap between 2.021 and 3.427 kHz, the cell with the added mass has a bandgap between 1.624 and 2.163 kHz. Thus, we can observe adding the mass to the unit cell shifts the bandgap toward lower frequencies. We leverage this shift and the overlap between the two bandgaps to design the metabarrier. To realize a metabarrier, a 2D plate configuration of N z = 54 by N y = 26 unit cells is designed (See Figure. 3a). Part of this plate is an interior area consisting of 10 by 24 unit cells along the y and z directions (as indicated with the red rectangle in Figure. 3a). Within this area rectangular portions of the unit cells, specified in red in the inset of Figure. 3a, have a higher mass density (8000 kg/m ) than 3 the standard cell.

Performance of the metabarrier is verified by exciting the out-of-plane mode ( x component) with a broadband Ricker wavelet centered at frequency of 5 kHz (See Figure. 3c). The transmissibility is computed before and after the metabarrier (the two points marked by stars in Figure. 3a) to evaluate its effect on the propagating wavefield. Figure. 3d shows that the transmissibility computed at the end of the metabarrier region unfolds a bandgap wider than the one of the standard unit cell. The starting frequency of this widened bandgap coincides with the opening frequency of the bandgap for the unit cell with the added mass (according to Figure. 3b), and the closing frequency matches the one of the bandgap of the standard unit cell.

Moreover, the bandgap computed at the position of the red star shows a stronger attenuation of the A 0 mode (65 dB as opposed to 52 dB). Thus, the metabarrier built by adding the mass to the unit cell leads to a wider and deeper attenuation. The performance of the metabarrier can be further validated by exciting a tapered sinusoid with the frequency of 1.850 kHz. Based on Figure. 3d, this frequency lies in a range where no propagation in the metabarrier region is expected. A snapshot of the results of such a simulation is shown in Figure. 3e. It can be clearly observed that at this frequency, no waves radiate in the metabarrier region and that the waves are reflected backwards. This observation guarantees the possibility of protecting a given region on the plate at a certain frequency range.

Next, we show the design of a metabarrier for the case, where we remove mass from the unit cell, i.e., by increasing the porosity. To do so, spheres with a given radius were subtracted from the unit cell at eight positions and by gradually varying the radius of the sphere, i.e., the resulting increase in porosity, a localized parametric analysis was performed. An example of the unit cell with the removed mass (equivalently, increased porosity) is shown in the inset of Figure. 4c. A comparison between dispersion relations of the bare unit cell and the cell with the increased porosity where the radius of the hole is equal to 5 mm are shown in Figure. 4a. It turns out that increasing the porosity pushes the modes toward lower frequencies, which can be explained by the stronger role of compliance over the mass. This finding enables the design of a metabarrier based on the increased porosity guaranteeing a much lighter metamaterial compared to the added mass instance. A metabarrier with the same size of the plate as the one for the case of the added mass is built and the configuration is illustrated in Figure. 4c. The widening of the bandgap and the higher attenuation within the metabarrier are confirmed in Figure. 4e. The results show the emergence of a bandgap between 1.431 and 3.427 kHz at the end of the metabarrier. Validation of

Х

FIG. 3: (a) Configuration of the metabarrier (indicated by the red square) where the modified cell is the unit cell with the added mass. (b) Comparison between the dispersion relation of the standard unit cell (blue) and the modified unit cell (red) with the added mass density of 8000 kg/m . The bandgap range is marked by shaded areas in each case. (c) Wavelet in time and 3 frequency domain used to compute the transmissibility. (d) Transmissibility measured before (blue) and at the end of the metabarrier (red) at points marked by blue and red stars in panel (a).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00590v1-figure-53e831c6-50d6-48b8-a6ca-86138b5a2fee.jpg)

- (e) Snapshot illustrating propagation of A 0 mode in the plate with an embedded metabarrier.

the metabarrier performance is achieved by exciting the plate with a Hann-windowed sinusoidal at the frequency 1.850 kHz (the wavelet is shown in Figure. 4b). The results confirm the protection of the metabarrier region at this frequency through reflecting the waves back to the bare plate. This result proves the possibility of wave mitigation by further lightening the Schwarz primitive cells.

## V. METADEVICES FOR ELASTIC WAVE FOCUSING

We have shown that by leveraging the variations in the dispersion relation after adding or removing mass to/from the unit cell, a metabarrier can be realized. Next, we leverage this effect

4000

FIG. 4: (a) Comparison between the dispersion relation of the bare unit cell (blue) and the unit cell with the increased porosity of radius 5 mm (red). (b) Sinusoidal source wavelet in time and frequency domain used to validate the performance of the metabarrier. (c) Configuration of the metabarrier where the design parameter is the unit cell with the increased porosity. (d) Snapshot of the A 0 mode of the propagation in the plate with an embedded metabarrier. (e)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00590v1-figure-7c5f0523-7b38-47ba-ad3c-1b3abbf1ee10.jpg)

Transmissibility measured before (blue) and at the end of the metabarrier (red) at points marked by blue and red stars in panel (c).

further in order to focus the incoming waves. Rather than designating an area in the plate with a constant ramp in the unit cell's mass or porosity, we now follow a curvilinear profile defined by an equation. In particular we are interested in the refractive index profile of a GRIN lens such as the Luneburg one [44]. The feasibility of this lens has been previously proved in different types of resonant structures: rods attached to a thin plate [45], octet lattice unit cells [16] in the kHz regime, I-shaped metallic structure embedded in a substrate in GHz regime [46] and also sound focusing based on acoustic metamaterials [47]. However, in these studies the refractive index profile is satisfied by adding mass, increasing the height of rods, thus yielding a heavier zone compared to the surrounding host surface. Here, we tailor the bandgap by increasing the porosity. This is of particular interest, because the refractive index distribution is tuned by the size of the pore on the unit cell resulting in a lightweight meta area.

## A. Design of a Luneburg lens with increased porosity

The Luneburg lens is spherically symmetric and satisfies a parabolic refractive index profile ( n ), where n must vary according to the Equation. 3 as a function of the radial distance:

$$n = \sqrt { 2 - ( \frac { r } { R } ) ^ { 2 } }.$$

Here, r is the radius of each layer and R is the radius of the full lens. In order to satisfy the requirement for the distribution of refractive index, the velocity of the waves need to vary in each spherical layer within the lens region. The refractive index n is defined as:

$$n = \frac { v _ { 0 } } { v }$$

where v 0 is the wave speed at a given frequency in the standard unit cell and v is the wave speed at each layer of the metalens.

Following the procedure outlined above for designing the Luneburg lens, we start by a localized parametric analysis of the unit cell with two different design parameters: namely as the added mass and the increased porosity, and the evolution of the dispersion relation of the A 0 mode in each case is tracked. Values of added mass density are between 5000 and 25000 kg/m 3 and the pore radius for the case of increased porosity varies between 1 to 6 mm, which translates to a porosity range between 90% and 95%. Phase speed is then calculated at the design frequency via v = ω/k , where ω is the angular frequency and k is the wavenumber. Figure. 5a and c show how the frequency range of the A 0 mode changes with increasing mass or porosity. We design the lens to be effective at the frequency of f = 1 kHz. Therefore, the phase speed at this frequency is calculated and its variations with the added mass are shown in Figure. 5b and increasing porosity in Figure. 5d. According to Equation. 3, a Luneburg lens needs a contrast in the refraction index profile that goes from 1 (outermost layer) to 1.4 (central layer). This requires a velocity profile illustrated in Figure. 5e.

By knowing the required velocity profile for the Luneburg lens (from Figure. 5e) and the wave speed variations caused by added mass or porosity (from Figs. 5b and d), we can determine the necessary porosity (or added mass) for each layer of the lens to match the desired refractive index profile, as illustrated in Figure. 5f. It should be noted that a rectangular layer design, as depicted in the inset of Figure. 5h, is preferred due to the cubic cell shape and size, which naturally distributes porosity (or mass) minimising aberrations. Consequently, the radius r in Equation. 3 can be translated into a distance from the center of the lens.

The procedures for both increased porosity and added mass are identical. To simplify and avoid repetition, we present results for the more relevant increased porosity model. The numerical model of the Luneburg lens involves replicating unit cells in two directions within a plate. The lens region spans from -L = -12 ∗ cell size to + L = +12 ∗ cell size , with a cell size of 3 cm. Value of the sphere radius for adding the porosity to each layer of the lens is determined by values reported in Figure. 5f. The out-of-plane mode is excited at one end of the plate with a sinusoidal signal at the design frequency (i.e., f = 1 kHz), resulting in plane wave propagation in the z direction. The source waveform in both the time and frequency domains is depicted in Figure. 5g. A snapshot of the propagating mode in Figure. 5h shows that waves inside the lens area (indicated by the white square) are slowed down and bent, so that they focus at the other end of the lens as expected, marked by the green star. The focusing spot is recovered by extracting the maximum out-of-plane

Lorem ipsum

FIG. 5: (a) localized parametric analysis of the A 0 mode with the added mass. (b) Calculated phase speed based at the design frequency. (c) localized parametric analysis of the A 0 mode with the increased porosity. (d) Calculated phase speed at the design frequency. (e) Velocity profile for Luneburg lens as a function of the distance from the center of the lens. (f) Porosity profile for the Luneburg lens based on the Schwarz primitive cell. (g) Source wavelet in the time and frequency domains. (h) Performance of the lens with square layers at frequency of 1 kHz and the focusing phenomenon. The colored inset shows the pore distribution inside the metalens area.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00590v1-figure-d982e62f-3deb-4e80-9dfe-19ef94254510.jpg)

displacement from -L to + L , as shown in the right panel of the Figure. 5h. Despite the fact that grading is based on a circular shape but applied to squared layers, we observe that focusing at f = 1 kHz is achieved.

## B. Experimental validation of the Luneburg lens with increased porosity

We now show the experimental realisation of the Luneburg lens. The 3D printed prototype shown in Figure. 6c is the top half of the lens validated numerically (i.e., the top half of the white rectangle in Figure. 5h) and therefore is composed of 13 by 25 unit cells. Due to restrictions on the size of the prototype to be 3D printed, the design was divided into three sub samples, which were then 3D printed via the SLS method and glued together. The resulting structure is shown in Figure. 6c. A plate with dimensions of 39 by 3 by 0.4 cm is added to the bottom of the sample and is attached to a shaker to facilitate the excitation. The sample is positioned horizontally on an absorbent sponge in order to minimize the interactions with the ground. Additionally, some playdough has been added to the sides of the plate attached to the shaker to weaken the propagation of unwanted modes of said plate into the structure (See Figure. 6a). Finally, the plate is excited with a sine wave of central frequency 1 kHz (i.e., the design frequency of the lens). An LDV is used to measure the resulting out-of-plane velocity of the propagating wavefield on a grid of 334 points. The resulting wavefield at t = 34 ms is shown in Fig. 6b, where for the sake of a better visualization the grid points are mirrored horizontally in order to mimic the wavefield in a full size lens. The measurement points on the bottom lines of the structure are excluded because as shown in the inset of Figure. 6c, some boundary cells have come out of the 3D printing broken. Bending of the out-of-plane mode within the graded area clearly occurs in Figure. 6b, thus experimentally validating the analysed Luneburg lens design.

## VI. DISCUSSION

Static properties of the Schwarz primitive cell are well known and research is ongoing to design porous and bone mimicking scaffolds based on such geometries for medical applications and bone implants. From a dynamic point of view, the performance of structures based on the Schwarz primitive cell for vibration and elastic wave guiding has been previously investigated [26, 28]. However, studies on the capacity of such structures in regard to wave mitigation and focusing in the context of metamaterials are missing. In this work, we showed that the Schwarz primitive cell can be used as the building block of metadevices for tailoring, mitigation and focusing of elastic waves. This validation paves the way for constructing further metadevices such as GRIN lenses with a hyperbolic secant profile for the refractive index profile, which are effective in a broader range of frequencies and not only at a single frequency. Interestingly, we observed that both adding and removing the mass shift the A 0 mode toward lower frequencies which is an effective means to designing metabarriers and focusing devices. This behaviour is expected for the case of adding mass to the unit cell. However, it seems that by increasing porosity the compliance and lower stiffness have a greater effect than the resulting mass removal, thus, causing the shift of the out-of-plane mode toward lower frequencies. The experimentally demonstrated case of increasing porosity to achieve wavefield focusing has further potential as it results in a much lighter structure compared to previously-realized metadevices in the literature, which are on the basis of adding the mass to the structure or increasing the height of the rod resonators, resulting in a heavier metastructure.. Devices based on the increased porosity are on one hand more feasible for 3D printing and therefore experimental validation. On the other hand, removing the mass changes the static properties of such cells and although the dynamic capacity of the cell based on the increased porosity is verified, research is needed to investigate the static properties of such lightweight structures. The material used in this study (i.e., versatile plastic) makes the structure fragile considering the delicacy of the unit cell, which limits application of metadevices for real world applications problems. In this

FIG. 6: (a) Experimental set-up for wave bending validation in the graded area of the Luneburg lens based on increased porosity. (b) Experimental measurement of the out-of-plane propagating mode at frequency of 1 kHz. The out-of-plane mode bending, marked by red dashed lines, is apparent during the propagation in the graded area. (c) 3D printed prototype of the Luneburg lens. The inset shows the cells on the border of the prototype, which are damaged through the 3D printing process and therefore the measured data at those points are excluded.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00590v1-figure-84c7460c-812d-4006-ad60-a0c0118adcfc.jpg)

regard, optimization of the pore shape and curvature as well as the use of an additional solid matrix and material optimization will be a necessary step to qualify metastructures based on the Schwarz primitive cell for a practical engineering usage.

## VII. CONCLUSION

In this study, we retrieved the band diagram of a Schwarz primitive unit cell made of versatile plastic based on numerical investigations both in the time and frequency domain. We also validated the dispersion relation through laboratory experiments on the small scale printed prototypes using a 1D laser Doppler vibrometer. Both the numerical and experimental analyses revealed a bandgap of Bragg type that matches very well the bandgap expected from the infinite model. Next, the variations of the dispersion relation when adding mass or increasing porosity were leveraged in

order to build metastructures such as metabarriers and GRIN metalenses. Both metabarriers show a wider and deeper attenuation performance, which indicates the capacity of the Schwarz primitive cell for further vibration mitigation. We also show that a Luneburg lens with a parabolic refractive index profile can be constructed for focusing the out-of-plane mode by both adding mass or porosity to the standard cell. Although the Luneburg lens is spherically symmetric and therefore is often designed based on circular layers of varying properties in the meta area, we show that the design still works with rectangular layers. We present a first experimental verification of wavefield focusing by increasing the porosity in an area of a plate constructed of Schwarz primitive unit cells. Our results encourage further application of bio-mimetic structures for tailoring and manipulating of elastic waves. Moreover, the proposed unit cell is scalable and therefore tunable to perform at a given frequency of interest leading to novel functionalities in elastic wave propagation including energy harvesting.

## VIII. ACKNOWLEDGEMENTS

A.H.N and A.C were supported by the H2020 FETOpen project BOHEME under Grant Agreement No. 863179, and A.C, H.R.T and B.Z acknowledge funding from the H2020 FET-proactive project MetaVEH, European Union under the Grant Agreement No. 952039. This work was supported by a grant from the Swiss National Supercomputing Centre (CSCS) under project ID s1185. Laboratory Equipment was provided to A.C. through the SNSF R'Equip grant 205418.

- [1] Yufan Tang, Shuwei Ren, Han Meng, Fengxian Xin, Lixi Huang, Tianning Chen, Chuanzeng Zhang, and Tian Jian Lu. Hybrid acoustic metamaterial as super absorber for broadband low-frequency sound. Scientific Reports , 2017. doi:10.1038/srep43340. URL https://doi.org/10.1038/srep43340 .
- [2] Jeonghoon Park, Dongwoo Lee, and Junsuk Rho. Recent advances in non-traditional elastic wave manipulation by macroscopic artificial structures. Applied Sciences , 10(2), 2020. ISSN 2076-3417. doi:10.3390/app10020547. URL https://www.mdpi.com/2076-3417/10/2/547 .
- [3] Zhongsheng Chen, Bin Guo, Yongmin Yang, and Congcong Cheng. Metamaterials-based enhanced energy harvesting: A review. Physica B: Condensed Matter , 438:1-8, 2014. ISSN 0921-4526. doi:https://doi.org/10.1016/j.physb.2013.12.040. URL https://www.sciencedirect.com/science/ article/pii/S092145261300834X .
- [4] Bao Zhao, Henrik R. Thomsen, Jacopo M. De Ponti, Emanuele Riva, Bart Van Damme, Andrea Bergamini, Eleni Chatzi, and Andrea Colombi. A graded metamaterial for broadband and high-capability piezoelectric energy harvesting. Energy Conversion and Management , 269:116056, 2022. ISSN 0196-8904. doi:https://doi.org/10.1016/j.enconman.2022.116056. URL https://www. sciencedirect.com/science/article/pii/S0196890422008445 .
- [5] Tena Dubˇ cek, Daniel Moreno-Garcia, Thomas Haag, Parisa Omidvar, Henrik R. Thomsen, Theodor S. Becker, Lars Gebraad, Christoph B¨rlocher, Fredrik Andersson, Sebastian D. Huber, Dirk-Jan van Maa nen, Luis Guillermo Villanueva, Johan O. A. Robertsson, and Marc Serra-Garcia. Binary classification of spoken words with passive phononic metamaterials, 2023.
- [6] Marco Miniaci, Anastasiia Krushynska, Federico Bosia, and Nicola M Pugno. Large scale mechanical metamaterials as seismic shields. New Journal of Physics , 18(8):083041, aug 2016. doi:10.1088/13672630/18/8/083041. URL https://dx.doi.org/10.1088/1367-2630/18/8/083041 .
- [7] Bao Zhao, Henrik R Thomsen, Xingbo Pu, Shitong Fang, Zhihui Lai, Bart Van Damme, Andrea Bergamini, Eleni Chatzi, and Andrea Colombi. A nonlinear damped metamaterial: Wideband attenuation with nonlinear bandgap and modal dissipation. arXiv preprint arXiv:2307.14165 , 2023.

- [8] Kyriakos Alexandros Chondrogiannis, Andrea Colombi, Vasilis Dertimanis, and Eleni Chatzi. Computational verification and experimental validation of the vibration-attenuation properties of a geometrically nonlinear metamaterial design. Phys. Rev. Appl. , 17:054023, May 2022. doi: 10.1103/PhysRevApplied.17.054023. URL https://link.aps.org/doi/10.1103/PhysRevApplied. 17.054023 .
- [9] L. Brillouin. Wave Propagation in Periodic Structures: Electric Filters and Crystal Lattices . Dover phoenix editions. Dover Publications, 2003. ISBN 9780486495569. URL https://books.google.ch/ books?id=m2WmGiU5nUwC .
- [10] C. Kittel. Introduction to Solid State Physics . Wiley, 2004. ISBN 9780471415268. URL https: //books.google.ch/books?id=kym4QgAACAAJ .
- [11] Andrea Colombi, Philippe Roux, Sebastien Guenneau, Philippe Gueguen, and Richard V Craster. Forests as a natural seismic metamaterial: Rayleigh wave bandgaps induced by local resonances. Scientific Reports , 6(1):1-7, 2016.
- [12] Xiaoming Zhou, Xiaoning Liu, and Gengkai Hu. Elastic metamaterials with local resonances: an overview. Theoretical and Applied Mechanics Letters , 2(4):041001, 2012. ISSN 2095-0349. doi: https://doi.org/10.1063/2.1204101. URL https://www.sciencedirect.com/science/article/pii/ S2095034915301586 .
- [13] Mattia Cenedese, Edoardo Belloni, and Francesco Braghin. Interaction of Bragg scattering bandgaps and local resonators in mono-coupled periodic structures. Journal of Applied Physics , 129(12), 03 2021. ISSN 0021-8979. doi:10.1063/5.0038438. URL https://doi.org/10.1063/5.0038438 . 124501.
- [14] A.O. Krushynska, M. Miniaci, F. Bosia, and N.M. Pugno. Coupling local resonance with bragg band gaps in single-phase mechanical metamaterials. Extreme Mechanics Letters , 12:30-36, 2017. ISSN 2352-4316. doi:https://doi.org/10.1016/j.eml.2016.10.004. URL https://www.sciencedirect.com/ science/article/pii/S2352431616300542 . Frontiers in Mechanical Metamaterials.
- [15] Taehwa Lee and Hideo Iizuka. Bragg scattering based acoustic topological transition controlled by local resonance. Phys. Rev. B , 99:064305, Feb 2019. doi:10.1103/PhysRevB.99.064305. URL https: //link.aps.org/doi/10.1103/PhysRevB.99.064305 .
- [16] Giulia Aguzzi, Constantinos Kanellopoulos, Richard Wiltshaw, Richard V Craster, Eleni N Chatzi, and Andrea Colombi. Octet lattice-based plate for elastic wave control. Scientific Reports , 12(1):1-14, 2022.
- [17] Matthieu Rupin, Fabrice Lemoult, Geoffroy Lerosey, and Philippe Roux. Experimental demonstration of ordered and disordered multiresonant metamaterials for lamb waves. Physical review letters , 112 (23):234301, 2014.
- [18] Srijan Datta, Antonello Tamburrino, and Lalita Udpa. Gradient index metasurface lens for microwave imaging. Sensors , 22(21), 2022. ISSN 1424-8220. doi:10.3390/s22218319. URL https://www.mdpi. com/1424-8220/22/21/8319 .
- [19] Rafael Fuentes-Dominguez, Mengting Yao, Andrea Colombi, Paul Dryburgh, Don Pieris, Alexander Jackson-Crisp, Daniel Colquitt, Adam Clare, Richard J. Smith, and Matt Clark. Design of a resonant luneburg lens for surface acoustic waves. Ultrasonics , 111:106306, 2021. ISSN 0041-624X. doi:https://doi.org/10.1016/j.ultras.2020.106306. URL https://www.sciencedirect.com/science/ article/pii/S0041624X20302419 .
- [20] Alfonso Climente, Daniel Torrent, and Jos´ S´nchez-Dehesa. e a Gradient index lenses for flexural waves based on thickness variations. Applied Physics Letters , 105(6):064101, 2014.
- [21] Hermann Amandus Schwarz. Gesammelte mathematische abhandlungen , volume 260. American Mathematical Soc., 1972.
- [22] H Montazerian, E Davoodi, M Asadi-Eydivand, J Kadkhodapour, and M Solati-Hashjin. Porous scaffold internal architecture design based on minimal surfaces: A compromise between permeability and elastic properties. Materials &amp; Design , 126:98-114, 2017. ISSN 0264-1275. doi:https://doi.org/10.1016/j.matdes.2017.04.009. URL https://www.sciencedirect.com/science/ article/pii/S0264127517303593 .
- [23] Sanjairaj Vijayavenkataraman, Lei Zhang, Shuo Zhang, Jerry Ying Hsi Fuh, and Wen Feng Lu. Triply periodic minimal surfaces sheet scaffolds for tissue engineering applications: An optimization approach toward biomimetic scaffold design. ACS Appl. Bio Mater , 1, 2018. URL https://doi.org/10.1021/

acsabm.8b00052 .

- [24] Chen Pan, Yafeng Han, and Jiping Lu. Design and optimization of lattice structures: A review. Applied Sciences , 10(18), 2020. ISSN 2076-3417. doi:10.3390/app10186374. URL https://www.mdpi. com/2076-3417/10/18/6374 .
- [25] Baris Sokollu, Orhan Gulcan, and Erhan Ilhan Konukseven. Mechanical properties comparison of strut-based and triply periodic minimal surface lattice structures produced by electron beam melting. Additive Manufacturing , 60:103199, 2022. ISSN 2214-8604. doi: https://doi.org/10.1016/j.addma.2022.103199. URL https://www.sciencedirect.com/science/ article/pii/S2214860422005887 .
- [26] N.V Viet and W Zaki. Free vibration and buckling characteristics of functionally graded beams with triply periodic minimal surface architecture. Composite Structures , 274:114342, 2021. ISSN 02638223. doi:https://doi.org/10.1016/j.compstruct.2021.114342. URL https://www.sciencedirect.com/ science/article/pii/S0263822321008047 .
- [27] Yuning Guo, Matheus Inguaggiato Nora Rosa, Mohit Gupta, Benjamin Emerson Dolan, Brandon Fields, Lorenzo Valdevit, and Massimo Ruzzene. Minimal surface-based materials for topological elastic wave guiding. Advanced Functional Materials , 32(30):2204122, 2022.
- [28] Yongbin Ma and Qingfeng Cheng. Semi-analytical solutions for wave propagation of periodically repetitive schwarz primitive triply periodic minimal surface based structure embedded with acoustic black hole resonators. International Journal of Acoustics &amp; Vibration , 28(3), 2023.
- [29] Silva Tarcisio, Lu Jin-You, Abu Al-Rub Rashid K., and Lee Dong-Wook. Investigation on tailoring the width and central frequency of bandgaps of tpms structures. International Journal of Mechanics and Materials in Design , 2023. ISSN 1573-8841. doi:10.1007/s10999-023-09677-2. URL https://doi. org/10.1007/s10999-023-09677-2 .
- [30] Lu Jin-You, Silva Tarcisio, Alzaabi Fatima, Abu Al-Rub Rashid K, and Lee Dong-Wook. Insights into acoustic properties of seven selected triply periodic minimal surfaces-based structures: A numerical study. Journal of Low Frequency Noise, Vibration and Active Control , 43, 2023. doi:doi: 10.1177/14613484231190986. URL https://doi.org/10.1177/14613484231190986 .
- [31] Waiel Elmadih, Wahyudin P. Syam, Ian Maskery, Dimitrios Chronopoulos, and Richard Leach. Mechanical vibration bandgaps in surface-based lattices. Additive Manufacturing , 25:421-429, 2019. ISSN 2214-8604. doi:https://doi.org/10.1016/j.addma.2018.11.011. URL https://www.sciencedirect.com/ science/article/pii/S2214860418302781 .
- [32] Min Lai, Alex N. Kulak, Daniel Law, Zhibing Zhang, Fiona C. Meldrum, and D. Jason Riley. Profiting from nature: macroporous copper with superior mechanical properties. Chem. Commun. , pages 35473549, 2007. doi:10.1039/B707469G. URL http://dx.doi.org/10.1039/B707469G .
- [33] Oraib Al-Ketan and Rashid K. Abu Al-Rub. Multifunctional mechanical metamaterials based on triply periodic minimal surface lattices. Advanced Engineering Materials , 21(10):1900524, 2019. doi:https://doi.org/10.1002/adem.201900524. URL https://onlinelibrary.wiley.com/doi/abs/10. 1002/adem.201900524 .
- [34] Paul J.F. Gandy and Jacek Klinowski. Exact computation of the triply periodic schwarz p minimal surface. Chemical Physics Letters , 322(6):579-586, 2000. ISSN 0009-2614. doi: https://doi.org/10.1016/S0009-2614(00)00453-X. URL https://www.sciencedirect.com/science/ article/pii/S000926140000453X .
- [35] Diab W. Abueidda, Iwona Jasiuk, and Nahil A. Sobh. Acoustic band gaps and elastic stiffness of pmma cellular solids based on triply periodic minimal surfaces. Materials &amp; Design , 145:2027, 2018. ISSN 0264-1275. doi:https://doi.org/10.1016/j.matdes.2018.02.032. URL https://www. sciencedirect.com/science/article/pii/S0264127518301163 .
- [36] Daniel Royer and Eug` ene Dieulesaint. Elastic waves in solids. 1980.
- [37] COMSOL Multiphysics ® v. 5.4. Structural Mechanics Module User's Guide .
- [38] D. Alleyne and P. Cawley. A two-dimensional Fourier transform method for the measurement of propagating multimode signals. The Journal of the Acoustical Society of America , 89(3):1159-1168, 03 1991. ISSN 0001-4966. doi:10.1121/1.400530. URL https://doi.org/10.1121/1.400530 .

- [39] Felix Bloch. Uber die quantenmechanik der elektronen in kristallgittern. ¨ Zeitschrift f¨r physik u , 52(7-8): 555-600, 1929.
- [40] Dimitri Komatitsch and Jeroen Tromp. Introduction to the spectral element method for threedimensional seismic wave propagation. Geophysical Journal International , 139(3):806-822, 12 1999. ISSN 0956-540X. doi:10.1046/j.1365-246x.1999.00967.x. URL https://doi.org/10.1046/j. 1365-246x.1999.00967.x .
- [41] Bart Van Damme and Armin Zemp. Measuring dispersion curves for bending waves in beams: A comparison of spatial fourier transform and inhomogeneous wave correlation. Acta Acustica united with Acustica , 104:228-234, 03 2018. doi:10.3813/AAA.919164.
- [42] Giulia Aguzzi, Henrik R. Thomsen, Aida Hejazi Nooghabi, Richard Wiltshaw, Richard V. Craster, Eleni N. Chatzi, and Andrea Colombi. Architected frames for elastic wave attenuation: Experimental validation and local tuning via affine transformation. Applied Physics Letters , 121(20), 11 2022. ISSN 0003-6951. doi:10.1063/5.0119903. URL https://doi.org/10.1063/5.0119903 . 201702.
- [43] Danilo Beli, Adriano Fabro, Massimo Ruzzene, and J.R.F. Arruda. Wave attenuation and trapping in 3d printed cantilever-in-mass metamaterials with spatially correlated variability. Scientific Reports , 9, 04 2019. doi:10.1038/s41598-019-41999-0.
- [44] R.K. Luneburg, M. Herzberger, and Brown University. Graduage School. Mathematical Theory of Optics: Supplementary Notes by M. Herzberger . Number Bd. 1 in Mathematical Theory of Optics: Supplementary Notes by M. Herzberger. 1944. URL https://books.google.ch/books?id=DAsJAQAAMAAJ .
- [45] A. Colombi. Resonant metalenses for flexural waves. J. Acoust. Soc. Am. , 140(5):EL423, 2016.
- [46] Haibing Chen, Qiang Cheng, Aihua Huang, Junyan Dai, Huiying Lu, Jie Zhao, Huifeng Ma, Weixiang Jiang, and Tiejun Cui. Modified luneburg lens based on metamaterials. International Journal of Antennas and Propagation , 2015, 2015.
- [47] Baoguo Yuan, Jiyu Liu, Houyou Long, Ying Cheng, and Xiaojun Liu. Sound focusing by a broadband acoustic Luneburg lens. The Journal of the Acoustical Society of America , 151(3):2238-2244, 03 2022. ISSN 0001-4966. doi:10.1121/10.0009909. URL https://doi.org/10.1121/10.0009909 .
Aida Hejazi Nooghabi Henrik R. Thomsen Bao Zhao Andrea Colombi
Elastic wave control in reticulated plates using Schwarz primitive cells
In this work, the Schwarz primitive unit cell is used as the building block
of different types of metastructures for steering and focusing elastic
vibrations. The emergence of a Bragg-type bandgap when constructing a
two-dimensional plate from such unit cells is experimentally validated. It is
demonstrated that increasing both mass and porosity of the Schwarz primitive
leads to a decrease in the frequency of the out-of-plane propagating wave
targeted in this study. By arranging these modified Schwarz primitive unit
cells in constant and graded layouts, two-dimensional plates with an embedded
metabarrier and a metalens are numerically designed. The metabarrier protects
an interior area of the plate from the propagating waves on a wide frequency
band (approx. 1.4-3.4 kHz). Equally, the refractive index profile necessary for
gradient index lenses is obtained via a progressive variation of the added mass
or, alternatively, the porosity of the unit cell over a rectangular area. For
the first time, bending of the out-of-plane mode towards the focusing point is
practically validated in a challenging mesoscale experiment requiring the
assembly of different three-dimensional printed sections of the plate. The
increased porosity design is advantageous not only in terms of overall
lightweight, but also towards additive manufacturing as it requires less
material.