Next Article in Journal
Multi-Frequency GPR Data Fusion through a Joint Sliding Window and Wavelet Transform-Weighting Method for Top-Coal Structure Detection
Next Article in Special Issue
A Single-Phase Lightweight Double-Leaf Multi-Stage Acoustic Black Hole Model of Metamaterial
Previous Article in Journal
Investigating the Out-of-Plane Bending Stiffness Properties in Hybrid Species Diagonal-Cross-Laminated Timber Panels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Active Tunable Elastic Metasurface for Abnormal Flexural Wave Transmission

1
School of Mechanical and Electrical Engineering College, Hainan University, Haikou 570228, China
2
College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150006, China
3
Qingdao Innovation and Development Center, Harbin Engineering University, Qingdao 266000, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(7), 2717; https://doi.org/10.3390/app14072717
Submission received: 11 January 2024 / Revised: 10 March 2024 / Accepted: 18 March 2024 / Published: 24 March 2024

Abstract

:
An active elastic metasurface has more flexibility than a passively modulated elastic metasurface, owing to the manipulation of the phase gradient that can be realized without changing the geometrical configuration. In this study, a negative proportional feedback control system was employed to provide positive active control stiffness for adaptive unit cells, with the aim of achieving the active modulation of the phase gradient. The relationship between the control gain and the phase velocity of the flexural wave was derived, and the transfer coefficients and phase shifts of the flexural wave through the adaptive unit cells were resolved using the transfer matrix method. Finite element simulations for wave propagations in the adaptive unit cells were conducted, and they verified the analytic solutions. Based on this theoretical and numerical work, we designed active elastic metasurfaces with adaptive unit cells with sub-wavelength thicknesses according to the generalized Snell’s law. These metasurfaces show flexibility in achieving abnormal functions for transmitted waves, including negative refraction and wave focusing, and transforming guided waves at different operating frequencies by manipulating the control gain. Therefore, the proposed active metasurface has great potential in the fields of the tunable manipulation of elastic waves and the design of smart devices.

1. Introduction

Metasurfaces [1,2,3,4,5] are artificially constructed materials with subwavelength thicknesses in the direction of wave propagation, exhibiting an excellent ability to manipulate wavefronts, as well as having a more compact physical space and simplicity compared to conventional metamaterials [6,7,8,9,10]. Metasurfaces have been widely used in various fields, such as medical imaging, holographic imaging, and signal processing, following studies on their extraordinary beam modulation ability in the fields of electromagnetics [11,12,13] and acoustics [14,15,16,17,18,19]. Recently, the research on metasurfaces has been expanded to the field of elastic waves [20,21,22,23,24,25,26,27,28]. However, elastic waves possess more degrees of freedom during solid propagation [29,30,31], which leads to a great challenge for researchers. Determining a method to extract and control the propagation of specific modes of guided waves from multimodal guided waves in a solid has become a hot research topic.
Researchers have designed elastic wave metasurfaces based on the generalized Snell’s law (GSL) for wavefront modulation behaviors such as abnormal refraction, wave focusing, and Bessel waves. Zhu et al. [21] first proposed metasurfaces for the abnormal deflection of A0 Lamb waves. Lee et al. [26] utilized the relationship between the effective mass and stiffness to achieve the modulation of transmitted traveling waves. Classical zigzag topology was used to alter the propagation path of flexural waves, taking advantage of the different propagation paths of the flexural wave to realize the phase shift of the basic unit cell across the full 2π [23]. A series of abnormal modulation behaviors of the flexural wave was realized by adjusting the phase gradient. A structure for the passive modulation of multi-mode guided waves was achieved by adding different lengths of cover layers on the unit element, and the theoretical analysis was verified by finite element simulations and experiments [25].
The above case can be identified as the passive modulation [32,33,34,35,36] of elastic waves. Passive modulation requires a change in the phase gradient by changing the geometry after fixing the phase distribution, which makes metasurfaces much less practical. Therefore, adaptive elastic metasurfaces have attracted considerable interest. Si-Min Yuan et al. [37] proposed a nuts-and-screws structure where the tunability of the phase shift is induced by the screw-in depth based on the dispersion theory of locally resonant phononic crystals. By adjusting the position of the nut on the screw to achieve the full 2π coverage of the phase shift of the functional element, the abnormal refraction and wave-focusing functions of the broadband plate wave were achieved based on the GSL without re-fabricating the metasurface. There are various ways to use piezoelectric materials for active modulation. Piezoelectric patches can be stacked together to form a unit cell, which can be adjusted via negative capacitance circuits connected to achieve a wide frequency range of elastic longitudinal waves with a variety of modulation functions. It is also possible to attach a piezoelectric sheet to the substrate material and use electromechanical tuning to modulate the propagation of asymmetric mode Lamb waves [38]. Shixuan Shao et al. [39] presented a novel metasurface that is proposed to simultaneously manipulate multi-mode guided waves in the plate, including shear horizontal waves, symmetric mode, and anti-symmetric mode Lamb waves. Monolithic in-plane and out-of-plane polarized piezoelectric metasurface sheets were staggered on both faces of the substrate, and each piezoelectric patch was individually connected to a negative capacitance circuit. The metasurface can be multifunctionally modulated in each mode by adjusting the shunt’s negative capacitance.
Currently, most metasurfaces are actively modulated by negative-capacitance circuits, but negative-capacitance circuits are difficult to achieve in real-word applications. In order to satisfy the requirement of modulating wave transmission in real time, the piezoelectric patches were attached at the upper and lower surfaces of the unit cells, with the upper surface’s piezoelectric patch acting as a brake, the lower surface piezoelectric patch acting as a sensor, and the piezoelectric patch externally connected to a negative proportional feedback control strategy [40]. This strategy was employed to provide positive active control stiffness to the piezoelectric sensor/actuator patch in the present work. This method can be utilized to implement anomalous wave manipulations in an active manner [37,41,42] for the flexural wave. In such systems, the key is to exploit the relationship between flexural wave phase velocity and stiffness by controlling the equivalent stiffness to determine the phase velocity. We derived the relationship between the equivalent stiffness and the phase velocity of the flexural wave and analyzed the transmission coefficients and phase shifts of each cell using the transfer matrix method. Based on the GSL, the adaptive unit cells were cleverly arranged to construct appropriate phase gradients, thus enabling the abnormal refraction [43,44], focusing [45,46,47,48], and transforming the guided wave of the metasurface to be realized at different frequencies.
The paper is organized as follows. Section 2 outlines the theoretical derivation of the adaptive unit cell and metasurface design. Section 3 demonstrates the abnormal modulation of flexural waves via metasurfaces through finite element simulations. Section 4 summarizes some conclusions.

2. Mechanisms of Actively Tunable Elastic Metasurface

2.1. Description of Adaptive Unit Cells

We used a negative proportional feedback control strategy to modulate the phase of the adaptive unit cells, which can control the flexural wave more efficiently and contribute to the fabrication of the metasurface. The proposed adaptive unit cell is shown in Figure 1, where Figure 1a shows the three-dimensional geometric model of the adaptive unit cell and Figure 1b shows the two-dimensional model. From Figure 1a, it can be seen that the adaptive unit cell contains two sub-cells. The gray and yellow regions indicate the base beam (aluminum) and the piezoelectric patch, respectively; the upper part of the piezoelectric patch is used as the piezoelectric sensor, and the lower part is applied as the piezoelectric actuator. Moreover, the negative proportional feedback control system is connected externally to the piezoelectric patch. The ratio of the sensed voltage generated by the deformation of the sensor to the external voltage applied by the control system is used to realize the positive active stiffness gain [40]. In the adaptive unit cell of length L = l1 + l2, the lengths of sub-cell 1 (l1) and sub-cell 2 (l2) are adjustable, and we define l1/L as p. The thickness of the base beam in sub-cell 1 is hs1 and that in sub-cell 2 is hs2, and the two substrates are covered with piezoelectric patches with a thickness of hp. The width of the piezoelectric patch is the same as the width of the base beam. The total thickness of sub-cell 1 is h1 = hs1 + 2hp and the total thickness of sub-cell 2 is h2 = hs2 + 2hp.

2.2. Derivations of the Governing Equation

The constitutive equation of piezoelectric materials can be expressed as follows:
S i = s ij T j - d in E n D m = d mj T j + ε mn E n ( i , j = 1 , 2 , 3 , 4 , 5 , 6 ; m , n = 1 , 2 , 3 )
where S is mechanical strain, T is mechanical stress, D is electrical displacement, E is the electric field, sij is the flexibility coefficient under constant electric field conditions, d is the piezoelectric stress constant, and ɛmn is the dielectric constant. The correspondence between the numerical subscripts and the coordinate axes is 1→x, 2→y, 3→z, 4→yz, 5→xz, 6→xy.
In the piezoelectric sheet used in this study, the polarization surface is perpendicular to the z-axis, and the action of the electric field in the z-axis direction is mainly considered, ignoring the electric field in the x-axis and y-axis directions, which implies that E1 = E2 = 0. The piezoelectric sheet is constrained in the x-y plane, with the stress relation T3 = T4 = T5 = 0. The piezoelectric material is transverse isotropic in the x-y plane, and it can be shown that s11 = s22.
Consider a piezoelectric material; its constitutive equations can be expressed as [42]
S 1 = s 11 T 1 - d 31 E 3 D 3 = d 31 T 1 + ε 33 E 3
It is worth noting that the potential D3 is constant at the electrode with area As [49,50]. According to the classical laminated beam theory, the equivalent density ρa and the equivalent stiffness Da of the laminated composite beam with the piezoelectric material can be written as [40]
ρ a = ρ s h s + 2 ρ p h p
D a = E s h s 3 12 ( 1 υ s 2 ) + 2 c 11 ( A s I 1 2 + w 1 h p 3 12 )
where ρs, Es, and υs are the mass density, Young’s modulus, and Poisson’s ratio of the substrate material, respectively; ρp is the density of the piezoelectric material; w1 is the width of the adaptive unit cell; hs is the thickness of the base beam; and I1 = (hp + hs)/2.
Following Euler’s beam theory, the vibration equation of a composite beam with a piezoelectric actuator and sensors is
D a 4 w x , t x 4 + ρ a h 2 w x , t t 2 + d 31 I 1 h p 2 Y a x 2 d 31 I 1 h p 2 Y s x 2 = 0
where Ya is the external voltage applied by the control system and Ys is the sensed voltage generated by the deformation of the sensor.
The sensed voltage of the piezoelectric sensor Ys produced by the deformation of the composite base beam can be given by [50]
Y s = d 31 I 1 ε 33 w 1 2 w x 2
where ɛ33 is the dielectric constant of the piezoelectric material and d31 is the piezoelectric stress constant of the piezoelectric material. The metasurface adaptive unit cells were designed using a negative proportional feedback control strategy. The sensing voltage Ys caused by the deformation of the piezoelectric sensor is measured and fed back to the piezoelectric actuator as an external control voltage. Therefore, the relationship between the external control voltage Ya and the sensing voltage Ys can be expressed as
Y a = g Y s = g d 31 I 1 ε 33 w 1 2 w x 2
where g is the proportional feedback control gain of the piezoelectric actuator. The active stiffness can be generated by an external control voltage applied to the piezoelectric actuator, which can manipulate the flexural wave propagation characteristics of the adaptive unit cell. The variation in g is also limited to 0~100 in order to maintain the stability of the constant voltage electric field.
Substituting Equations (6) and (7) into Equation (5), the coupled equation can be obtained as
D a 1 4 w ( x , t ) x 4 + ρ a h 2 w ( x , t ) t 2 = 0
where Da1 is the equivalent stiffness of the beam with the piezoelectric actuator and sensor acting together and can be written as
D a 1 = D a + ( 1 + g ) d 31 2 I 1 2 ε 33 w 1 h p
The separable solution to Equation (8) can be written as w(x,t) = w(x)eiωt, where ω is 2πf. Substituting this solution and Equation (9) into Equation (8), the control equation for the flexural wave can be obtained as
4 w x x 4 k 4 2 w x t 2 = 0
where k is the wavenumber and k4 is expressed as
k 4 = ρ a h ω 2 D a 1
The phase velocity c for the flexural wave is
c = D a 1 ω 2 ρ a h
Equation (12) indicates that flexural waves are dispersive, and their phase velocity depends on the flexural stiffness. In designing the metasurface, we can vary the magnitude of the gain using a negative proportional feedback control system to obtain the desired phase velocity.

2.3. Generalized Snell’s Law

According to the generalized Snell’s law, the relationship between the incident and transmitted angles of a flexural wave passing through the metasurface is as follows
sin ( θ t ) λ t - sin ( θ i ) λ i = 1 2 π d φ y d y
where λi = c/f, θi and θt are the incident and transmitted angles, respectively; λi and λt are the wavelength of the incident and transmitted flexural waves, respectively; and y/dy is the phase gradient along the y-axis.
The wavelengths of the incident and transmitted regions are the same when the same material is applied, and hence λi = λt. Equation (13) can be organized as
sin ( θ t ) - sin ( θ i ) λ i = 1 2 π d φ y d y
In particular, if the wave is incident vertically, the magnitude of the angle of refraction is calculated according to the following formula:
θ t = arcsin λ t 2 π × d φ y d y

2.4. Focusing Principle

For a certain focal position (x0, y0), the phase distribution along the metasurface is
φ y = k i ( ( y y 0 ) 2 + x 0 2 x 0 ) + φ y 0
where y0 is the transverse coordinate of the focal point O, x0 is the distance from the focal point O, φy is the phase, and ki is the wavenumber in the incident region plate. When designing the focal point y0 = 0, the focusing equation can be simplified as
φ y = k i ( y 2 + x 0 2 x 0 )

2.5. Theoretical Formulations of the Transmission Matrix Method

In this subsection, we will use the transmission matrix method (TMM) to resolve the transmission coefficients and phase shifts of the unit cell. It has been shown that the transmission coefficient and phase shift of a three-dimensional model is equivalent to those of a two-dimensional model [44]. Therefore, the wave propagation can be considered as a plane strain problem in the xz plane, which is independent of the y-axis direction, as shown by the basic unit in Figure 2. The general solution to Equation (10) is
w x = a + e i k x + a e i k x + b + e k x + b e k x
where a+e−ikx and b+ekx are positive propagating and attenuating waves, respectively, and a_eikx and b_ekx are negative propagating and attenuating waves, respectively. a+, a_, b+, and b_ represent the complex coefficients. The displacement w(x), slope θ(x), bending moment M(x), and shear force V(x) compose the state vector V = {w, θ, M, V}T.
The basic unit in Figure 2 can be divided into four regions, where region I is the incident region; regions II and III are sub-cell 1 and sub-cell 2, respectively; and region IV is the transmission region. Each region shares the same central axis. The materials of region I and IV are the same as the base beam material with thickness h0 = h1. As shown in Figure 2, the state vectors on the left and right sides of the connection between the regions working together can be defined as VR1, VL2, VR2, VL3, and VR3 and VL4. The connectivity at the boundaries can be determined from the fact that VR1 = VL2, VR2 = VL3 and VR3 = VL4.
In this way, VL2 can be expressed as
V L 2 = 1 1 1 1 i k 2 i k 2 i k 2 k 2 D 2 k 2 2 D 2 k 2 2 D 2 k 2 2 D 2 k 2 2 D 2 k 2 3 i D 2 k 2 3 D 2 k 2 3 D 2 k 2 3 a + 2 a 2 b + 2 b 2
where k2 is the wavenumber in region II and D2 is the flexural rigidity of region II. According to the above equation, VL2 can be expressed as VL2 = N2A, which can be obtained as A = N 2 1 VL2.
Similarly, VR2 can be written as
V R 2 = e i k 2 l 1 e i k 2 l 1 e k 2 l 1 e k 2 l 1 i k 2 e i k 2 l 1 i k 2 e i k 2 l 1 k 2 e k 2 l 1 k 2 e k 2 l 1 D 2 k 2 2 e i k 2 l 1 D 2 k 2 2 e i k 2 l 1 D 2 k 2 2 e k 2 l 1 D 2 k 2 2 e k 2 l 1 i D 2 k 2 3 e i k 2 l 1 i D 2 k 2 3 e i k 2 l 1 D 2 k 2 3 e k 2 l 1 D 2 k 2 3 e k 2 l 1 a + 2 a 2 b + 2 b 2
For brevity, VR2 can be expressed as VR2 = M2A. Combining Equations (19) and (20) leads to VR2 = M2  N 2 1 VL2 = U2 VL2.
With respect to VR3, there is the following mathematical relationship.
V L 3 = e i k 3 l 1 e i k 3 l 1 e k 3 l 1 e k 3 l 1 i k 3 e i k 3 l 1 i k 3 e i k 3 l 1 k 3 e k 3 l 1 k 3 e k 3 l 1 D 3 k 3 2 e i k 3 l 1 D 3 k 3 2 e i k 3 l 1 D 3 k 3 2 e k 3 l 1 D 3 k 3 2 e k 3 l 1 i D 3 k 3 3 e i k 3 l 1 i D 3 k 3 3 e i k 3 l 1 D 3 k 3 3 e k 3 l 1 D 3 k 3 3 e k 3 l 1 a + 3 a 3 b + 3 b 3
where k3 is the wavenumber in region III and D3 is the flexural rigidity of region III. According to the above equation, VL3 can be expressed as VL2 = N3A, which can be obtained as A = N 3 1 VL3.
The state vector VR3 can be expressed by the following equation:
V R 3 = e i k 3 ( l 1 + l 2 ) e i k 3 ( l 1 + l 2 ) e k 3 ( l 1 + l 2 ) e k 3 ( l 1 + l 2 ) i k 3 e i k 3 ( l 1 + l 2 ) i k 3 e i k 3 ( l 1 + l 2 ) k 3 e k 3 ( l 1 + l 2 ) k 3 e k 3 ( l 1 + l 2 ) D 3 k 3 2 e i k 3 ( l 1 + l 2 ) D 3 k 3 2 e i k 3 ( l 1 + l 2 ) D 3 k 3 2 e k 3 ( l 1 + l 2 ) D 3 k 3 2 e k 3 ( l 1 + l 2 ) i D 3 k 3 3 e i k 3 ( l 1 + l 2 ) i D 3 k 3 3 e i k 3 ( l 1 + l 2 ) D 3 k 3 3 e k 3 ( l 1 + l 2 ) D 3 k 3 3 e k 3 ( l 1 + l 2 ) a + 3 a 3 b + 3 b 3
According to the above equation, VR3 can be expressed as VR3 = M3A, and from Equations (21) and (22), it can be expressed as VR3 = M3 N 3 1 VL2= U3 VL2.
The above equation gives the following relationship
V L 4 = V R 4 = U 3 V L 3 = U 3 V R 2 = U 3 U 2 V L 2 = U 3 U 2 V R 1
It is possible to solve for the transmission coefficients and phase shifts of the flexural waves from Equation (23). For this purpose, the wavefields in regions I, Ⅱ, III and IV are prepared as follows:
w 1 x = e i k 1 x + r e i k 1 x + r e i k 1 x w 2 x = a + 2 e i k 2 x + a 2 e i k 2 x + b + 2 e k 2 x + b 2 e k 2 x w 3 x = a + 3 e i k 3 x + a 3 e i k 3 x + b + 3 e k 3 x + b 3 e k 3 x w 4 x = t e i k 4 x + t e k 4 x
where r, r*, t, and t* are the amplitude ratios of the reflected propagating wave, the reflected attenuating wave, the transmitted propagating wave, and the transmitted attenuating wave to the incident wave, respectively. k4 is the wavenumber in region IV, and k1 is the wavenumber in region I.
The equation of the VR1 relationship is as follows:
V R 1 = 1 1 i k 1 k 1 D 1 k 1 2 D 1 k 1 2 i D 1 k 1 3 D 1 k 1 3 r r + 1 i k 1 D 1 k 1 2 i D 1 k 1 3
where D1 is the flexural rigidity of region I. According to the above equation, the mathematical relationship of VR1 can be written as VR1 = U1r + H1.
VL4 can be expressed as
V L 4 = e i k 4 ( l 1 + l 2 ) e k 4 ( l 1 + l 2 ) i k 4 e i k 4 ( l 1 + l 2 ) k 4 e k 4 ( l 1 + l 2 ) D 4 k 4 2 e i k 4 ( l 1 + l 2 ) D 4 k 4 2 e k 4 ( l 1 + l 2 ) i D 4 k 4 3 e i k 4 ( l 1 + l 2 ) D 4 k 4 3 e k 4 ( l 1 + l 2 ) t t
where D4 is the flexural rigidity of region IV. According to the above equation, VL4 can be expressed as VL4 = U4t.
Combining Equations (23) and (25) leads to
U 3 U 2 U 1 r + U 3 U 2 H 1 = U 4 t
According to the above equation, the following relationship can be obtained:
U 4 U 3 U 2 U 1 t r = U 3 U 2 H 1
Based on this equation, the solutions for r, r*, t, and t* can be solved computationally. Since the same bending stiffness is present in regions I and IV, |t|2 is taken to calculate the transmission coefficient. The phase can be calculated with the following equation [8].
φ = arctan ( Im ( t ) Re ( t ) ) + π 2 arctan ( Im ( t ) Re ( t ) ) + 3 π 2   Re t < 0
where “Im” and “Re” are the functions to extract the imaginary and real parts, respectively.

3. Numerical Results and Discussion

3.1. The Validation of Transmittance and Phase Shift

The effects of control gain (g) and p-values on the phase shift and transmission coefficient were investigated separately for an operating frequency of 6 kHz. The two-dimensional basic cell shown in Figure 2 was constructed in COMSOL (plane strain module) to calculate the transmission coefficient and phase shift of the adaptive unit cell, and the structural parameters are shown in Table 1. The substrate of the adaptive unit cell is made of aluminum materials whose surfaces are covered with piezoelectric materials. The density, elastic modulus, and Poisson’s ratio of aluminum materials are ρs = 2700 kg/m3, Ep = 70 GPa, and υ = 0.33. A piezoelectric material with mass density ρp = 7700 kg/m3, elasticity constant c11 = 70.6 Gpa, piezoelectric constant d31 = −12.6374 C/m2, and dielectric constant ɛ33 = 1.59 × 10−8 F/m was adopted in the present work [50]. A unit stress is applied in the −z direction to the leftmost interface of region I. The grid size is one 140th of the wavelength. The choice of meshing size has been validated.
The effects of g on the wave velocity of sub-cell 1 and sub-cell 2 were investigated according to Equation (12) and are shown in Figure 3a. It can be observed from Figure 3a that the velocity of flexural waves rises as g increases. Therefore, the value of g was chosen to be 100 for sub-cell 1 and 0 for sub-cell 2, respectively. The phase shift and transmittance of adaptive unit cells when p varied from 0 to 1 were also examined with this prerequisite. The results obtained by the TMM were verified by comparison with those derived from the finite element method (FEM). To perform the FEM calculations, a location was selected for the phase and transmittance coefficient calculation in region IV, which was 240 mm away from the line connecting regions Ⅲ and IV. It is worth noting that the absolute value of the phase can be altered if another point at a different location in the transmittance region is selected. However, the difference between the phase shift for the adaptive unit cell is always the same.
Figure 3b shows the transmission coefficient and phase shift versus the value of g when p = 0.5, and Figure 3c shows the phase shift and transmittance with respect to varying p as g was chosen as 100 for sub-cell 1 and 0 for sub-cell 2. In Figure 3b,c, the red solid line indicates the transmission coefficient resolved by the TMM, and the cyan dashed line indicates the numerical solution to the FEM. The black solid line and the dark yellow dashed line indicate the phase shifts obtained from the TMM and FEM calculations, respectively. As can be seen in Figure 3b, there is a very small error in the transmission coefficients calculated using the TMM and FEM; in fact, the maximum difference in values is within 0.05, while the phase gradient variation is achieved by adjusting the phase. We think that this produces a small effect on the tectonic metasurface. As can be seen in Figure 3b,c, the adaptive unit cell fails to cover a phase shift ranging from 0 to 2π when changing the control gain, but it can achieve this range with the change in p. According to the GSL, we used the adaptive unit cells to compose the column elastic metasurface by varying the parameters of the adaptive cells in order to effectively change the phase of each adaptive unit cell. We exploited this feature to construct a passive modulation metasurface by varying the adaptive unit cells of p and further varying g for active modulation based on it.

3.2. Results of Abnormal Transmittance

Based on the generalized Snell’s law, different abnormal refraction strategies can be conducted by adjusting the adaptive unit cells of the metasurface. Additionally, it can be seen from Section 3.1 that the control gain (g) and p-values have an effect on the propagation of the flexural wave. Finite element simulation of flexural wave propagation was conducted using COMSOL Multiphysics software 6.0. Three-dimensional elements are adopted to discretize the incident field, the metasurface, and the transmitted field, whose size accounts for one 180th of the wavelength. The choice of meshing size has been validated. In order to achieve the abnormal refraction of flexural waves, 40 adaptive unit cells were utilized to form the metasurface, with a slit width w2 of 1.98 mm between adaptive unit cells. Perfectly matched layers (PMLs) were arranged around the boundaries of the whole plate to eliminate the influence of reflected waves from boundaries (PML acts as a near-ideal wave absorber). A z-direction displacement −1 × exp(−((y − 0)/cos0)2/602) (mm) was applied at the boundary between the left PML and the incident plate to excite the flexural-wave Gaussian beam, and the wave in the plane surrounded by the black dashed line on the left side in Figure 4 represents the Gaussian beam. The wavelengths of flexural waves corresponding to f = 6 kHz and f = 7 kHz are 69.94 mm and 64.75 mm, respectively, within the 3 mm thick aluminum plate, satisfying the sub-wavelength requirement.
Firstly, the abnormal transmission of passively modulated metasurfaces was realized by changing p. Two different abnormal refraction angles, θt =30° and 38.7°, with operating frequency f = 6 kHz were designed as the targets. The g values for sub-cells 1 and 2 were assumed as 100 and 0, respectively, and the designed values of p were set to match the requirement of phase shift ranging from 0 to 2π. The p-values of the Rith adaptive cell units (Ri = 1 – 40) for different refraction angles in the case of f = 6 kHz are listed in Table A1 (θt = 30°, 38.7°). Full-wave simulations were performed, as shown in Figure 4a,b, from which it can be seen that the vertically incident Gaussian beam can be deflected on demand via the metasurface, realizing the abnormal transmission of the flexural waves. Obviously, the propagation of the transmitted wave matches very well with the theoretical design route.
Then, regarding the active modulation of flexural waves, the geometrical parameters of the unit cells were fixed, and abnormal transmission was achieved by varying the g from the negative proportional feedback control strategy applied in the sub-cells. The target angles for transmitted waves are assumed to be 20° and 38.7° when the incident frequency is 6 kHz, and at 7 kHz, these angles are assumed to be 20°, 30°, and 38.7°, respectively. Additionally, the used values of g for the active modulation are provided in Table A2, and the active tunability performance is shown in Figure 5. The simulated angles and target angles basically coincide with each other, which proves the effectiveness and reliability of the active metasurface based on the piezoelectrical effects and negative proportional feedback control technique. Meanwhile, the angles we designed are coincident at different operating frequencies, which suggests the reproducibility of the study. Moreover, it can be observed in Figure 4 and Figure 5 that the main energy of the beam is concentrated in the designed deflection route, but there is still a small amount of energy appearing outside the designed angle, which belongs to the higher-order diffraction [51].

3.3. Realization of Focus Functionality

In this subsection, the focusing functionality of flexural waves via the passively modulated and actively modulated metasurface is detailed separately. For a positively incident flexural wave, the phase distribution along the x-axis is calculated using Equation (17) within different operating frequencies and focusing point locations. The phase distribution is symmetrically distributed about the y-axis, and thus only 20 adaptive unit cells at one side of the y-axis need to be calculated.
In the FEM simulation, a unit displacement along the -z direction is applied to the left side of the metasurface for the excitation of the flexural wave. We designed passively modulated metasurfaces, which are expected to realize two focal points, (100 mm, 0) and (120 mm, 0), with different focal lengths at a functional frequency of 6 kHz by varying the p and their design values, as shown in Table A1. Figure 6a,b show the numerically calculated out-of-plane displacement plots (z-direction) with the designed focal points (100 mm, 0) and (120 mm, 0) marked with black dashed crosses. It can be clearly seen that the flexural wave is incident vertically on the left side, and the transmitted wave on the right side produces a focusing phenomenon, with a clear focusing phenomenon at the position of the design focal point. To further demonstrate the flexural wave focusing effect obtained with our designed metasurface, the normalized displacement amplitude on the straight line x = x0 where the focal point is located was calculated. Figure 7a,b show the normalized displacement amplitude on the line of x = x0 where the design focus is located, in which the blue solid line is the displacement amplitude at x = x0 and the horizontal blue dashed line denotes the displacement amplitude of the incident wave. From Figure 7, it can be seen that the vibration amplitude at the focal point of the design target is several times larger than the displacement amplitude of the incident wave, and the highest point of the amplitude is at y = 0 mm.
Apart from passive adjustment, an actively metasurface was also designed to achieve the focusing of flexural waves. In this case, different focuses at different operating frequencies were required without changing the geometry of the metasurface. Table A3 presents the design values of g used in the sub-cells for achieving active wave focusing. From Figure 8a,b, it can be seen that the displacements of transmitted wave fields are concentrated in the positions (80 mm, 0) and (120 mm, 0), which is in accordance with the design focal points at 6 kHz. Similarly, at another frequency of 7 kHz, three sets of focal points, (80 mm, 0), (100 mm, 0), and (120 mm, 0), were selected, and the verified performance is shown in Figure 8c–e. Figure 9a–e correspond sequentially to the normalized displacement amplitude on the line x = x0, where the focus is located in Figure 8a–e. As shown in Figure 9, the displacement amplitude at the focusing point far exceeds that of the surrounding region.

3.4. Transforming Guided Wave

In addition to the functions described above, the metasurface can also transform circular waves into plane waves. The transformation of the guided wave was investigated using a focused metasurface model (80 mm, 0) with a displacement in the -z direction applied at coordinates (x,y) = (−80 mm, 0) to excite the circular wave. Figure 10 shows the out-of-plane amplitude plot (z-direction) of a circular wave converted to a plane wave at operating frequencies of (a) f = 6 kHz and (b) 7 kHz. From Figure 10, it can be seen that the circular wave formed by the point source excitation as the incident applied on the left side can be successfully transformed into a plane wave when propagating through the transmitted field, completing the waveform transformation as designed. Numerical simulation results demonstrate that the metasurface not only achieves flexural wave focusing but also transforms circular waves into plane waves.

4. Conclusions

In the present study, a tunable metasurface consisting of adaptive unit cells attached with piezoelectrical actuators and sensors was achieved by applying the negative proportional feedback control strategy. The flexural motion of the metasurface was derived based on the equivalent parameters of the adaptive unit cells. Employing the TMM, the analytical model was established to predict the transmittance and the phase shift of flexural waves propagating through the metasurface, and the validity was verified by comparing the results to those from FEM simulations. Based on the generalized Snell’s law, a metasurface with 40 adaptive unit cells was designed. The passively modulated metasurface was constructed by changing the proportion of the adaptive unit cell occupied by sub-cells, and abnormal modulation flexural waves were achieved at 6 kHz. With the fixed geometrical parameters of the passively modulated metasurface, different phase gradients were generated by adjusting the negative proportional feedback system control gain. Full-wave simulations were performed in finite elements, and the metasurface achieved functions such as the abnormal transmission of multiple angles and the planar lensing of multiple positions at different operating frequencies. A comparison of the theoretical design and finite element simulations shows that this study is reproducible and reliable, which proves the feasibility of the design structure and method. The introduction of a negative proportional positive feedback control strategy provides an active method for adjusting the phase gradient to modulate flexural waves.

Author Contributions

Conceptualization, B.L. and J.L.; methodology, B.L. and J.L.; software, B.L. and J.L.; validation, B.L., J.L. and W.L.; formal analysis, B.L. and J.L.; investigation, J.L., W.L. and Q.M.; resources, B.L. and J.L.; data curation, B.L. and J.L.; writing—original draft preparation, B.L. and J.L.; writing—review and editing, B.L. and J.L.; visualization, B.L., J.L., W.L. and Q.M.; supervision, J.L.; project administration, J.L.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 52365009, the Natural Science Foundation of Hainan Province, grant number 122MS004, and the Fundamental Research Funds for the Central Universities, grant number 3072022JC2601.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author due to privacy reasons.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The following three tables show the p and g with values for the adaptive unit cells at different frequencies. Among them, Table A1 reflects the passively modulated metasurfaces to regulate the flexural wave by changing the geometrical structure, and Table A2 and Table A3 reflect the use of control gain to regulate the flexural wave with fixed metasurface geometrical structure.
Table A1. The p of the Rith adaptive unit cell for abnormal refraction and focusing in passive modulation.
Table A1. The p of the Rith adaptive unit cell for abnormal refraction and focusing in passive modulation.
RiThe p of the Adaptive Unit Cell at 6 kHz
Abnormal TransmittanceFocusing
θt = 30°θt = 38.7°x0 = 100 mmx0 = 120 mm
10.0470.0470.0470.047
20.1220.1530.0550.055
30.2130.2330.0690.064
40.2470.2650.1020.091
50.270.2920.1690.141
60.2920.330.2160.197
70.320.4230.2420.234
80.3750.5030.2630.25
90.4670.5360.2850.27
100.5110.5590.3130.291
110.5360.5880.3660.319
120.5550.6380.4730.375
130.5750.7590.5220.477
140.6020.8160.550.52
150.6580.8610.5780.548
160.7590.9560.6270.573
170.8060.9940.7640.613
180.8410.9940.8280.736
190.8830.9940.8950.813
200.9750.9940.9920.866
210.0470.047
    
220.1220.047
    
230.2130.047
    
240.2470.047
    
250.270.153
    
260.2920.233
    
270.320.265
    
280.3750.292
    
290.4670.33
    
300.5110.423
    
310.5360.503
    
320.5550.536
    
330.5750.559
    
340.6020.588
    
350.6580.638
    
360.7590.759
    
370.8060.816
    
380.8410.861
    
390.8830.956
    
400.9750.994
    
Table A2. The g of the Rith adaptive unit cell with abnormal refraction during active modulation.
Table A2. The g of the Rith adaptive unit cell with abnormal refraction during active modulation.
RiThe g of the Adaptive Unit Cell
6 kHz7 kHz
θt = 20°θt = 38.7°θt = 20°θt = 30°θt = 38.7°
Sub-Cell1Sub-Cell2Sub-Cell1Sub-Cell2Sub-Cell1Sub-Cell2Sub-Cell1Sub-Cell2Sub-Cell1Sub-Cell2
110016.80502.9991007.7491008.7311009.572
210014.50705.81005.3421003.811008.58
310010.64502.0691001.78610013.61006.251
410011.58500.7051001.8671004.5151008.263
510013.75802.3631002.751006.3010.60
610016.4604.4781003.8651008.53316.840
710018.87706.7721004.60210010.54653.0010
810018.22507.9861003.09110010.03854.470
910012.31806.82609.6361005.22322.5180
1010012.59409.37010.8111005.48432.9280
1110016.023014.891014.4721008.34460.3580
1210021.612023.931020.04710013.0371000.456
1310028.6291008.643027.22710019.0061004.62
1410036.43210013.586036.00910025.7561001.162
1510037.54810015.148042.48310026.90110010.512
1610021.8389.43100040.16810013.7751001.579
1710023.97427.056100059.86510015.9261003.479
1810036.20551.7971001.53310010027.18810012.509
1910055.87282.2271008.47310010046.71610030.409
2010000012.36910010001000
2110016.805001007.7491008.7311004.123
2210014.50700.8841005.3421003.811002.966
2310010.64500.0351001.78610013.61000.695
2410011.58501.8571001.8671004.5151001.945
2510013.75804.7431002.751006.3011004.264
2610016.4608.5021003.8651008.5331007.227
2710018.877012.8941004.60210010.54610010.274
2810018.225016.3951003.09110010.03810011.102
2910012.318017.32709.6361005.2231007.587
3010012.594024.902010.8111005.4841009.73
3110016.023039.7014.4721008.34410015.647
3210021.612065.535020.04710013.03710025.399
3310028.62910027.918027.22710019.00610039.729
3410036.43210044.143036.00910025.75610060.514
3510037.54810060.131042.48310026.90110082.342
3610021.83810064.353040.16810013.77521.271100
3710023.97463.580059.86510015.92610059.331
3810036.2054.68801.53310010027.18810035.139
3910055.87224.71408.47310010046.7161008.848
4010001.95012.369100100080.7310
Table A3. The g of the Rith adaptive unit cell with focusing during active modulation.
Table A3. The g of the Rith adaptive unit cell with focusing during active modulation.
RiThe g of the Rith Adaptive Unit Cell
6 kHz7 kHz
x0 = 80 mmx0 = 120 mmx0 = 80 mmx0 = 100 mmx0 = 120 mm
Sub-Cell1Sub-Cell2Sub-Cell1Sub-Cell2Sub-Cell1Sub-Cell2Sub-Cell1Sub-Cell2Sub-Cell1Sub-Cell1
101.43503.53301.7301.9303.534
201.49703.52701.81601.95803.529
301.67403.50801.96202.07703.509
401.48502.98401.74601.73602.985
500.19201.25400.52500.16701.25
600.0600.55900.31100.06500.56
701.03800.9701.12700.66200.971
802.75203.00205.62302.32801.946
901.59903.28304.66304.61203.285
1004.5304.81903.17607.41904.82
1106.83505.32204.99709.41205.324
1205.71403.503.51507.63702.414
13010.6982.94007.347012.012.9470
14021.95718.2750016.2022.06818.2850
15042.82239.7940032.079039.88139.8080
16083.01751.7910060.175071.12751.8080
1710.03310016.46401.45710083.942016.4740
1834.21410019.949021.69210089.628019.9590
1966.64610022.789048.81710093.326022.80
2098.76310017.278076.56510074.579017.2880

References

  1. Yi, K.; Collet, M.; Ichchou, M.; Li, L. Flexural waves focusing through shunted piezoelectric patches. Smart Mater. Struct. 2016, 25, 075007. [Google Scholar] [CrossRef]
  2. Cao, L.; Yang, Z.; Xu, Y.; Fan, S.-W.; Zhu, Y.; Chen, Z.; Li, Y.; Assouar, B. Flexural wave absorption by lossy gradient elastic metasurface. J. Mech. Phys. Solids 2020, 143, 104052. [Google Scholar] [CrossRef]
  3. Shen, Y.; Xu, Y.; Liu, F.; Yang, Z. Metasurface-guided flexural waves and their manipulations. Int. J. Mech. Sci. 2023, 257, 108538. [Google Scholar] [CrossRef]
  4. Chen, A.-L.; Wang, Y.-S.; Wang, Y.-F.; Zhou, H.-T.; Yuan, S.-M. Design of Acoustic/Elastic Phase Gradient Metasurfaces: Principles, Functional Elements, Tunability, and Coding. Appl. Mech. Rev. 2022, 74, 020801. [Google Scholar] [CrossRef]
  5. Wang, W.; Iglesias, J.; Jin, Y.; Djafari-Rouhani, B.; Khelif, A. Experimental realization of a pillared metasurface for flexural wave focusing. APL Mater. 2021, 9, 051125. [Google Scholar] [CrossRef]
  6. Feng, L.; Liu, X.-P.; Lu, M.-H.; Chen, Y.-B.; Chen, Y.-F.; Mao, Y.-W.; Zi, J.; Zhu, Y.-Y.; Zhu, S.-N.; Ming, N.-B. Acoustic Backward-Wave Negative Refractions in the Second Band of a Sonic Crystal. Phys. Rev. Lett. 2006, 96, 014301. [Google Scholar] [CrossRef]
  7. Lu, Y.; Srivastava, A. Level Repulsion and Band Sorting in Phononic Crystals. J. Mech. Phys. Solids 2017, 111, 100–112. [Google Scholar] [CrossRef]
  8. Liu, Z.; Zhang, X.; Mao, Y.; Zhu, Y.Y.; Yang, Z.; Chan, C.T.; Sheng, P. Locally Resonant Sonic Materials. Science 2000, 289, 1734–1736. [Google Scholar] [CrossRef] [PubMed]
  9. Oh, J.H.; Kwon, Y.E.; Lee, H.J.; Kim, Y.Y. Elastic metamaterials for independent realization of negativity in density and stiffness. Sci. Rep. 2016, 6, 23630. [Google Scholar] [CrossRef] [PubMed]
  10. Reynolds, M.; Daley, S. An active viscoelastic metamaterial for isolation applications. Smart Mater. Struct. 2014, 23, 045030. [Google Scholar] [CrossRef]
  11. Sun, S.; He, Q.; Xiao, S.; Xu, Q.; Li, X.; Zhou, L. Gradient-index meta-surfaces as a bridge linking propagating waves and surface waves. Nat. Mater. 2012, 11, 426–431. [Google Scholar] [CrossRef]
  12. Yu, N.; Capasso, F. Flat optics with designer metasurfaces. Nat. Mater. 2014, 13, 139–150. [Google Scholar] [CrossRef]
  13. Zhang, J.; Kosugi, Y.; Otomo, A.; Nakano, Y.; Tanemura, T. Active metasurface modulator with electro-optic polymer using bimodal plasmonic resonance. Opt. Express 2017, 25, 30304–30311. [Google Scholar] [CrossRef]
  14. Qi, S.; Li, Y.; Assouar, B. Acoustic Focusing and Energy Confinement Based on Multilateral Metasurfaces. Phys. Rev. Appl. 2017, 7, 054006. [Google Scholar] [CrossRef]
  15. Xie, B.; Tang, K.; Cheng, H.; Liu, Z.; Chen, S.; Tian, J. Coding Acoustic Metasurfaces. Adv. Mater. 2017, 29, 1603507. [Google Scholar] [CrossRef] [PubMed]
  16. Liang, Z.; Li, J. Extreme Acoustic Metamaterial by Coiling Up Space. Phys. Rev. Lett. 2012, 108, 114301. [Google Scholar] [CrossRef] [PubMed]
  17. Fan, S.-W.; Zhao, S.-D.; Cao, L.; Zhu, Y.; Chen, A.-L.; Wang, Y.-F.; Donda, K.; Wang, Y.-S.; Assouar, B. Reconfigurable curved metasurface for acoustic cloaking and illusion. Phys. Rev. B 2020, 101, 024104. [Google Scholar] [CrossRef]
  18. Li, Y.; Liang, B.; Gu, Z.-M.; Zou, X.-Y.; Cheng, J.-C. Reflected wavefront manipulation based on ultrathin planar acoustic metasurfaces. Sci. Rep. 2013, 3, 2546. [Google Scholar] [CrossRef] [PubMed]
  19. Schwan, L.; Umnova, O.; Boutin, C. Sound absorption and reflection from a resonant metasurface: Homogenisation model with experimental validation. Wave Motion 2017, 72, 154–172. [Google Scholar] [CrossRef]
  20. Shen, Y.; Xu, Y.; Liu, F.; Wang, F.; Yang, Z. 3D-printed meta-slab for focusing flexural waves in broadband. Extrem. Mech. Lett. 2021, 48, 101410. [Google Scholar] [CrossRef]
  21. Zhu, H.; Semperlotti, F. Anomalous Refraction of Acoustic Guided Waves in Solids with Geometrically Tapered Metasurfaces. Phys. Rev. Lett. 2016, 117, 034302. [Google Scholar] [CrossRef]
  22. Zeng, L.; Zhang, J.; Liu, Y.; Zhao, Y.; Hu, N. Asymmetric transmission of elastic shear vertical waves in solids. Ultrasonics 2019, 96, 34–39. [Google Scholar] [CrossRef] [PubMed]
  23. Hu, Y.; Li, M.; Liu, H.; Li, B. Broadband manipulation of flexural waves based on phase-modulated elastic metasurfaces. Eng. Struct. 2023, 275, 115209. [Google Scholar] [CrossRef]
  24. Kim, M.S.; Lee, W.R.; Kim, Y.Y.; Oh, J.H. Transmodal elastic metasurface for broad angle total mode conversion. Appl. Phys. Lett. 2018, 112, 241905. [Google Scholar] [CrossRef]
  25. Tian, Z.; Yu, L. Elastic Phased Diffraction Gratings for Manipulation of Ultrasonic Guided Waves in Solids. Phys. Rev. Appl. 2019, 11, 024052. [Google Scholar] [CrossRef]
  26. Lee, H.; Lee, J.K.; Seung, H.M.; Kim, Y.Y. Mass-stiffness substructuring of an elastic metasurface for full transmission beam steering. J. Mech. Phys. Solids 2018, 112, 577–593. [Google Scholar] [CrossRef]
  27. Xu, Y.; Zhou, X.; Wang, W.; Wang, L.; Peng, F.; Li, B. On natural frequencies of non-uniform beams modulated by finite periodic cells. Phys. Lett. A 2016, 380, 3278–3283. [Google Scholar] [CrossRef]
  28. Cao, L.; Yang, Z.; Xu, Y. Steering elastic SH waves in an anomalous way by metasurface. J. Sound Vib. 2018, 418, 1–14. [Google Scholar] [CrossRef]
  29. Xu, Y.; Cao, L.; Peng, P.; Zhou, X.; Assouar, B.; Yang, Z. Beam splitting of flexural waves with a coding meta-slab. Appl. Phys. Express 2019, 12, 097002. [Google Scholar] [CrossRef]
  30. Zhu, H.; Walsh, T.F.; Semperlotti, F. Total-internal-reflection elastic metasurfaces: Design and application to structural vibration isolation. Appl. Phys. Lett. 2018, 113, 221903. [Google Scholar] [CrossRef]
  31. Liu, Y.; Liang, Z.; Liu, F.; Diba, O.; Lamb, A.; Li, J. Source Illusion Devices for Flexural Lamb Waves Using Elastic Metasurfaces. Phys. Rev. Lett. 2017, 119, 034301. [Google Scholar] [CrossRef]
  32. Cao, L.; Yang, Z.; Xu, Y.; Fan, S.-W.; Zhu, Y.; Chen, Z.; Vincent, B.; Assouar, B. Disordered Elastic Metasurfaces. Phys. Rev. Appl. 2020, 13, 014054. [Google Scholar] [CrossRef]
  33. Su, X.; Lu, Z.; Norris, A.N. Elastic metasurfaces for splitting SV- and P-waves in elastic solids. J. Appl. Phys. 2018, 123, 091701. [Google Scholar] [CrossRef]
  34. Lin, Z.; Xu, W.; Xuan, C.; Qi, W.; Wang, W. Modular elastic metasurfaces with mass oscillators for transmitted flexural wave manipulation. J. Phys. D Appl. Phys. 2021, 54, 255303. [Google Scholar] [CrossRef]
  35. Xu, W.; Zhang, M.; Lin, Z.; Liu, C.; Qi, W.; Wang, W. Anomalous refraction manipulation of Lamb waves using single-groove metasurfaces. Phys. Scr. 2019, 94, 105807. [Google Scholar] [CrossRef]
  36. Zhou, W.; Wang, S.; Wu, Q.; Xu, X.; Huang, X.; Huang, G.; Liu, Y.; Fan, Z. An inverse design paradigm of multi-functional elastic metasurface via data-driven machine learning. Mater. Des. 2023, 226, 111560. [Google Scholar] [CrossRef]
  37. Yuan, S.-M.; Chen, A.-L.; Du, X.-Y.; Zhang, H.-W.; Assouar, B.; Wang, Y.-S. Reconfigurable flexural waves manipulation by broadband elastic metasurface. Mech. Syst. Signal Process. 2022, 179, 109371. [Google Scholar] [CrossRef]
  38. Li, S.; Xu, J.; Tang, J. Tunable modulation of refracted lamb wave front facilitated by adaptive elastic metasurfaces. Appl. Phys. Lett. 2018, 112, 021903. [Google Scholar] [CrossRef]
  39. Shao, S.; Xia, R.; Li, Z. Tunable piezoelectric metasurface for manipulating multi-mode guided waves in plate. Eng. Struct. 2022, 270, 114917. [Google Scholar] [CrossRef]
  40. Ren, T.; Liu, C.; Li, F.; Zhang, C. Active tuning of the vibration band gap characteristics of periodic laminated composite metamaterial beams. J. Intell. Mater. Syst. Struct. 2020, 31, 843–859. [Google Scholar] [CrossRef]
  41. Liu, F.; Shi, P.; Chen, Z.; Shen, Y.; Xu, Y.; Yang, Z. Tunable reflection and broadband absorption of flexural waves by adaptive elastic metasurface with piezoelectric shunting circuits. Smart Mater. Struct. 2023, 32, 055018. [Google Scholar] [CrossRef]
  42. Yaw, Z.; Zhou, W.; Chen, Z.; Lim, C. Stiffness tuning of a functional-switchable active coding elastic metasurface. Int. J. Mech. Sci. 2021, 207, 106654. [Google Scholar] [CrossRef]
  43. Li, B.; Hu, Y.; Chen, J.; Su, G.; Liu, Y.; Zhao, M.; Li, Z. Efficient Asymmetric Transmission of Elastic Waves in Thin Plates with Lossless Metasurfaces. Phys. Rev. Appl. 2020, 14, 054029. [Google Scholar] [CrossRef]
  44. Xu, Y.; Cao, L.; Yang, Z. Deflecting incident flexural waves by nonresonant single-phase meta-slab with subunits of graded thicknesses. J. Sound Vib. 2019, 454, 51–62. [Google Scholar] [CrossRef]
  45. Qiu, H.; Chen, M.; Huan, Q.; Li, F. Steering and focusing of fundamental shear horizontal guided waves in plates by using multiple-strip metasurfaces. Europhys. Lett. 2019, 127, 46004. [Google Scholar] [CrossRef]
  46. Lin, Z.; Tol, S. Electroelastic metasurface with resonant piezoelectric shunts for tunable wavefront control. J. Phys. D Appl. Phys. 2023, 56, 164001. [Google Scholar] [CrossRef]
  47. Chen, Y.; Li, X.; Nassar, H.; Hu, G.-K.; Huang, G. A programmable metasurface for real time control of broadband elastic rays. Smart Mater. Struct. 2018, 27, 115011. [Google Scholar] [CrossRef]
  48. Wang, G.; Chen, S.; Wen, J. Low-frequency locally resonant band gaps induced by arrays of resonant shunts with Antoniou’s circuit: Experimental investigation on beams. Smart Mater. Struct. 2010, 20, 015026. [Google Scholar] [CrossRef]
  49. Zhang, H.; Wen, J.-H.; Chen, S.-B.; Wang, G.; Wen, X.-S. Flexural wave band-gaps in phononic metamaterial beam with hybrid shunting circuits. Chin. Phys. B 2015, 24, 036201. [Google Scholar] [CrossRef]
  50. Li, F.; Zhang, C.; Liu, C. Active tuning of vibration and wave propagation in elastic beams with periodically placed piezoelectric actuator/sensor pairs. J. Sound Vib. 2017, 393, 14–29. [Google Scholar] [CrossRef]
  51. Xu, Y.; Yang, Z.; Cao, L. Deflecting Rayleigh surface acoustic waves by a meta-ridge with a gradient phase shift. J. Phys. D Appl. Phys. 2018, 51, 175106. [Google Scholar] [CrossRef]
Figure 1. (a) A 3D model of the adaptive unit cell with a negative proportional feedback control strategy. (b) A 2D model of the 3D adaptive unit cell in (a).
Figure 1. (a) A 3D model of the adaptive unit cell with a negative proportional feedback control strategy. (b) A 2D model of the 3D adaptive unit cell in (a).
Applsci 14 02717 g001
Figure 2. The basic unit for deriving the transmission coefficient and the phase shift.
Figure 2. The basic unit for deriving the transmission coefficient and the phase shift.
Applsci 14 02717 g002
Figure 3. (a) The phase velocity curve varying with g for f = 6 kHz. (b) The transmission coefficient and phase shift as g increases while p = 0.5. (c) The transmission coefficient and phase shift with respect to varying p at 6 kHz when g is 100 in sub-cell 1 and 0 in sub-cell 2.
Figure 3. (a) The phase velocity curve varying with g for f = 6 kHz. (b) The transmission coefficient and phase shift as g increases while p = 0.5. (c) The transmission coefficient and phase shift with respect to varying p at 6 kHz when g is 100 in sub-cell 1 and 0 in sub-cell 2.
Applsci 14 02717 g003
Figure 4. The out-of-plane displacement diagrams (z-direction) for flexural wave propagation with passive modulation by changing p at f = 6 kHz: the design refraction angles are (a) 30° and (b) 38.7°.
Figure 4. The out-of-plane displacement diagrams (z-direction) for flexural wave propagation with passive modulation by changing p at f = 6 kHz: the design refraction angles are (a) 30° and (b) 38.7°.
Applsci 14 02717 g004
Figure 5. The out-of-plane displacement diagrams (z-direction) for flexural wave propagation with active modulation by changing g at f = 6 and 7 kHz: the design refractive angles are (a) 20°, (b) 38.7°, (c) 20°, (d) 30°, and (e) 38.7°, respectively.
Figure 5. The out-of-plane displacement diagrams (z-direction) for flexural wave propagation with active modulation by changing g at f = 6 and 7 kHz: the design refractive angles are (a) 20°, (b) 38.7°, (c) 20°, (d) 30°, and (e) 38.7°, respectively.
Applsci 14 02717 g005
Figure 6. The out-of-plane displacement diagrams (z-direction) for flexural wave propagation with passive modulation by changing p at f = 6 kHz. The design focal positions were (a) (100 mm, 0) and (b) (120 mm, 0).
Figure 6. The out-of-plane displacement diagrams (z-direction) for flexural wave propagation with passive modulation by changing p at f = 6 kHz. The design focal positions were (a) (100 mm, 0) and (b) (120 mm, 0).
Applsci 14 02717 g006
Figure 7. The normalized amplitude (z-component) of the displacement field along the straight line (a) x0 = 100 mm and (b) x0 = 120 mm when performing focusing (the vertical coordinate is the displacement amplitude).
Figure 7. The normalized amplitude (z-component) of the displacement field along the straight line (a) x0 = 100 mm and (b) x0 = 120 mm when performing focusing (the vertical coordinate is the displacement amplitude).
Applsci 14 02717 g007
Figure 8. The out-of-plane displacement diagrams (z-direction) for flexural wave propagation with active modulation by changing g at f = 6 and 7 kHz. The design focal positions were (a) (80 mm, 0), (b) (120 mm, 0) (c) (80 mm, 0), (d) (100 mm, 0), and (e) (100 mm, 0).
Figure 8. The out-of-plane displacement diagrams (z-direction) for flexural wave propagation with active modulation by changing g at f = 6 and 7 kHz. The design focal positions were (a) (80 mm, 0), (b) (120 mm, 0) (c) (80 mm, 0), (d) (100 mm, 0), and (e) (100 mm, 0).
Applsci 14 02717 g008
Figure 9. The normalized amplitude (z-component) for Figure 8. Values along the straight line are (a) x0 = 100 mm, (b) x0 = 120 mm, (c) x0 = 80 mm, (d) x0 = 100 mm, and (e) x0 = 120 mm.
Figure 9. The normalized amplitude (z-component) for Figure 8. Values along the straight line are (a) x0 = 100 mm, (b) x0 = 120 mm, (c) x0 = 80 mm, (d) x0 = 100 mm, and (e) x0 = 120 mm.
Applsci 14 02717 g009
Figure 10. The out-of-plane displacement diagrams of the transforming guided wave at (a) f = 6 kHz and (b) 7 kHz.
Figure 10. The out-of-plane displacement diagrams of the transforming guided wave at (a) f = 6 kHz and (b) 7 kHz.
Applsci 14 02717 g010
Table 1. Geometry of the adaptive unit cell.
Table 1. Geometry of the adaptive unit cell.
h0h1h2hs1hs2hpLw1
3 mm3 mm1 mm2.6 mm0.6 mm0.2 mm64 mm5 mm
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lin, B.; Li, J.; Lin, W.; Ma, Q. Active Tunable Elastic Metasurface for Abnormal Flexural Wave Transmission. Appl. Sci. 2024, 14, 2717. https://doi.org/10.3390/app14072717

AMA Style

Lin B, Li J, Lin W, Ma Q. Active Tunable Elastic Metasurface for Abnormal Flexural Wave Transmission. Applied Sciences. 2024; 14(7):2717. https://doi.org/10.3390/app14072717

Chicago/Turabian Style

Lin, Bizun, Jingru Li, Wei Lin, and Qingfen Ma. 2024. "Active Tunable Elastic Metasurface for Abnormal Flexural Wave Transmission" Applied Sciences 14, no. 7: 2717. https://doi.org/10.3390/app14072717

APA Style

Lin, B., Li, J., Lin, W., & Ma, Q. (2024). Active Tunable Elastic Metasurface for Abnormal Flexural Wave Transmission. Applied Sciences, 14(7), 2717. https://doi.org/10.3390/app14072717

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop