Next Article in Journal
Influence of Natural Fractures and Laminae on Fracture Propagation and Failure Mode of Continental Shale
Previous Article in Journal
A Lab-Scale Evaluation of Parameters Influencing the Mechanical Activation of Kaolin Using the Design of Experiments
Previous Article in Special Issue
Bandgap Calculation and Experimental Analysis of Piezoelectric Phononic Crystals Based on Partial Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Cross Tetrachiral Honeycomb Metamaterial with Designable Static and Dynamic Performances

1
State Key Laboratory for Turbulence and Complex Systems, Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China
2
Department of Applied Mechanics and Engineering, School of Aeronautics and Astronautics, Sun Yat-sen University, Shenzhen 518107, China
3
Mechanical and Materials Engineering, University of Nebraska-Lincoln, Lincoln, NE 68588, USA
*
Authors to whom correspondence should be addressed.
Materials 2024, 17(18), 4652; https://doi.org/10.3390/ma17184652
Submission received: 8 August 2024 / Revised: 14 September 2024 / Accepted: 16 September 2024 / Published: 23 September 2024
(This article belongs to the Special Issue Acoustic and Mechanical Metamaterials: Recent Advances)

Abstract

:
A novel cross tetrachiral honeycomb metamaterial is proposed, which not only possesses the negative Poisson’s ratio property, but also has a wide-frequency bandgap. The effective elastic parameters of the cross tetrachiral honeycomb are first theoretically analyzed; then, its designable performances for negative Poisson’s ratio and elastic modulus are studied by varying geometric parameters. The dynamic properties of the cross tetrachiral honeycomb metamaterial are investigated by analyzing the band structure. It is shown that without the addition of external mass to the structure, a designable wide bandgap can be generated to isolate the in-plane waves effectively by selecting the ligament angles and the radius of central cylinder. In addition, an effective approach is proposed for tuning the bandwidth without changing the geometric parameters of the structure. Compared to classical negative Poisson’s ratio metamaterials, the proposed cross tetrachiral honeycomb metamaterial is designable and tunable for achieving a specific static or dynamic performance, and has potential applications in engineering practice.

1. Introduction

Mechanical metamaterials are artificially designed materials with mesoscale or microscale structural arrangements, and have extraordinary properties that are not found in natural materials. Their mechanical properties do not just depend on their material properties, but can be designed by the mesoscale or microscale structures [1,2,3]. Therefore, metamaterials have a variety of structural types obtained through designs to meet broad engineering requirements in defense and military, aerospace, civil construction, vehicle manufacturing, shipbuilding engineering, and other industries. Metamaterials have been used in some practical applications in specific projects such as the skin filling of deformable airfoils [4], automotive comfort optimization [5], clinical diagnostics, and biomedical testing [6].
Based on elaborate structural designs, mechanical metamaterials are adjustable [7] and programmable [8,9,10], and can achieve unique physical and mechanical properties [11,12,13,14], such as negative compressibility [15], negative thermal expansion [16], negative Poisson’s ratio [17], and mechanical vibration isolation [18]. Among these, metamaterials with negative Poisson’s ratio, i.e., auxetic materials [19], have attracted more attention, and many structures with this property have been reported, including trichiral and anti-trichiral, tetrachiral and anti-tetrachiral, and hexachiral honeycomb structures [20,21,22]. Beyond the negative Poisson’s ratio, the dynamic behaviors of these metamaterials have also been explored. They have shown good performances in impact resistance [23], fracture tolerance [24], energy absorption, and vibration isolation [25].
Yajun Xin et al. [26] chose a star-shaped negative Poisson’s ratio structure, and proposed a novel re-entrant method to adjust low frequency bandgaps by adding rubber-coated masses. As a typical type of auxetic metamaterial, the chiral honeycomb metamaterial has been developed to obtain both excellent mechanical properties and vibration isolation. Luyun Chen et al. [27] have analyzed the geometric parameters and material properties of hexachiral metamaterials to obtain a wide and complete bandgap. However, most of the studies in this area focus on the generation of bandgaps by additional mass. Xiaoning Liu et al. [28] added a softly coated heavy cylinder into the two-dimensional periodic chiral lattice to achieve low-frequency bandgaps. Pei Sun et al. [29] proposed an internal resonance unit with cruciate ligament features to replace the metamaterial inclusion for extreme low-frequency bandgaps. Similarly, based on a chiral lattice, Rui Zhu et al. [18] first designed an elastic metamaterial beam with distributed multiple inner resonators for achieving wide-band vibration suppression without sacrificing the load-bearing capacity of the beam. Furthermore, four-ligament chiral honeycomb structures have attracted more attention for their better specific energy absorption than other honeycomb structures [30]. Andrea Bacigalupo et al. [31] numerically optimized the designs of metamaterial beams with tetrachiral, tetrachiral backhander, and hexachiral chiral honeycombs to obtain the best dispersion spectra, and used adaptive surrogate-based optimization technology to design a double quad-ligamentous chiral metamaterial with broadband low-frequency filtering [32]. Although beyond negative Poisson’s ratio, the chiral honeycomb metamaterials can effectively reduce vibration at locally resonant bandgaps, most designs use composite materials for complex configurations, resulting in localized resonance by adding extra mass, which leads to higher manufacturing costs, narrower bandgap width, and difficulty in flexible adjustment. As such, a design using a single material and integrated configuration is of significance to engineering applications. In addition, the metamaterial design methodology to satisfy both static and dynamic performance criterion is significant.
In this paper, inspired by the tetrachiral honeycomb and cross chiral structures, a novel cross tetrachiral honeycomb metamaterial with negative Poisson’s ratio is presented that is capable of wave isolation. Different from the usual approach of metamaterials to form a localized resonant bandgap by introducing additional mass, the cross tetrachiral honycomb metamaterial forms a wide bandgap by introducing more geometrical parameters including tile angle θ. The metamaterial can be geometrically and mechanically tuned to adjust the broadband bandgap while ensuring excellent static properties. By introducing more geometric parameters, the metamaterial can adjust the broadband bandgap through geometric regulation and mechanical regulation, while ensuring an excellent static performance. The remainder of the paper is arranged as follows. In Section 2, the effective elastic parameters of the proposed metamaterials are theoretically analyzed. In Section 3, three dimensionless quantities are introduced to describe the designability of static machinal response. In Section 4, the dynamic performance of the proposed metamaterial is studied through band structure analysis. Bandgaps have been designed through the selection of geometric parameters and can be tuned by external loading. In Section 5, conclusions are provided to summarize the design and analysis of cross tetrachiral honeycomb metamaterials.

2. Cross Tetrachiral Honeycomb Metamaterials

2.1. Model

In order to improve the deformability of tetrachiral honeycomb structures, a cross-chiral honeycomb structure [33,34,35] with folded corners, as shown in Figure 1a, is introduced into the classical tetrachiral honeycomb structure with the unit cell shown in Figure 1b. For classical traditional tetrachiral honeycomb structures, there is a circular ring at the center of the ligaments to achieve a light weight structure, but it also reduces the stiffness of the unit cell. This weakness affects the negative Poisson’s ratio and the structural stability in large deformations. To overcome this drawback, we replace the circular ring with a solid cylinder, and attach four identical ligaments tangentially to the cylinder to form an elementary unit, as marked by the blue dashed rectangular box in Figure 1c. By translating and mirroring the elementary unit, a unit cell of cross tetrachiral honeycomb metamaterial is formed, as shown in Figure 1c. The geometrical configuration of the unit cell can be described by the length of the elementary unit l, the tilt angle of the ligament θ, the radius of the solid cylinder r, the in-plane thickness t, and the cross-section width b of the ligaments, as shown in Figure 1c.

2.2. Effective Elastic Parameters

The symmetry of the structure formed from the unit cell (Figure 1c) indicates that the mechanical properties of the unit cell along the x and y directions are identical. Therefore, only the elastic properties in one direction (either x or y) need to be analyzed to describe the structure’s behavior.
Here, we consider small deformations consistent with linear elasticity, a rigid solid cylinder, and the ligaments analyzed using Euler–Bernoulli beam theory, as well as assuming that the deformations are in plane. It is worth noting that when the metamaterial undergoes an overall large deformation, it is actually a result of the cumulative effect of the bending deformations of the ligaments and the rotations of the cylinders. The strain generated locally in each ligament remains as a small deformation that satisfies the Euler–Bernoulli theory. Based on the geometric symmetry, when the unit cell is subject to a uniformly distributed load in the x or y directions obtained by an action on boundary points A1~A4 or B1~B4, respectively, only surface translations are observed. Suppose the unit cell is loaded in the y direction and eliminates the rigid body displacement of the center of the unit cell. In order to simplify the analysis, only one elementary unit in the unit cell in Figure 1c is considered, as shown in Figure 2.
The loading of the elementary unit is shown in Figure 2a, where F1 and F2 are externally applied loads along the y direction, respectively, on endpoints A and C; M1, M2, M3, and M4 are moments applied, respectively, on endpoints A, B, C, and D. Because the ligaments are tangentially attached to the cylinder, its lengths of two sides are not equal and can be separately named l11 and l1, as shown in Figure 2b.
Using the center point O of the solid cylinder as the fixed reference point with the assumed displacements u O = v O = 0 where u denotes displacement in the x direction and v denotes the displacement in the y direction, under the action of the external loads, the solid cylinder will have a rigid body rotation with a rotation angle φ around point O.
According to the force equilibrium along the y direction, we obtain F 1 = F 2 , and through the symmetry, the moment has the relations of M 1 = M 2 and M 3 = M 4 . And then, the moment equilibrium equations in plane can be derived using the following equation:
F 1 l 11 sin θ + r cos θ = M 1 + M 4 .
It is clear that the deformations of ligaments AA’ and CC’ are the same due to the symmetry, and so are the ligaments BB’ and DD’. Therefore, the deformation of the elementary unit can be divided into two parts for analysis. The deformation of ligament AA’, as shown in Figure 3, is first taken into consideration.
Under the actions of force F1 and moment M1, the displacements of ligament AA’ include axial tension/compression and bending deformations, as well as the rotation φ of the rigid cylinder. According to Figure 2b, the lengths of ligament AA’, named l11 and l1, have the following relations:
l 1 = 2 r sin θ + l 2 cos θ = r q tan θ + 1 2 cos θ l = α l ,
and
l 11 = l 1 r sin arccos 1 b r b tan θ = α r q sin arccos 1 b q r q b q tan θ l = β l ,
where α ( r q , θ ) and β r q , b q , θ can be represented by the dimensionless parameters r q = r l and b q = b l with the length of elementary unit l. If we set l 11 to be the equivalent length of the ligament, the displacements at point A can be expressed as follows:
u A = F 1 l 11 3 sin θ cos θ 3 E I M 1 l 11 2 cos θ 2 E I F 1 l 11 sin θ cos θ E A + φ l 0 cos θ 1 ,
and
v A = F 1 l 11 3 sin 2 θ 3 E I M 1 l 11 2 sin θ 2 E I + F 1 l 11 cos 2 θ E A + φ l 0 sin θ 1 ,
where E is the elastic modulus, I = t b 3 / 12 is the moment of inertia of the ligament, A = t b is the area of the cross-section of the ligament, and l 0 is the length from the center O of the cylinder to the endpoint A to be l 0 = l 1 2 + r 2 . As shown in Figure 4a, the rotation of the cylinder can cause the displacement of the ligament at the endpoint A to be φ l 0 . If we introduce a new tilt angle θ 1 between OA and ligament AA’, as shown in Figure 4b, its relation with θ is θ 1 = θ + arctan r l 1 = θ + arctan r q α .
The rotation of endpoint A should be equal to 0 due to the symmetric conditions, i.e.,
F 1 l 11 2 sin θ 2 E I M 1 l 11 E I + φ = 0 .
Then, the displacements at point A can be expressed in terms of external load F 1 , as follows:
u A = F 1 l 11 r cos θ 2 E I l 0 cos θ 1 l 11 cos θ 2 + l 11 2 sin θ 4 E I l 0 cos θ 1 l 11 cos θ 6 l 11 sin θ cos θ E A ,
and
v A = F 1 l 11 r cos θ 2 E I l 0 sin θ 1 l 11 sin θ 2 + l 11 2 sin θ 4 E I l 0 sin θ 1 l 11 sin θ 6 + l 11 cos 2 θ E A .
Similarly, for the deformation of ligament BB’, as illustrated in Figure 5, the displacements of endpoint B can be obtained using the following equations:
u B = M 4 l 11 2 sin θ 2 E I φ l 0 sin θ 1 ,
v B = M 4 l 11 2 cos θ 2 E I + φ l 0 cos θ 1 .
Considering Equations (1) and (6) and the symmetry conditions at the endpoint B with M 2 l 11 E I φ = 0 , the displacements of endpoint B can be expressed by the external load F 1 as an independent variable, i.e.,
u B = F 1 l 11 2 E I r cos θ + l 11 sin θ 2 l 0 sin θ 1 l 11 sin θ 2 ,
v B = F 1 l 1 2 E I r cos θ + l 1 sin θ 2 l 0 cos θ 1 l 1 cos θ 2 .
Thus, combining the displacements of endpoints A and B in Equations (9)–(12), for the elementary unit, the equivalent Poisson’s ratio and elastic modulus can be obtained using the following equations:
ν q = 2 u B l / 2 v A l ,
and
E s = 2 F 1 l t / 4 v A l .
Furthermore, we can define the equivalent Poisson’s ratio and elastic modulus in a dimensionless expression as follows:
ν q r q , b q , θ = r q cos θ + β sin θ 2 α 2 + r q 2 sin θ 1 β sin θ 2 r q cos θ α 2 + r q 2 sin θ 1 β sin θ 2 + β sin θ 2 α 2 + r q 2 sin θ 1 β sin θ 6 + b q 2 cos 2 θ 6 ,
and
E q r q , b q , θ = E s E = b q 3 2 β 6 r q cos θ α 2 + r q 2 sin θ 1 β sin θ 2 + 3 β sin θ α 2 + r q 2 sin θ 1 β sin θ 6 + b q 2 cos 2 θ .
Therefore, elastic parameters, i.e., equivalent Poisson’s ratio and elastic modulus, can be determined by three geometric dimensionless quantities—bq, rq, and θ—and the elastic properties of the whole structure along the x or y axes can be determined using Equations (15) and (16) by considering the symmetrical properties.

3. Static Performance

In order to verify the effectivity and accuracy of the equivalent elastic properties calculated using Equations (15) and (16), the finite element method is introduced to provide numerical simulation results for comparison. If the unit cell in Figure 1c is made of aluminum, its finite element model is shown in Figure 6a with the material properties and geometric parameters listed in Table 1. The unit cell produces 1% strain in the simulation. For the unit cell subjected to extension loading along the y direction, as shown in Figure 1c, its deformation can be numerically calculated by using the commercial software COMSOL Multiphysics to obtain the equivalent elastic properties. The finite element model in Figure 6a is meshed by the second-order tetrahedral element with size 1.2 × 10−2 mm. Among the geometric parameters in Equations (15) and (16), the tilt angle of ligament θ is more flexible to alter the equivalent elastic properties. Therefore, to change the tilt angle of the ligament θ from 5° to 35°, the varied equivalent Poisson’s ratios and elastic moduli are shown in Figure 6b,c, respectively.
For the equivalent Poisson’s ratios in Figure 6b, the analytical results (blue line) obtained using Equation (15) match very well with the numerical results (red square). The relative errors between them are shown in Figure 6d, which obtain a maximum of about 2% at a tilt angle θ = 18 ° , demonstrating the high accuracy of Equation (15). For the equivalent elastic modulus in Figure 6c and relative errors in Figure 6e, the analytical results from Equation (16) are less than the analytical results at small tilt angles, but still obtain a high accuracy (relative errors of less than 5%) when the tilt angles are larger than θ = 16 ° . Therefore, the proposed analytical equations can be used to theoretically predict the static mechanical properties of cross tetrachiral honeycomb metamaterials. It is noteworthy that the proposed cross tetrachiral honeycomb metamaterial has a good performance of negative Poisson’s ratio, which can be adjusted by altering the tilt angle θ.
In Figure 6b, the absolute value of equivalent negative Poisson’s ratio decreases as the tilt angle θ increases. As shown in Figure 2a, the small tilt angle of ligament θ can introduce a smaller vertical deformation along the tensile direction caused by ligaments AA’ and CC’, and the deformation can rotate the cylinder to induce a bending moment acting on ligaments BB’ and DD’ for providing horizontal deformation. Therefore, the equivalent negative Poisson’s ratio depends on the rotation of the cylinder to turn the vertical and horizontal displacements of deformed ligaments. With the gradually increasing tilt angle of the ligament, the bending deformations of ligaments AA’ and CC’ increase to cause a larger vertical displacement of the unit cell; however, the influence of the cylinder rotation on the horizonal displacement caused by ligaments BB’ and DD’ is limited. Therefore, the equivalent negative Poisson’s ratio can be adjusted by selecting a proper tilt angle θ.
For the equivalent elastic modulus, as shown in Figure 6c, there is a better consistency between the analytical results from Equation (16) and the numerical results at larger tilt angles. As the tilt angle θ increases, the equivalent elastic modulus becomes lower. This is because for a small tilt angle θ, the vertical displacement is mostly contributed by the longitudinally tensile deformations of ligaments AA’ and CC’, for which their stiffness is much larger than their bending stiffness, so that the equivalent elastic modulus is relatively big. On the contrary, the larger tilt angle can induce the vertical displacement with the bending-dominated deformations of ligaments AA’ and CC’, resulting in a larger deformation producing a smaller equivalent elastic modulus.
From Figure 6d,e, it can be seen that there is a difference between the theoretical predictions and the simulation results for Poisson’s ratio and elastic modulus, especially for the elastic modulus at smaller tilt angles θ. When θ is small, the ligaments are mainly subjected to tensile or compressive deformations by external loads. This causes the deformation of the ligament to be more dominated by axial deformation than bending deformation. However, the theoretical prediction is based on beam assumption, and only focuses on the bending deformation, neglecting the axial deformation. This is the reason why an error is caused between the theory and the simulation results at small tilt angles θ, while the axial deformation of the ligament is important and non-negligible. In addition, it demonstrates that the tilt angle θ is a useful parameter to change the equivalent static properties of cross tetrachiral honeycomb metamaterials.
Since the analytical results are in good agreement with the numerical results, the equivalent elastic parameters Eq and νq of cross tetrachiral honeycomb metamaterials can be predicted using Equations (15) and (16). There are three independent dimensionless geometric parameters—rq, bq, and θ—which can be used to design the equivalent elastic parameters of the proposed metamaterial, as shown in Figure 7. If we set the tilt angle θ = 8°, Eq and νq can be altered by changing geometric parameters rq and bq, as shown in Figure 7a,b, respectively. Similarly, if we set rq = 0.1 or bq = 0.05 and change the other two parameters, Eq and νq will correspondingly vary, as shown in Figure 7c,d or Figure 7e,f. Therefore, the static mechanical properties of cross tetrachiral honeycomb metamaterials can be predicted by selecting suitable geometric parameters.
Furthermore, it should be noticed that there is a stress concentration at the folding angle of two adjacent ligaments. In order to alleviate the stress concentration, the sharp folding angle is optimized into an arc curve with small curvature, as shown in Figure 8a. The stress distributions around the folding angle are calculated by using the commercial software COMSOL Mutiphysics 6.0, as shown in Figure 8b,c. The geometric parameters used in Figure 8a are the same as those in Table 1 and r1 = 20 mm. The unit cells in Figure 8b,c are subjected to a 1% tensile strain along the vertical direction. From Figure 8b,c, with the optimization, the maximum stress can be reduced from 166 MPa to 157 MPa with a reduction ratio of the maximum stress of about 5.42%, but the equivalent elastic modulus and Poisson’s ratio only change by about 1.05% and 0.24%, respectively. The results indicate that the small curvature at the folding angle of two adjacent ligaments does not make a big difference in the equivalent elastic parameters, but it can reduce the stress concentration at the folding angle to a certain extent. Actually, the folding angle curvature is closer to a straight line for meeting the requirements of machining in practical engineering, and can make the manufacturing process easier.

4. Dynamic Performance

As chiral metamaterials are highly applied in dynamic environments, it is very significant to investigate and adjust their dynamic properties. For the classical tetrachiral honeycomb metamaterial, as shown in Figure 1b, it is hard to design bandgaps in its band structure without any additional mass. Even if the metamaterial is made to produce a locally resonant bandgap by adding additional mass, the bandwidth will be narrow, and it will also be difficult to conveniently tune the bandgap. Therefore, how to manipulate the dynamic performance of cross tetrachiral honeycomb metamaterials is important to explore its potential in engineering applications.

4.1. Band Structures

For a linear elastic solid, its governing equation of wave motion with neglected body forces is as follows:
λ + μ u + μ 2 u = ρ u ¨ ,
where ρ is the mass density, λ and μ are the Lame constants, and u is the displacement vector, respectively. Considering a harmonic wave with angular frequency ω propagating in a two-dimensional Oxy plane, the solution of Equation (17) can be written as follows:
u = u k e i k r ω t ,
where r = x , y , k = k x , k y is the wave vector, and u k is the Bloch displacement vector.
By substituting Equation (18) into Equation (17), the wave propagation equation in the periodic structure can be rewritten as follows:
K λ , μ , k ω M ρ u k = 0 .
To obtain a non-trivial solution of u k from Equation (19), an eigenvalue problem should be solved, i.e.,
K λ , μ , k ω M ρ = 0
and the band structures of the periodic structure can be obtained.
For the same unit cell considered in a static case with the material properties and geometric parameters in Table 1, the dynamic performance of the cross tetrachiral honeycomb metamaterial is calculated using Equation (20) through the finite element commercial software COMSOL Multiphysics. To ensure the accuracy of the dynamic simulation, the finite element model is meshed by the second-order tetrahedral element with size 1.2 × 10−2 mm, which is less than one-tenth of the wavelength. Consider the in-plane case as the unit cell and its irreducible Brillouin zone shown in Figure 9a; if the tilt angle of ligament is selected as θ = 8 ° , the band structures can be calculated and are shown in Figure 9b. For the dynamic case, even if the metamaterial undergoes an overall large deformation, the end part of the metamaterial retains a small deformation to satisfy the assumptions of the Euler–Bernoulli theory. In the simulation, during the unit cell with a 2% strain, the strain of each ligament is less than 0.56%.
As shown in Figure 9b, there is a wide bandgap in the band structures from a frequency of 10.8 kHz to 13.4 kHz, as marked in the blue shade. To verify the band structures, a two-dimensional cross tetrachiral honeycomb structure with 13 × 6 unit cells, as shown in Figure 9c, is introduced to calculate the transmission ratio. A longitudinal plane wave is excited along the left side with the displacement amplitude U i , and the transmission signal along the right side with displacement amplitude U i is calculated using the finite element commercial software COMSOL Multiphysics. Perfect-matched layers (PMLs) are attached to both sides of the structure to eliminate the influence of reflected waves. The transmission ratio of the wave energy is calculated by U t 2 / U i 2 , and is shown in Figure 9d. The blue shade in Figure 9d marks the frequency range of the bandgap in Figure 9b, where a sharp drop in transmission ratios is observed. At other frequencies, the transmitted waves can propagate with a higher transmission ratio, which is consistent with the band structures.
In order to analyze the generation of the bandgap, the wave modes at the up and down boundaries of the bandgap with a frequency of 10.0 kHz and 13.6 kHz (marked 1 and 2 at Γ in Figure 9b) are calculated using COMSOL, and the displacement amplitude of the unit cell is amplified and shown in Figure 10. In Figure 10, it can be noticed that the cylinders rotate rigidly and the ligaments move like cantilever beams with in-plane flexural deformations. Therefore, our simplified model in Section 2 can also be used in analyzing the dynamic case. The longitudinally incident wave propagates through the cross tetrachiral honeycomb structure to induce the flexural deformations of horizontal ligaments, which rotate cylinders to drive the flexural deformations of vertical ligaments; then, the wave moment transfers from one to another.
Comparing the mode shapes in Figure 10a,b, there is no big difference to be found in the vertical ligaments, but there is a phase difference π between the horizontal ligaments. The lowest bandgap is formed by a combination of locally resonant bandgaps, as shown in Figure 10. In Figure 10a, if we focus on one elementary unit, the flexural deformations of the horizontal and vertical ligaments are all following the rotation direction of the cylinder to present an easier excited wave mode at a lower frequency (10.0 kHz). However, for the elementary unit in Figure 10b, an opposite flexural deformation between the horizontal and vertical ligaments can be observed, which is a more difficult mode to be excited and has a higher frequency (13.6 kHz). As a result, a bandgap is generated between both wave modes.

4.2. Bandgap Influenced by Geometric Parameters

To further explore the possibility of designing a cross tetrachiral honeycomb structure with specific dynamic properties, the influences of the geometric parameters of the unit cell are investigated one by one. Firstly, when the tilt angle of ligament θ is designed as θ = 5 ° , 6 ° , 12 ° , and 35 ° , the band structures are calculated, respectively, as shown in Figure 11a–d.
From Figure 11, it can be seen that there is no bandgap at θ = 5 ° . The bandgap generates at θ = 6 ° . Then, the bandwidth increases with the increasing θ until reaching 6.3 kHz (from 11.1 kHz to 17.4 kHz) at θ = 12 ° . After that, the bandwidth will decrease. Based on the band structures in Figure 11, the frequency of the mode shape at the upper boundary is the main reason for the change in bandgap, which indicates that the mode shape of Figure 10b is greatly affected by the tilt angle of the ligament. In contrast, the frequency of the lower boundary is relatively stable, which varies slightly with increasing θ. Therefore, the bandgap of the cross tetrachiral honeycomb structure can be adjusted by changing the tilt angle of the ligament.
Although the bandgap of the cross tetrachiral honeycomb structure is not generated by the local resonance of the cylinder, the change in cylinder radius can also vary the equivalent elastic modulus, as shown in Figure 7e; therefore, it also has an influence on the dynamic performance. For the same unit cell shown in Figure 6a, if we fix the tilt angle of the ligament to θ = 8 ° and change the radius r of the cylinder from 1.9 mm to 4.5 mm, the band structures can be calculated and they are separately shown in Figure 12a,b. Comparing to the case shown in Figure 9a for r = 3 mm with bandwidth 2.6 kHz, the smaller radius r = 1.9 mm can induce a narrower bandgap with bandwidth 0.2 kHz, as shown in Figure 12a, and the larger one r = 4.5 mm can provide a wider bandgap with bandwidth 4.6 kHz, as shown in Figure 12b. In addition, as shown in Figure 7e, the equivalent elastic modulus will decrease when the radius r of the cylinder becomes larger, so the lower frequency of the bandgap will accordingly become lower from 13.0 kHz (r = 1.9 mm) to 9.5 kHz (r = 4.5 mm). However, the upper frequency of the bandgap is more stable, growing from 13.0 kHz to 14.1 kHz; therefore, the bigger radius r of the cylinder provides a wider bandgap. Nevertheless, for the radius r = 4.5 mm, if the tilt angle of the ligament θ is changed from 8 ° to 5 ° , the band structures, as shown in Figure 12c, demonstrate that the bandwidth of the bandgap is influenced by the tilt angle of the ligament θ and the radius r together.
Since the tilt angle of the ligament θ and the radius of the cylinder r both have great influence on modulating the band structures, the relationships between the bandgap and selected geometric parameters should be investigated thoroughly. If the radius is fixed as r = 3 mm, the bandgap changed with the tilt angle of the ligament θ, as shown in Figure 13a. The bandgap does not exist at small tilt angles ( θ < 8 ° ). It should be noted that the relationship between bandwidth and tilt angle are not monotonic and there is a maximum bandwidth (6.3 kHz) at a certain angle ( θ = 12 ° ). The change in bandwidth is strongly dependent on the upper frequency of the bandgap with a rapid increase and gradual decrease, while the lower frequency of the bandgap alters smoothly. If the tilt angle of the ligament is fixed as θ = 8 ° , the bandgap change with radius r is shown in Figure 13b. With the increase in cylinder radius r, the bandwidth also continuously increases. When r = 1.9 mm, the bandwidth is 0.2 kHz, and when the radius r reaches 4.5 mm, the bandwidth increases to 4.6 kHz. Therefore, the cross tetrachiral honeycomb metamaterial is very sensitive to the changes in geometric parameters, and its bandgap can be adjusted by the proper selection of these parameters.

4.3. Bandgap Changed by External Load

In addition to the adjustment of the geometric parameters, the external load can induce a large deformation of the cross tetrachiral honeycomb structure, and can also cause a change in the dynamic performance. This characteristic will be more significant in engineering because it does not need to change the structural configuration and can meet a variety of engineering requirements in dynamic behaviors of cross tetrachiral honeycomb structures.
For the same model shown in Figure 9a, when vertical compression acts on the unit cell to induce a −2% strain, as shown in Figure 14a, its band structure is calculated and shown in Figure 14b. Comparing with the band structure in Figure 9b without an external load, the bandgap can be changed from 10.8~13.4 kHz to 11.1~16.9 kHz, and the bandwidth correspondingly varies from 2.6 kHz to 5.8 kHz, which is 3.2 kHz wider than that without external load. This is because the compressive deformation can increase the tilt angle of the vertical ligament of the unit cell. Accordingly, the mode shape at the upper frequency of the bandgap, as marked 1 in Figure 14b, is calculated and is shown in Figure 14c; it has the same mode shape as Figure 10b. This means that no matter if the external compression exists or not, the mechanism of producing the bandgap is the same, but the frequency and bandwidth of the bandgap can be changed.
If the vertical tension is acting on the unit cell and induces a 2% strain, as shown in Figure 14d, its band structure is calculated and is shown in Figure 14e. It is clear in Figure 14e that there is no bandgap, because the tensive deformation can make the tilt angle of the vertical ligament become smaller. For comparison with the bandgap in Figure 9b, the mode shape marked as 2 in Figure 14e is calculated and shown in Figure 14f. It can be noticed that the mode shape in Figure 14f is the same as Figure 10a, which is the mode shape of the lower frequency of the bandgap without external tension. This demonstrates that the mechanism of bandgap modulation through deformation is the same as that through changing the tilt angle of the ligament.
From Figure 14, it can be seen that the external load can effectively change the band structure and bandwidth of the bandgap, and can even convert the bandgap into the passband. In order to investigate the change in bandwidth with external load, the relationship between strain and bandgap is studied and is shown in Figure 15. In Figure 15, the results illustrate that with the increase in tensive strain, the bandgap gradually decreases and becomes narrower until it becomes passband at 2% strain; however, when the compressive strain works, the bandgap increases and becomes wider in order to reach the widest bandwidth at nearly 7 kHz. This confirms that the dynamic performance of the cross tetrachiral honeycomb structure has a very high sensitivity to external load.
In combination with the geometry of the unit cell after compression or tension in Figure 14a or Figure 14d, the reasons for the change in bandgap with strain can be analyzed. When the structure is vertically compressed, as in Figure 14a, the adjacent cylinders rotate harder in relative directions due to the external load, which makes the ligament between the cylinders more curved to form an arch. As a result, more energy is needed to stimulate the ligaments to produce the conversion of mode shape in Figure 14c, so the bandgap becomes wider at this time. When the structure is vertically stretched, the external load causes the adjacent cylinders to rotate oppositely, so that the ligaments between the cylinders become straight and gradually tend towards a straight line. As a result, less energy is required to stimulate the ligaments to keep the same mode shape as in Figure 14f, so the bandgap becomes narrower until it becomes a passband, which is a peculiar property of this structure.

5. Conclusions

A novel cross tetrachiral honeycomb metamaterial is proposed in this paper, and the combination of the cross chiral honeycomb structure and the tetrachiral honeycomb structure can be easily designed to reach the desired static and dynamic properties, particularly to achieve both excellent functions of negative Poisson’s ratio and wave isolation. Compared to the existing chiral honeycomb structures, the cross tetrachiral honeycomb metamaterial can break through the limitations of narrower and higher frequency bandgaps by introducing more tunable geometric parameters, and the integral design without additional masses can make the manufacturing process easy and precise.
Through theoretical analysis, the equivalent elastic modulus and Poisson’s ratio can be accurately predicted. In addition, these equivalent elastic parameters can be determined and adjusted by three geometric dimensionless parameters, which offer a significant convenience for the designable static performance of cross tetrachiral honeycomb structures in potential engineering applications.
The investigation of the dynamic properties of cross tetrachiral honeycomb metamaterials provides a tunable way to alter the bandgap for wave isolation. The influence of geometric parameters on the bandgap is studied, and the results indicate that the selection of the tilt angle of the ligament and the cylinder radius are important for adjusting the bandgap. Additionally, the external loading can also effect the change in bandgap from nonexistence to a wider one. With a small strain, the width of the bandgap can be changed greatly and linearly, and the bandgap can even be transformed into a passband. Therefore, the dynamic performance of cross tetrachiral honeycomb metamaterials is also designable.
In fact, the proposed metamaterial attempts to answer the research question about how to simultaneously satisfy both requirements of negative Poisson’s ratio and wide-frequency bandgap through a simple and low-cost design. Under some complex working conditions in engineering, the negative Poisson’s ratio and vibration isolation are required simultaneously.
In summary, the cross tetrachiral honeycomb metamaterial not only achieves a negative Poisson’s ratio effect, but can also result in excellent functions of vibration reduction and wave isolation. Therefore, the cross tetrachiral honeycomb metamaterial has designable static and dynamic performances by changing the geometric parameters or adjusting small linear deformations, and has broad prospects for engineering applications.

Author Contributions

Conceptualization, F.L. and S.S.; methodology, F.L. and W.W.; formal analysis, F.L. and S.S.; investigation, F.L. and S.S.; writing—original draft preparation, F.L.; writing—review and editing, S.S., R.X. and M.N.; visualization, F.L. and W.W.; supervision, Z.L. and S.S.; project administration, Z.L. and S.S.; funding acquisition, Z.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 numbers 11991033 and 12232001.

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.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Chen, Y.; Jia, Z.; Wang, L. Hierarchical honeycomb lattice metamaterials with improved thermal resistance and mechanical properties. Compos. Struct. 2016, 152, 395–402. [Google Scholar] [CrossRef]
  2. Walser, R.M. Electromagnetic metamaterials. In Complex Mediums II: Beyond Linear Isotropic Dielectrics, Proceedings of the International Symposium on Optical Science and Technology, San Diego, CA, USA, 29 July–3 August 2001; SPIE: Bellingham, WA, USA, 2001; Volume 4467, pp. 1–15. [Google Scholar]
  3. Lincoln, R.L.; Scarpa, F.; Ting, V.P.; Trask, R.S. Multifunctional composites: A metamaterial perspective. Multifunct. Mater. 2019, 2, 043001. [Google Scholar] [CrossRef]
  4. Airoldi, A.; Bettini, P.; Panichelli, P.; Oktem, M.F.; Sala, G. Chiral topologies for composite morphing structures—Part I: Development of a chiral rib for deformable airfoils. Phys. Status Solidi B 2015, 252, 1435–1445. [Google Scholar] [CrossRef]
  5. Jung, J.; Kim, H.-G.; Goo, S.; Chang, K.-J.; Wang, S. Realisation of a locally resonant metamaterial on the automobile panel structure to reduce noise radiation. Mech. Syst. Signal Process. 2019, 122, 206–231. [Google Scholar] [CrossRef]
  6. Akbari, M.; Shahbazzadeh, M.J.; La Spada, L.; Khajehzadeh, A. The graphene field effect transistor modeling based on an optimized ambipolar virtual source model for DNA detection. Appl. Sci. 2021, 11, 8114. [Google Scholar] [CrossRef]
  7. Liu, Y.; Xiong, F.; Yang, K.; Chen, Y. A novel omnidirectional self-locked energy absorption system inspired by windmill. J. Appl. Mech.-Trans. Asme 2020, 87, 085001. [Google Scholar] [CrossRef]
  8. Coulais, C.; Teomy, E.; de Reus, K.; Shokef, Y.; van Hecke, M. Combinatorial design of textured mechanical metamaterials. Nature 2016, 535, 529–532. [Google Scholar] [CrossRef]
  9. Fu, K.; Zhao, Z.; Jin, L. Programmable granular metamaterials for reusable energy absorption. Adv. Funct. Mater. 2019, 29, 1901258. [Google Scholar] [CrossRef]
  10. Pan, F.; Li, Y.; Li, Z.; Yang, J.; Liu, B.; Chen, Y. 3D pixel mechanical metamaterials. Adv. Mater. 2019, 31, 1900548. [Google Scholar] [CrossRef]
  11. Wenxuan, T.; Zhonglei, M.; Tiejun, C. Theory, experiment and applications of metamaterials. Sci. China-Phys. Mech. Astron. 2015, 58, 127001. [Google Scholar]
  12. Yu, X.; Zhou, J.; Liang, H.; Jiang, Z.; Wu, L. Mechanical metamaterials associated with stiffness, rigidity and compressibility: A brief review. Prog. Mater. Sci. 2018, 94, 114–173. [Google Scholar] [CrossRef]
  13. Zadpoor, A.A. Mechanical meta-materials. Mater. Horiz. 2016, 3, 371–381. [Google Scholar] [CrossRef]
  14. Bertoldi, K.; Vitelli, V.; Christensen, J.; van Hecke, M. Flexible mechanical metamaterials. Nat. Rev. Mater. 2017, 2, 17066. [Google Scholar] [CrossRef]
  15. Qu, J.; Gerber, A.; Mayer, F.; Kadic, M.; Wegener, M. Experiments on metamaterials with negative effective static compressibility. Phys. Rev. X 2017, 7, 041060. [Google Scholar] [CrossRef]
  16. Evans, J.S.O. Negative thermal expansion materials. J. Chem. Soc.-Dalton Trans. 1999, 19, 3317–3326. [Google Scholar] [CrossRef]
  17. Lipton, J.I.; MacCurdy, R.; Manchester, Z.; Chin, L.; Cellucci, D.; Rus, D. Handedness in shearing auxetics creates rigid and compliant structures. Science 2018, 360, 632–635. [Google Scholar] [CrossRef]
  18. Zhua, R.; Liub, X.N.; Hub, G.K.; Sunc, C.T.; Huang, G.L. A chiral elastic metamaterial beam for broadband vibration suppression. J. Sound Vib. 2014, 333, 2759–2773. [Google Scholar] [CrossRef]
  19. Novak, N.; Vesenjak, M.; Ren, Z. Auxetic cellular materials—A Review. Stroj. Vestn. J. Mech. Eng. 2016, 62, 485–493. [Google Scholar] [CrossRef]
  20. Lakes, R. Deformation mechanisms in negative poisson ratio materials—Structural aspects. J. Mater. Sci. 1991, 26, 2287–2292. [Google Scholar] [CrossRef]
  21. Grima, J.N.; Gatt, R.; Farrugia, P.-S. On the properties of auxetic meta-tetrachiral structures. Phys. Status Solidi B-Basic Solid State Phys. 2008, 245, 511–520. [Google Scholar] [CrossRef]
  22. Alderson, A.; Alderson, K.L.; Attard, D.; Evans, K.E.; Gatt, R.; Grima, J.N.; Miller, W.; Ravirala, N.; Smith, C.W.; Zied, K. Elastic constants of 3-, 4-and 6-connected chiral and anti-chiral honeycombs subject to uniaxial in-plane loading. Compos. Sci. Technol. 2010, 70, 1042–1048. [Google Scholar] [CrossRef]
  23. Yang, S.; Qi, C.; Wang, D.; Gao, R.; Hu, H.; Shu, J. A comparative study of ballistic resistance of sandwich panels with aluminum foam and auxetic honeycomb cores. Adv. Mech. Eng. 2013, 5, 589216. [Google Scholar] [CrossRef]
  24. Lakes, R. Foam structures with a negative poissons ratio. Science 1987, 235, 1038–1040. [Google Scholar] [CrossRef] [PubMed]
  25. May, P.R.A.; Fuster, J.M.; Haber, J. Woodpecker drilling behavior—endorsement of the rotational theory of impact brain injury. Arch. Neurol. 1979, 36, 370–373. [Google Scholar] [CrossRef] [PubMed]
  26. Xin, Y.; Wang, H.; Wang, C.; Cheng, S.; Zhao, Q.; Sun, Y.; Gao, H.; Ren, F. Properties and tunability of band gaps in innovative reentrant and star-shaped hybrid honeycomb metamaterials. Results Phys. 2021, 24, 104024. [Google Scholar] [CrossRef]
  27. Chen, L.; Guo, Y.; Yi, H. Optimization study of bandgaps properties for two-dimensional chiral phononic crystals base on lightweight design. Phys. Lett. A 2021, 388, 127054. [Google Scholar] [CrossRef]
  28. Liu, X.N.; Hu, G.K.; Sun, C.T.; Huang, G.L. Wave propagation characterization and design of two-dimensional elastic chiral metacomposite. J. Sound Vib. 2011, 330, 2536–2553. [Google Scholar] [CrossRef]
  29. Sun, P.; Guo, H.; Jin, F.; Zhang, Z.; Liu, N.; Yuan, T.; Ma, L.; Wang, Y. Mechanics and extreme low-frequency band gaps of auxetic hexachiral acoustic metamaterial with internal resonant unit. Appl. Acoust. 2022, 200, 109046. [Google Scholar] [CrossRef]
  30. Alomarah, A.; Masood, S.H.; Sbarski, I.; Faisal, B.; Gao, Z.; Ruan, D. Compressive properties of 3D printed auxetic structures: Experimental and numerical studies. Virtual Phys. Prototyp. 2020, 15, 1–21. [Google Scholar] [CrossRef]
  31. Bacigalupo, A.; Lepidi, M.; Gnecco, G.; Vadalà, F.; Gambarotta, L. Optimal design of the band structure for beam lattice metamaterials. Front. Mater. 2019, 6, 2. [Google Scholar] [CrossRef]
  32. Bacigalupo, A.; Gnecco, G.; Lepidi, M.; Gambarotta, L. Computational design of innovative mechanical metafilters via adaptive surrogate-based optimization. Comput. Methods Appl. Mech. Eng. 2021, 375, 113623. [Google Scholar] [CrossRef]
  33. Grima, J.N.; Cassar, R.N.; Gatt, R. On the effect of hydrostatic pressure on the auxetic character of NAT-type silicates. J. Non-Cryst. Solids 2009, 355, 1307–1312. [Google Scholar] [CrossRef]
  34. Smith, C.W.; Grima, J.N.; Evans, K.E. A novel mechanism for generating auxetic behaviour in reticulated foams: Missing rib foam model. Acta Mater. 2000, 48, 4349–4356. [Google Scholar] [CrossRef]
  35. Lu, Z.; Wang, Q.; Li, X.; Yang, Z. Elastic properties of two novel auxetic 3D cellular structures. Int. J. Solids Struct. 2017, 124, 46–56. [Google Scholar] [CrossRef]
Figure 1. (a) The cross chiral honeycomb structure; (b) tetrachiral honeycomb; (c) unit cell of cross tetrachiral honeycomb metamaterial.
Figure 1. (a) The cross chiral honeycomb structure; (b) tetrachiral honeycomb; (c) unit cell of cross tetrachiral honeycomb metamaterial.
Materials 17 04652 g001
Figure 2. (a) Force diagram of elementary unit; (b) geometric parameters elementary unit.
Figure 2. (a) Force diagram of elementary unit; (b) geometric parameters elementary unit.
Materials 17 04652 g002
Figure 3. Deformations of ligament AA’ under actions of external load F1 and moment M1.
Figure 3. Deformations of ligament AA’ under actions of external load F1 and moment M1.
Materials 17 04652 g003
Figure 4. Displacement of ligament AA’ induced by cylinder rotation: (a) schematic diagram; (b) geometric relations between deformed and undeformed ligaments.
Figure 4. Displacement of ligament AA’ induced by cylinder rotation: (a) schematic diagram; (b) geometric relations between deformed and undeformed ligaments.
Materials 17 04652 g004
Figure 5. Deformation of ligament BB’ under the action of moment M4.
Figure 5. Deformation of ligament BB’ under the action of moment M4.
Materials 17 04652 g005
Figure 6. (a) Finite element model of unit cell. The analytical and numerical results of equivalent (b) Poisson’s ratio and (c) elastic modulus. Relative errors between analytical and numerical results of equivalent (d) Poisson’s ratios and (e) elastic modulus.
Figure 6. (a) Finite element model of unit cell. The analytical and numerical results of equivalent (b) Poisson’s ratio and (c) elastic modulus. Relative errors between analytical and numerical results of equivalent (d) Poisson’s ratios and (e) elastic modulus.
Materials 17 04652 g006
Figure 7. (a) Diagram of Eq with respect to bq and rq when θ = 8 ° . (b) Diagram of υq with respect to bq and rq when θ = 8 ° . (c) Diagram of Eq with respect to bq and θ when rq = 0.1. (d) Diagram of υq with respect to bq and θ when rq = 0.1. (e) Diagram of Eq with respect to rq and θ when bq = 0.05. (f) Diagram of υq with respect to rq and θ when bq = 0.05.
Figure 7. (a) Diagram of Eq with respect to bq and rq when θ = 8 ° . (b) Diagram of υq with respect to bq and rq when θ = 8 ° . (c) Diagram of Eq with respect to bq and θ when rq = 0.1. (d) Diagram of υq with respect to bq and θ when rq = 0.1. (e) Diagram of Eq with respect to rq and θ when bq = 0.05. (f) Diagram of υq with respect to rq and θ when bq = 0.05.
Materials 17 04652 g007
Figure 8. (a) Optimization of the folding angle of two adjacent ligaments. (b) Stress cloud image of the unit cell without optimization. (c) Stress cloud image of the unit cell with optimization.
Figure 8. (a) Optimization of the folding angle of two adjacent ligaments. (b) Stress cloud image of the unit cell without optimization. (c) Stress cloud image of the unit cell with optimization.
Materials 17 04652 g008
Figure 9. (a) The unit cell and the Brillouin zone. (b) The band structures. (c) The cross tetrachiral honeycomb structure. (d) The transmission ratio of the in-plane propagating wave in the structure.
Figure 9. (a) The unit cell and the Brillouin zone. (b) The band structures. (c) The cross tetrachiral honeycomb structure. (d) The transmission ratio of the in-plane propagating wave in the structure.
Materials 17 04652 g009
Figure 10. The mode shapes of unit cell at the (a) down and (b) up boundaries of the bandgap.
Figure 10. The mode shapes of unit cell at the (a) down and (b) up boundaries of the bandgap.
Materials 17 04652 g010
Figure 11. Band structures of different tilt angles: (a) θ = 5 ° , (b) θ = 6 ° , (c) θ = 12 ° , (d) θ = 35 ° .
Figure 11. Band structures of different tilt angles: (a) θ = 5 ° , (b) θ = 6 ° , (c) θ = 12 ° , (d) θ = 35 ° .
Materials 17 04652 g011
Figure 12. Band structures of different cylinder radii and tilt angles: (a) r = 1.9 mm, θ = 8 ° . (b) r = 4.5 mm, θ = 8 ° . (c) r = 4.5 mm, θ = 5 ° .
Figure 12. Band structures of different cylinder radii and tilt angles: (a) r = 1.9 mm, θ = 8 ° . (b) r = 4.5 mm, θ = 8 ° . (c) r = 4.5 mm, θ = 5 ° .
Materials 17 04652 g012
Figure 13. Bandgaps change with (a) the tilt angle of ligament θ and (b) radius r.
Figure 13. Bandgaps change with (a) the tilt angle of ligament θ and (b) radius r.
Materials 17 04652 g013
Figure 14. (a) Unit cell compressed by 2% strain. (b) Band structure after compression. (c) Mode shape of marked 1 at Γ, frequency 16.9 kHz. (d) Unit cell stretched by 2% strain. (e) Band structure after tension. (f) Mode shape of marked 2 at Γ, frequency 10.2 kHz.
Figure 14. (a) Unit cell compressed by 2% strain. (b) Band structure after compression. (c) Mode shape of marked 1 at Γ, frequency 16.9 kHz. (d) Unit cell stretched by 2% strain. (e) Band structure after tension. (f) Mode shape of marked 2 at Γ, frequency 10.2 kHz.
Materials 17 04652 g014
Figure 15. Bandwidth of bandgap varying with strain.
Figure 15. Bandwidth of bandgap varying with strain.
Materials 17 04652 g015
Table 1. Material properties and geometric parameters.
Table 1. Material properties and geometric parameters.
E (GPa)υρ (kg/m3)t (m)b (m)l (m)r (m)A (m2)I (m4)
700.3327000.0010.0010.060.00310−6 8.33 × 10 14
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

Liu, F.; Shao, S.; Wang, W.; Xia, R.; Negahban, M.; Li, Z. A Novel Cross Tetrachiral Honeycomb Metamaterial with Designable Static and Dynamic Performances. Materials 2024, 17, 4652. https://doi.org/10.3390/ma17184652

AMA Style

Liu F, Shao S, Wang W, Xia R, Negahban M, Li Z. A Novel Cross Tetrachiral Honeycomb Metamaterial with Designable Static and Dynamic Performances. Materials. 2024; 17(18):4652. https://doi.org/10.3390/ma17184652

Chicago/Turabian Style

Liu, Fengming, Shixuan Shao, Weihan Wang, Rongyu Xia, Mehrdad Negahban, and Zheng Li. 2024. "A Novel Cross Tetrachiral Honeycomb Metamaterial with Designable Static and Dynamic Performances" Materials 17, no. 18: 4652. https://doi.org/10.3390/ma17184652

APA Style

Liu, F., Shao, S., Wang, W., Xia, R., Negahban, M., & Li, Z. (2024). A Novel Cross Tetrachiral Honeycomb Metamaterial with Designable Static and Dynamic Performances. Materials, 17(18), 4652. https://doi.org/10.3390/ma17184652

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