Research Article  Open Access
Vibration and Sound Transmission Performance of Sandwich Panels with Uniform and Gradient Auxetic Double Arrowhead Honeycomb Cores
Abstract
Auxetic mechanical metamaterials that exhibit a negative Poisson’s ratio (NPR) can be artificially designed to exhibit a unique range of physical and mechanical properties. Novel sandwich structures composed of uniform and gradient auxetic double arrowhead honeycomb (DAH) cores were investigated in terms of their vibration and sound transmission performance stimulated by nonhomogeneous metamaterials with nonperiodic cell geometries. The spectral element method (SEM) was employed to accurately evaluate the natural frequencies and dynamic responses with a limited number of elements at high frequencies. The results indicated that the vibrating mode shapes and deformations of the DAH sandwich models were strongly affected by the patterned gradient metamaterials. In addition, the sound insulation performance of the considered DAH sandwich models was investigated regarding the sound transmission loss (STL) from 1 Hz to 1500 Hz under a normal incident planar wave, and this performance was compared with that for hexagonal honeycomb sandwich panels. A programmable structuralacoustic optimization was implemented to maximize the STL while maintaining a constant weight and high strength. The results showed that the uniform DAH sandwich models with larger NPRs generally exhibited better vibration and acoustic attenuation behaviors and that the optimized gradient increasing NPR models yielded higher STL values than the optimized gradient decreasing NPR models for two specified frequency cases, with improvements of 6.52 dB and 2.52 dB and a higher bending stiffness but a lower overall STL. Thus, sandwich panels consisting of auxetic DAHs can achieve desirable vibroacoustic performance with a higher bending stiffness than conventional hexagonal honeycomb sandwich structures, and the design of gradient DAHs can be extended to obtain optimized vibration and noisecontrol capabilities.
1. Introduction
Metamaterials are artificial structures engineered to achieve unusual properties, and metamaterials that behave mechanically and have a negative Poisson’s ratio (NPR) are called “auxetic” mechanical metamaterials [1–4]. These metamaterials expand rather than contract when stretched under uniaxial loading, differing from most natural materials. There have been important developments in many fields of auxetic metamaterials, including auxetic cellular solids of foams and honeycombs [5–7], microporous polymers [8], 3D reentrants [9], chiral structures [10], and rotating rigid structures [11]. Compared to traditional materials, these auxetic structures have superior mechanical properties such as enhanced shear resistance, indentation resistance, and fracture toughness in various automotive, aerospace, biomedicine, and intelligent system applications [1–4].
In recent decades, cellular solids have been progressively employed [6, 7], and auxetic foams and honeycombs have been used in numerous studies to explore the multiple functionalities of the cellular solid structures. Among the practical engineering applications involving cellular materials, sandwich structures (usually a simplified 2D sandwich beam) with inplane (trusslike or prismatic) honeycomb cores have been explored for vibration and acoustic attenuation purposes, with the sound transmission loss (STL) frequently used as an evaluation metric [12–23]. These materials can be artificially tailored to create highstiffness structures with low volumes and weights in contrast to cellular foams [13]. As a basis for the design of sandwiched core arrangements, previous scholars have reported that sandwich structures with embedded inplane cellular cores transmit less sound than those with outofplane cellular cores [14], and core anisotropy can lead to a higher STL [15]. ElRaheb and Wagner [16] previously investigated the structuralacoustic performance of sandwich panels with periodic square trusslike cores. Ruzzene [17] analyzed the vibration and sound radiation performance of sandwich panels with inplane hexagonal and reentrant honeycomb cores, and Griese et al. [18] studied the effects of prismatic conventional and auxetic hexagonal honeycomb cell geometries on the sound transmission properties of sandwich panels. These studies indicated that sandwich panels with auxetic reentrant cellular cores exhibit better sound insulation behavior than the conventional honeycomb sandwich structures with positive Poisson’s ratios. Spadoni and Ruzzene [19] designed an auxetic chiral trusscore sandwich beam and explored its potential benefit in terms of vibration isolation and sound transmission reduction. To obtain an optimized STL, Denli and Sun [20] performed structuralacoustic optimizations of prismatic random cellular cores for minimum noise radiation in different frequency regions but without considering their manufacturability. Franco et al. [21] optimized sandwich panels with innovative cellular core configurations to minimize the sound radiation responses to various sources of excitation. As an extension of the study mentioned in [18], Galgalikar and Thompson [22] optimized hexagonal honeycomb sandwich panels for maximum STL and discovered that the STL response strongly depends on the numbers of unit cells in the horizontal and vertical directions in addition to the cell internal angles. Li and Yang [23] optimized the STL responses of sandwich panels with hybrid hexagonal honeycomb cores and suggested that hybrid core configurations constitute a broader platform for optimal acoustic designs. Apart from these direct structuralacoustic optimizations to obtain an optimized STL, topology optimization based on the band gaps of periodic metamaterials has been applied for sound and vibration control [13] because the stop bands through which plane waves cannot propagate can be obtained via the corresponding dispersion curves [24, 25].
The aforementioned sandwiched cellular structures are usually characterized by repeating unit cells with a fixed geometry. However, most current metamaterials are generally classified into two types based on their structural arrangements: homogenous (with uniform periodic structures) and inhomogeneous (with nonuniform nonperiodic structures) [26]. Lim [27] theoretically pioneered a nonperiodic functionally gradient beam to attain Poisson curving in 2002. Then, sandwich structures with combined conventional and reentrant hexagonal metamaterials were investigated based on various characteristics, such as the specified objective deformations [28]; the strength and failure properties under bending, compressive and impact loads [29–31]; the dynamic behavior [32, 33]; and the sound transmission performance [23, 34]. Moreover, functionally gradient auxetic metamaterials based on other topologies, such as double arrowhead honeycombs (DAHs), have been developed [35–38] in recent years. These nonuniform metamaterials have attracted considerable interest due to their highorder functionality and broad designability.
Accordingly, a novel sandwich panel with cellular cores of auxetic DAHs is proposed and investigated with respect to its vibration and acoustic characteristics. To the best of our knowledge, such structures have not been previously investigated. In addition, gradient double arrowhead cellular cores are designed and implemented in the direction of the panel thickness based on the features of the nonuniform metamaterials. Programmable structuralacoustic optimization processes are then implemented to design a modifiedgradient DAH for an optimized sandwich panel that radiates less sound at different frequency regions. To precisely predict the structural vibroacoustic performance, the spectral element method (SEM) [17, 19] was employed to accurately compute the natural frequencies and dynamic responses of cellular structures at arbitrary frequencies with a largely reduced number of elements and few degrees of freedom (DOFs) [39, 40].
In the remainder of the paper, Section 2 describes the unit cell geometry and the configurations of the sandwich structures based on DAHs. Section 3 presents SEM formulations to compute the dynamic eigenfrequencies and responses and gives the computational procedures for sound radiation and transmission loss. Next, the vibroacoustic performance of the considered structures is discussed in Section 4. Moreover, a programmable structuralacoustic optimization process for maximum STLs is presented in Section 5. Finally, Section 6 concludes this study.
2. Geometries of the Sandwich Panels with Uniform and Gradient DAH Cores
2.1. Cell Geometry of Auxetic DAH Arrays
As shown in Figure 1, four parameters define the cell geometry of a standard DAH, where θ_{1} and θ_{2} denote the respective internal angles between the two inclined cell ribs and the vertical axis, l is the height of the base triangle in the x direction, and t is the cell wall thickness. L_{x} and L_{y} are the cell dimensions, and each unit cell is associated with a rectangle that measures L_{x} × L_{y}, where L_{x} = 2l. Within a continuous honeycomb array, the cell dimensions can be related to the parameters bywhere the superscript denotes the layer and α = L_{y}/L_{x} is the cell aspect ratio. The cell aspect ratio α determines the entire honeycomb core dimension, and the internal angle sequences Θ_{1} and Θ_{2} define the cell shape. For loading in the y direction, the mechanical indices of Poisson’s ratio and effective elastic modulus are given by Qiao and Chen [35]:where is Young’s modulus of the solid. To design a DAH cellular core, the angle sequence Θ_{2} is given, and Θ_{1} can then be derived through equation (1). The geometry of the honeycomb array is then determined and is scaleindependent with a constant L_{y}/L_{x}. In addition, the feasible constraint of must be imposed to attain the auxetic characteristics, yielding . The cell aspect ratio α utilized in this investigation is 0.2, so that < 68.2°. Here, the thicknesses of the cell wall and face sheets must be further considered to avoid excessive element distortion and overlap.
(a)
(b)
2.2. Geometric Description of Sandwich Panels with Uniform and Gradient DAH Cores
The considered sandwich panels with DAH cores are depicted in Figure 2. A DAH core is sandwiched between two constraining face sheets. Each model has an identical total length of 2 m, and each cellular core comprises 40 cells in the length direction and 5 cells across the thickness of the core. L_{x} is calculated to be 50 mm and L_{y} is 10 mm with α = 0.2. However, the total thickness of each model differs because the unit cell in the top layer has a varying cell angle . The constraint of > arctan(1/5α) is then set to limit the influence resulting from the distinction of the panel thicknesses. Therefore, varies from 45 to 65° without element overlap, and the total thickness ranges from 51.7 to 65 mm.
To maintain a constant weight, the thickness of the two face sheets is constrained by a constant value of t_{s} = 2 mm. In addition, the thickness of the cellular cores is varied to maintain the weight equal to that of a baseline cellular sandwich panel with dimensions of 2 m × 50 mm and a relative density of 0.1. To this point in the study, both the overall length and weight have been held constant to compare the various properties of the sandwich models. The base material for the sandwich panel is aluminum with a Young’s modulus E_{s} = 71.9 GPa, Poisson’s ratio ν = 0.33, and mass density ρ_{s} = 2700 kg/m^{3}. The total weights of the two face sheets and core are 21.6 kg/m and 27.0 kg/m, respectively.
The first three models shown in Figure 2(a) are defined as uniform DAH sandwich panels with identical unit cells. The mechanical indices are plotted in Figure 3(a) based on equations (2) and (3) as θ_{2} increases from 45 to 65°. Both of these values monotonously decrease as θ_{2} increases. Here, we stipulate that a large NPR signifies a Poisson’s ratio with a large absolute value. Thus, the DAHs with larger NPRs yield large effective elastic modulus values. In contrast, as illustrated in Figures 2(b) and 2(c), the last two structures exhibit a varying pattern of DAHs across the panel thickness and are the socalled gradient DAH sandwich panels. As plotted in Figures 3(b) and 3(c), functionally increasing and decreasing sequences of Θ_{2} can generate sandwich models with DAHs of a decreasing NPR (DNPR) and an increasing NPR (INPR), respectively, particularly in this scenario. Regarding the gradient sandwich panels, the DAHs with large NPRs mainly yield large values, as do the uniform sandwich panels.
(a)
(b)
(c)
3. Theories of Vibration and Acoustic Analysis for DAH Cores
3.1. SEM for the Vibration of Cellular Structures
The considered sandwich panel is assumed to be infinite along the zaxis as a simplified 2D planar strain problem, following the method of modeling prismatic cellular cores in references [17–20, 22, 23, 40]. Therefore, the inplane cellular structures are composed of rigidly connected beam elements treated as simplified 2D frames, and the long and thin frames are modeled as Bernoulli–Euler beams with an offplane unit depth [17, 23, 40]. As shown in Figure 4, these elements are oriented in the xOy plane, and the element displacement is described in the local coordinate system (ξ, ψ) that rotates with respect to the global coordinate system (x, y) at an angle . The DOFs of each element are u and , which represent the longitudinal and transverse displacements, respectively.
The partial differential equations of motion that describe the longitudinal and transverse vibrations of a uniform Bernoulli–Euler beam are represented in [17, 18, 23, 39, 40]:where E and ρ are Young’s modulus and mass density of the solid materials and A and I are the sectional area and moment of inertia of the beam, respectively. Note that E = E_{s}/(1 ‒ ν^{2}) here under the planar strain assumption. Equation (4) can be simplified to an ordinary differential equation of harmonic motion at a circular frequency ω as
The homogeneous solution to equation (5) can be obtained in the following form:
In this case, a_{i} (i = 1, …, 6) are integration constants andwhere k_{L} and k_{F} are the wavenumbers of the longitudinal and transverse elastic waves, respectively. The spectral nodal displacement vector d comprises the longitudinal displacements, transverse displacements, and slopes of the frame elements and can be represented bywhere i and j denote the initial and final nodes of the beam element, as shown in Figure 4. Based on the displacements at the element boundaries (ξ = 0 and ξ = L^{(e)}), equation (6) expresses the integration constant a_{i} with respect to d, which finally yieldswhere U(ξ; ω) is the continuous displacement field and N(ξ;ω) is the matrix of the transcendental dynamic shape functions. Following the sign convention shown in Figure 5, the spectral components of the longitudinal tensile forces N, bending moments M, and transverse shear forces Q can be connected with U(x) and W(x) by
The spectral nodal longitudinal tensile force, transverse shear force, and bending moment can be related to the corresponding force and moment matrix f bywhere the nodal force vector f_{c}(ω) is associated with the concentrated dynamic forces that are directly applied at the spectral nodes, and f_{d}(ω) is associated with the distributed dynamic forces. A spectral distributed load P(x; ω) must be transferred to each node of the element by
Substituting equation (9) into (10) and implementing the results into the righthand side of equation (11), we obtainwhere is the symmetric matrix of general spectral framed elements. Equation (13) must be transformed from the local coordinate system to the global coordinate system using the coordinate transformation matrix T:where
Finally, the spectral equation for each element is assembled, and the spectral equation of the entire system is given bywhere is the overall spectral element matrix, is the overall nodal displacement vector, and is the overall nodal force vector.
Then, the eigenvalue problem for a spectral element structural system can be reduced from equation (16) toand the infinite eigenfrequencies ω_{i} can be determined using the condition thatat ω = ω_{i}. Here, is a transcendental function of frequency ω, and thus equation (18) is a transcendental eigenvalue problem with exact eigensolutions. The Wittrick–Williams (W–W) algorithm [41–43] is employed here to search for the eigenfrequencies without missing any eigenvalues. In this algorithm, the total number of eigenfrequencies below an arbitrarily chosen frequency can be given bywhere is the total number of eigenfrequencies below , is the sign count of and is computed as the number of negative elements on the diagonal of the upper triangular matrix of that is obtained by Doolittle decomposition. is the total number of eigenfrequencies below under the condition that all displacements are constrained to zero. for simplified planar frame elements can be obtained fromafter dividing the complete structure into N frame members with fixed ends, where and are the numbers of the longitudinal and flexural eigenfrequencies below for the frame member with fixed boundaries, respectively. The iteration procedure implemented to compute the eigenfrequency below resembles the dichotomies algorithm, which was detailed in the study of Li and Yang [23].
For a significantly minimized system scale, the total numbers of nodes and elements are 447 and 881 for each sandwich model, respectively. Since the SEM can provide exact frequencydomain solutions without refined mesh discretization as the frequency increases [39], the mesh connections remain unchanged regardless of the material properties, structure size, or targeted frequencies. Although the spectral matrices in equation (16) need repeated assembly and computation at each frequency step, the SEM still provides an extremely accurate and timesaving alternative for dynamic problems involving cellular sandwich structures.
3.2. Sound Radiation and Transmission Loss Analysis
The considered sandwich model with DAH cores is exposed to a planar incident wave of acoustic pressure on the bottom face sheet as shown in Figure 6. The model is located between two rigid baffles of infinite length. The vibration of the bottom face sheet is transmitted through the cellular core to the upper face sheet; then, the sandwich structure radiates sound in the fluid domain. The ends of the structure, both ends of the face sheets, and the cellular core ribs adjacent to the ends, are all constrained in the x and y directions with pinned boundary conditions. The contribution of the airborne transmission to the total radiated sound power is neglected because of the high structural stiffness [40, 44], and thus only the structureborne transmission is considered by following most of the relevant studies referred in Section 1.
The acoustic incident wave is harmonic with a unit amplitude and an incidence angle α with respect to the yaxis. The incident sound power over the bottom face sheet is given bywhere and represent the fluid density and sound speed, respectively, and L_{B} represents the bottom face sheet. The spectral nodal force can be obtained by equation (12). The transmitted sound pressure at an observation point can be computed by Rayleigh’s first integral for onedimensional baffled radiators in two dimensions [45]:where j is an imaginary unit, k = ω/c_{a} is the acoustic wavenumber, r is the position of the fluid particles on the upper face sheet L_{U} with a normal velocity , is the distance between the fluid particle and observation point, and denotes a Hankel function of the first kind. The transmitted sound power W_{t} can be derived via integration over the upper face sheet [17–19, 21–23] by substituting equation (22) intowhere denotes the complex conjugate. The spectral element shape functions are transcendental in the integrands, and a fifthorder Gaussian quadrature is used to compute the numerical integral of the radiated power with a compromise between computational accuracy and efficiency. The angle of wave incidence (α = 0°) is specified for simplicity in this study. Accordingly, the transmission loss is defined as
The STL characteristics of a panel can be generally divided into four distinct regions: stiffness, resonance, mass, and coincidence regions from low to high frequencies [18, 25]. Moreover, the theoretical STL of an infinitely long and thin panel was estimated by Griese et al. [18]:where m_{p} is the mass per unit surface area of the panel and f is the frequency. The final expression is known as the mass law for sound transmission loss under a normal planar incident wave when α = 0° and m_{p} = 48.6 kg/m for the considered models that are assumed as single panels in this scenario.
4. Vibration and Sound Transmission Performance of the Sandwich Panels with DAH Cores
4.1. Natural Frequencies
The considered models are treated as undamped linear elastic systems because the effects of damping on natural frequency searching are generally small [39]. The sandwich panels are simplysupported by pin joints, as illustrated in Figure 6. The natural frequencies that were computed via the W–W algorithm (SEM) are listed in Table 1 and plotted in Figure 7 for each considered model.

(a)
(b)
For the uniform models 1 to 5, the fundamental frequency that represents the structural bending stiffness decreases as the feature angle θ_{2} of the DAHs increases. This result is associated with the low of the unit cells with a large feature angle, as shown in Figure 3(a), and the deviations in the panel thicknesses also have an effect. However, the highorder natural frequencies of the uniform models with large θ_{2} values are larger than those of the models with small θ_{2} values until the feature angle exceeds 60°. This result implies that the static mechanical indices cannot completely reflect the structural dynamic behavior at highorder natural frequencies. Concerning the gradient models, the natural frequencies of model 7 are always higher than the frequencies of model 6, mainly due to the large differences in the panel thickness.
The vibrating modes of the corresponding natural frequencies for the representative models are depicted in Figure 8. The modes of a global vibrating mode shape correspond to global natural frequencies, and the others correspond to local natural frequencies. The figures in the left column illustrate the third global bending modes when the entire sandwich panel bends in the form of a curved beam in two dimensions. The figures in the right column present the foremost local frequencies when the cellular cores are partly distorted. It can be observed that the foremost local frequencies occur when the initially straight curves become bended as plotted in Figure 7. The varying core configurations of the considered models yield diverse mode shapes, and the gradient DAHs within nonuniform cellular cores further provide flexible vibrating deformation mechanisms.
4.2. Dynamic Responses and Sound Radiation
For dynamic responses, damping can be introduced by a complex modulus E^{∗} = E_{s}(1 + jη) with a structural damping factor η = 0.01; in this case, aluminum damping is considered in addition to the damping effects of the hinged connections of the cell ribs [17]. The acoustic loading and boundary conditions are the same as illustrated in Figure 6, and the dynamic deformation of the considered models is depicted in Figure 9 at 500, 1000, and 1400 Hz.
To visualize the sound radiation distributions after vibrating deformations as shown in Figure 9, the radiated sound pressure levels (SPLs) are also presented in Figure 10 for models 3 and 6 at 500, 1000, and 1400 Hz. The amplitude of the incident pressure is 1 Pa (94 dB), and the reference pressure p_{0} is 20 μPa. As plotted in Figure 10, the fluid observation area ranges from nondimensionalized values of x = –0.5 to x = 0.5 and from nondimensionalized values of y = 0 to y = 1. Additionally, the sandwich panels are located in the middle between x = –0.1 and x = 0.1 at y = 0.
(a)
(b)
4.3. Sound Transmission Loss
In this study, the STL responses of the considered models with varying tailored DAH cores were primarily investigated in the region from 1 to 1500 Hz, which covers the stiffness and resonance regions based on the natural frequencies discussed in Section 4.1. Thus, the acoustic performance at higher frequencies falls beyond the scope of this research. Two uniform hexagonal honeycomb sandwich models based on the configurations in [18, 22] are depicted in Figure 11, for comparison with the DAH sandwich panels. The compared hexagonal honeycomb sandwich models have the same unit cell arrangement of 40 × 5, overall dimensions of 2 m × 50 mm, and total weight of m_{p} = 48.6 kg/m for the baseline panel introduced in Section 2.2.
To verify the effectiveness and accuracy of the vibroacoustic calculation method in this paper, the numeric results compared with those of [18] are shown in Figure 12, which indicates that the sound transmission results via the proposed method agree quantitatively with those of the published materials. Moreover, the STLs of the considered models from 1 to 1500 Hz are shown in Figures 13(a)–13(c) where the mass law trends are also plotted. The STL dips as the low frequencies align with the oddorder bending modes listed in Table 1 because the considered symmetric models strongly vibrate under symmetric excitation when odd bending resonance is encountered. The effective Young’s modulus is selected as a function of the panel stiffness following [18, 22]. Although the DAHs in the considered models yield a smaller than the hexagonal honeycombs as shown in Figures 3 and 11, the modal density of the DAH sandwich models is much larger than the modal density of the sandwich panels with hexagonal cellular cores. This finding implies that sandwich panels with a DAH core yield a larger inplane bending stiffness.
(a)
(b)
The overall average STL (STL_{o}) of the considered models is shown in Figure 13(d) to assess the STL performance as an acoustic evaluation metric. A model that radiates less sound yields a larger STL_{o}. It can be concluded that a uniform sandwich model with a small θ_{2} or large NPR yields the best sound insulation performance, as shown in Figures 13(a) and 13(d). The STL_{o} of model 6 is larger than the STL of model 7, as shown in Figure 13(b), and the STL curve shapes of models 6 and 7 resemble the STL shapes of the uniform models (models 1 and 5) that feature the same θ_{2} for the bottom layer, respectively. As plotted in Figure 13(c), the STL_{o} of model 9 with auxetic hexagonal honeycombs is larger than that of model 8, which is in agreement with the findings in [17, 18, 22]. Moreover, although both models 8 and 9 with conventional hexagonal honeycomb cores transmit less sound than the considered DAH sandwich models, their fundamental and highorder natural frequencies are far lower than those of the sandwich models with DAHs, as illustrated in Figure 13(c). This result indicates that the DAH sandwich panels have a stronger bending stiffness with desirable noise attenuation performance.
5. Programmable StructuralAcoustic Optimization of the Gradient DAHs for Maximum STLs
5.1. Mathematical Formulations for StructuralAcoustic Shape Optimizations
The structuralacoustic optimization problem is defined aswhere and are the lower and upper bounding vectors of Θ_{2}, respectively; STL_{a}(Θ_{2}) is the average STL between the specified frequencies of f_{1} and f_{2}; f_{0}(Θ_{2}) is the fundamental frequency of the optimized model; f_{0} is the minimum allowable natural frequency; and β is the interval of set to avoid element overlap. In this scenario, = 45°, = 65°, and β = 10°. The second and third constraints limit the DAHs with gradient Θ_{2} sequences, where the positive sign represents a decreasing gradient Θ_{2} (a gradient INPR) and the negative sign represents an increasing gradient Θ_{2} (a gradient DNPR). The baseline model was employed with model 3 due to its moderate mechanical and acoustic performance, with values of f_{0} = 83.7 Hz and STL_{o} = 39.67 dB.
The programmable design was achieved via a global optimization method based on the MultiStart algorithm. As a workflow [23], the MultiStart algorithm can automatically distribute the tasks in a problem and simultaneously generate a series of stochastic starting points for multiple processes or processors with parallel computing. Subsequently, the problems run independently from the corresponding start points, and the MultiStart algorithm finally combines the distinct local minima into a global vector to obtain a relatively optimized solution, where the chosen local solver involves sequential quadratic programming (SQP). In this scenario, the procedure was executed with 1 × 10^{2} random start points, and the iteration tolerance for the local solver was 1 × 10^{−3} to maintain computational efficiency and reasonability. The optimization procedures were executed for specified tonal and frequency band cases.
5.2. Optimal Design for Specified Tonal and Frequency Band Cases
Here, a tonal excitation at 1400 Hz, which was close to an STL dip in the baseline model, was chosen for the tonal case. The STLs of the baseline and optimized gradient DAH sandwich models are compared in Figure 14(a), and the iterative curves converging to the optimum solutions are plotted in Figure 14(b). The configurations of the optimized models are shown in Figure 14(c), and the mechanical properties of the gradient DAH cores are given in Figure 14(d). In addition, the key parameters of the optimal design are listed and compared in Table 2, and both of the optimized models satisfy the high strength constraint and yield feasible core geometries. The STL_{a} increments are approximately 9.72 dB for the optimized gradient DNPR model and 16.24 dB for the optimized gradient INPR model. The SPL distributions and dynamic deformations at 1400 Hz are illustrated in Figure 15 for the optimized models. As illustrated in Figures 15(b) and 14(a) for the optimized gradient INPR model, localized resonance occurs and the STL exceeds the mass law curve value at the target frequency. The double arrowhead unit cells adjacent to the bottom face sheet absorb the majority of the vibrating energy; hence, the optimized gradient INPR model effectively transmits the least sound as a result.
(a)
(b)
(c)
(d)

(a)
(b)
Here, the chosen computational frequency interval was 20 Hz for the frequency band from 1000 to 1500 Hz. The optimized information is depicted in Figure 16, and the characteristic parameters are listed in Table 3. The results demonstrate that the optimized solutions strictly satisfy the constraints of the bending stiffness and avoid element overlap. Furthermore, the reduction in average sound radiation power is approximately 2.51 dB for the optimized gradient DNPR model and 5.03 dB for the optimized gradient INPR model compared to the STL_{a} of the baseline model. More optimized solutions could be obtained by increasing the number of start points, decreasing the iteration tolerance, or using larger integral orders.
(a)
(b)
(c)
(d)

5.3. Comparison and Discussion of the Sound Transmission Performance of the Optimized and Considered Models
The optimized gradient INPR models can achieve better insulation of radiated sound power than can the optimized gradient DNPR models, with improvements of 6.52 dB and 2.52 dB for the tonal and frequency band cases in this scenario, respectively. As plotted in Figures 14(d) and 16(d), the double arrowhead unit cells adjacent to the fluid domain feature a large and yield a low in the gradient INPR models, and vice versa. The gradient INPR models enable a programmable DAH core for a high STL_{a} that even exceeds the mass law curve values and yield a high bending stiffness; however, the use of gradient INPRs generally sacrifices overall sound insulation, with a lower STL_{o}. It can be concluded that the DAHs with large NPRs in upper layers improve the potential for obtaining an optimized STL within a specified frequency region, following the phenomenon that a uniform DAH sandwich model with a large NPR yields high sound insulation performance. Referring to the STL curves in Figures 14(a) and 16(a), the DAHs in the bottom layers mainly affect the STL curve shapes, as discussed in Section 4.3, and therefore, the overall noise insulation properties are limited.
6. Conclusions
Sandwich structures composed of uniform and gradient DAH cores across the panel thickness were presented, and the vibration and sound transmission properties of the proposed structures were evaluated by comparing uniform and gradient DNPR and INPR models, as well as uniform hexagonal honeycomb sandwich panels. The following conclusions were drawn from the study:(1)Regarding the vibration properties, the various core configurations of the gradient sandwich DAH panels yield more diverse mode shapes and provide more flexible vibrating deformation mechanisms than do the uniform models.(2)With respect to the sound transmission performance with normal incidence, the uniform DAH sandwich models with large NPRs generally yield better sound insulation behaviors. Moreover, the considered DAH sandwich models insulate less sound but are far stiffer than the conventional hexagonal honeycomb sandwich models.(3)The programmable optimized gradient INPR models yield higher STL values than the optimized DNPR models, with improvements of 6.52 dB and 2.52 dB for the specified frequency cases and a large bending stiffness; however, a lower overall sound insulation was observed. For the gradient sandwich models, the DAHs with large NPRs in the upper layers improve the potential for obtaining an optimized target STL, and the DAHs in the bottom layers affect mainly the STL curve shapes and govern the overall noise insulation properties.(4)Further consideration could be given to the vibroacoustic performance and design of the proposed structures at high frequencies and under other excitation conditions. Furthermore, the design variables can also be extended to consider more than a single feature angle. Moreover, further studies could be conducted to investigate gradient auxetic metamaterials with varying topologies.
Data Availability
The programming data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The support for this work provided by the National Natural Science Foundation of China (51479115), HighTech Ship Research Projects by MIIT ([2014]148 and [2016]548), and Opening Project by the State Key Laboratory of Ocean Engineering (GKZD010071) is gratefully acknowledged.
References
 C. Huang and L. Chen, “Negative Poisson’s ratio in modern functional materials,” Advanced Materials, vol. 28, no. 37, pp. 8079–8096, 2016. View at: Publisher Site  Google Scholar
 K. E. Evans and A. Alderson, “Auxetic materials: functional materials and structures from lateral thinking!,” Advanced Materials, vol. 12, no. 9, pp. 617–628, 2000. View at: Publisher Site  Google Scholar
 Y. Prawoto, “Seeing auxetic materials from the mechanics point of view: a structural review on the negative Poisson’s ratio,” Computational Materials Science, vol. 58, pp. 140–153, 2012. View at: Publisher Site  Google Scholar
 H. M. A. Kolken and A. A. Zadpoor, “Auxetic mechanical metamaterials,” RSC Advances, vol. 7, no. 9, pp. 5111–5129, 2017. View at: Publisher Site  Google Scholar
 R. Lakes, “Foam structures with a negative Poisson’s ratio,” Science, vol. 235, no. 4792, pp. 1038–1040, 1987. View at: Publisher Site  Google Scholar
 L. J. Gibson and M. F. Ashby, Cellular Solids: Structure and Properties, Cambridge University, Cambridge, UK, 1999.
 Q. Zhang, X. Yang, P. Li et al., “Bioinspired engineering of honeycomb structure  using nature to inspire human innovation,” Progress in Materials Science, vol. 74, pp. 332–400, 2015. View at: Publisher Site  Google Scholar
 K. E. Evans, B. D. Caddock, and M. J. Nobes, “Microporous materials with negative Poisson’s ratios. II. Mechanisms and interpretation,” Journal of Physics D: Applied Physics, vol. 22, no. 12, pp. 1883–1887, 1989. View at: Publisher Site  Google Scholar
 S. Babaee, J. Shim, J. C. Weaver, E. R. Chen, N. Patel, and K. Bertoldi, “3D soft metamaterials with negative Poisson's ratio,” Advanced Materials, vol. 25, no. 36, pp. 5044–5049, 2013. View at: Publisher Site  Google Scholar
 R. Lakes, “Deformation mechanisms in negative Poisson’s ratio materials: structural aspects,” Journal of Materials Science, vol. 26, no. 9, pp. 2287–2292, 1991. View at: Publisher Site  Google Scholar
 J. N. Grima and K. E. Evans, “Auxetic behavior from rotating squares,” Journal of Materials Science Letters, vol. 19, no. 17, pp. 1563–1565, 2000. View at: Publisher Site  Google Scholar
 V. D’Alessandro, G. Petrone, F. Franco, and S. De Rosa, “A review of the vibroacoustics of sandwich panels: models and experiments,” Journal of Sandwich Structures and Materials, vol. 15, no. 5, pp. 541–582, 2013. View at: Publisher Site  Google Scholar
 J. Wu, “Topology optimization studies for light weight acoustic panels,” University of Toronto, Toronto, Canada, 2016, M.S. thesis. View at: Google Scholar
 J. A. Moore and R. H. Lyon, “Sound transmission loss characteristics of sandwich panel constructions,” Journal of the Acoustical Society of America, vol. 89, no. 2, pp. 777–791, 1991. View at: Publisher Site  Google Scholar
 P. Thamburaj and J. Q. Sun, “Optimization of anisotropic sandwich beams for higher sound transmission loss,” Journal of Sound and Vibration, vol. 254, no. 1, pp. 23–36, 2002. View at: Publisher Site  Google Scholar
 M. ElRaheb and P. Wagner, “Transmission of sound across a trusslike periodic panel; 2D analysis,” Journal of the Acoustical Society of America, vol. 102, no. 4, pp. 2176–2183, 1997. View at: Publisher Site  Google Scholar
 M. Ruzzene, “Vibration and sound radiation of sandwich beams with honeycomb truss core,” Journal of Sound and Vibration, vol. 277, no. 45, pp. 741–763, 2004. View at: Publisher Site  Google Scholar
 D. Griese, J. D. Summers, and L. Thompson, “The effect of honeycomb core geometry on the sound transmission performance of sandwich panels,” Journal of Vibration and Acoustics, vol. 137, no. 2, Article ID 021011, 2014. View at: Google Scholar
 A. Spadoni and M. Ruzzene, “Structural and acoustic behavior of chiral trusscore beams,” Journal of Vibration and Acoustics, vol. 128, no. 5, pp. 616–626, 2006. View at: Publisher Site  Google Scholar
 H. Denli and J. Q. Sun, “Structuralacoustic optimization of sandwich structures with cellular cores for minimum sound radiation,” Journal of Sound and Vibration, vol. 301, no. 12, pp. 93–105, 2007. View at: Publisher Site  Google Scholar
 F. Franco, K. A. Cunefare, and M. Ruzzene, “Structuralacoustic optimization of sandwich panels,” Journal of Vibration and Acoustics, vol. 129, no. 3, pp. 330–340, 2007. View at: Publisher Site  Google Scholar
 R. Galgalikar and L. L. Thompson, “Design optimization of honeycomb core sandwich panels for maximum sound transmission loss,” Journal of Vibration and Acoustics, vol. 138, no. 5, Article ID 051005, 2016. View at: Publisher Site  Google Scholar
 Q. Li and D. Yang, “Mechanical and acoustic performance of sandwich panels with hybrid cellular cores,” Journal of Vibration and Acoustics, vol. 140, no. 6, Article ID 061016, 2018. View at: Publisher Site  Google Scholar
 A. S. Phani, J. Woodhouse, and N. A. Fleck, “Wave propagation in twodimensional periodic lattices,” Journal of the Acoustical Society of America, vol. 119, no. 4, pp. 1995–2005, 2016. View at: Publisher Site  Google Scholar
 S. Ehsan Moosavimehr and A. Srikantha Phani, “Sound transmission loss characteristics of sandwich panels with a truss lattice core,” Journal of the Acoustical Society of America, vol. 141, no. 4, pp. 2921–2932, 2017. View at: Publisher Site  Google Scholar
 T. J. Cui, M. Q. Qi, X. Wan, J. Zhao, and Q. Cheng, “Coding metamaterials, digital metamaterials and programmable metamaterials,” Light: Science and Applications, vol. 3, no. 10, p. e218, 2014. View at: Publisher Site  Google Scholar
 T. C. Lim, “Functionally graded beam for attaining Poissoncurving,” Journal of Materials Science Letters, vol. 21, no. 24, pp. 1899–1901, 2002. View at: Publisher Site  Google Scholar
 Y. Han and W. Lu, “Evolutionary design of nonuniform cellular structures with optimized Poisson's ratio distribution,” Materials and Design, vol. 141, pp. 384–394, 2018. View at: Publisher Site  Google Scholar
 Y. Hou, Y. H. Tai, C. Lira, F. Scarpa, J. R. Yates, and B. Gu, “The bending and failure of sandwich structures with auxetic gradient cellular cores,” Composites Part A: Applied Science and Manufacturing, vol. 49, pp. 119–131, 2013. View at: Publisher Site  Google Scholar
 Y. Hou, R. Neville, F. Scarpa, C. Remillat, B. Gu, and M. Ruzzene, “Graded conventionalauxetic Kirigami sandwich structures: flatwise compression and edgewise loading,” Composites Part B: Engineering, vol. 59, pp. 33–42, 2014. View at: Publisher Site  Google Scholar
 C. Yang, H. D. Vora, and Y. Chang, “Behavior of auxetic structures under compression and impact forces,” Smart Materials and Structures, vol. 27, no. 2, Article ID 025012, 2018. View at: Publisher Site  Google Scholar
 C. Lira, F. Scarpa, and R. Rajasekaran, “A gradient cellular core for aeroengine fan blades based on auxetic configurations,” Journal of Intelligent Material Systems and Structures, vol. 22, no. 9, pp. 907–917, 2011. View at: Publisher Site  Google Scholar
 L. Boldrin, S. Hummel, F. Scarpa et al., “Dynamic behaviour of auxetic gradient composite hexagonal honeycombs,” Composite Structures, vol. 149, pp. 114–124, 2016. View at: Publisher Site  Google Scholar
 M. S. Mazloomi, M. Ranjbar, L. Boldrin, F. Scarpa, S. Patsias, and N. Ozada, “Vibroacoustics of 2D gradient auxetic hexagonal honeycomb sandwich panels,” Composite Structures, vol. 187, pp. 593–603, 2018. View at: Publisher Site  Google Scholar
 J. X. Qiao and C. Q. Chen, “Impact resistance of uniform and functionally graded auxetic double arrowhead honeycombs,” International Journal of Impact Engineering, vol. 83, pp. 47–58, 2015. View at: Publisher Site  Google Scholar
 T.C. Lim, “A 3D auxetic material based on intersecting double arrowheads,” Physica Status Solidi (b), vol. 253, no. 7, pp. 1252–1260, 2016. View at: Publisher Site  Google Scholar
 N. Novak, M. Vesenjak, and Z. Ren, “Computational simulation and optimization of functionally graded auxetic structures made from inverted tetrapods,” Physica Status Solidi (b), vol. 254, no. 12, p. 1600753, 2017. View at: Publisher Site  Google Scholar
 W. Hou, X. Yang, W. Zhang, and Y. Xia, “Design of energydissipating structure with functionally graded auxetic cellular material,” International Journal of Crashworthiness, vol. 23, no. 4, pp. 366–376, 2017. View at: Publisher Site  Google Scholar
 U. Lee, Spectral Element Method in Structural Dynamics, John Wiley & Sons, New York, USA, 2009.
 H. Yang, H. Li, and H. Zheng, “A structuralacoustic optimization of twodimensional sandwich plates with corrugated cores,” Journal of Vibration and Control, vol. 23, no. 18, pp. 3007–3022, 2017. View at: Publisher Site  Google Scholar
 W. H. Wittrick and F. W. Williams, “A general algorithm for computing natural frequencies of elastic structures,” The Quarterly Journal of Mechanics and Applied Mathematics, vol. 24, no. 3, pp. 263–284, 1971. View at: Publisher Site  Google Scholar
 F. W. Williams, “An algorithm for exact eigenvalue calculations for rotationally periodic structures,” International Journal for Numerical Methods in Engineering, vol. 23, no. 4, pp. 609–622, 1986. View at: Publisher Site  Google Scholar
 S. Yuan, K. Ye, and F. W. Williams, “Second order modefinding method in dynamic stiffness matrix methods,” Journal of Sound and Vibration, vol. 269, no. 3–5, pp. 689–708, 2004. View at: Publisher Site  Google Scholar
 L. Cheng, Y. Y. Li, and J. X. Gao, “Energy transmission in a mechanicallylinked doublewall structure coupled to an acoustic enclosure,” Journal of the Acoustical Society of America, vol. 117, no. 5, pp. 2742–2751, 2005. View at: Publisher Site  Google Scholar
 E. G. Williams, Fourier acoustics: Sound radiation and Nearfield Acoustical Holography, Academic, London, UK, 1999.
Copyright
Copyright © 2019 Qing Li and Deqing Yang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.