Next Article in Journal
Application of Dimming Compensation Technology Via Liquid Crystal Lens for Non-Imaging Projection Laser Systems
Next Article in Special Issue
Analysis of the Transmission Characteristic and Stress-Induced Birefringence of Hollow-Core Circular Photonic Crystal Fiber
Previous Article in Journal
Focusing of Microcrystals and Liquid Condensates in Acoustofluidics
Previous Article in Special Issue
Design of Polarization Splitter via Liquid and Ti Infiltrated Photonic Crystal Fiber
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Elastic Wave Propagation of Two-Dimensional Metamaterials Composed of Auxetic Star-Shaped Honeycomb Structures

1
Department of Mechanical Engineering, National Cheng Kung University, University Road, Tainan City 701, Taiwan
2
Department of Digital Design and Information Management Chung Hwa University of Medical Technology, Tainan City 717, Taiwan
*
Author to whom correspondence should be addressed.
Crystals 2019, 9(3), 121; https://doi.org/10.3390/cryst9030121
Submission received: 29 January 2019 / Revised: 18 February 2019 / Accepted: 25 February 2019 / Published: 26 February 2019
(This article belongs to the Special Issue Sonic and Photonic Crystals)

Abstract

:
In this paper, the wave propagation in phononic crystal composed of auxetic star-shaped honeycomb matrix with negative Poisson’s ratio is presented. Two types of inclusions with circular and rectangular cross sections are considered and the band structures of the phononic crystals are also obtained by the finite element method. The band structure of the phononic crystal is affected significantly by the auxeticity of the star-shaped honeycomb. Some other interesting findings are also presented, such as the negative refraction and the self-collimation. The present study demonstrates the potential applications of the star-shaped honeycomb in phononic crystals, such as vibration isolation and the elastic waveguide.

1. Introduction

Phononic crystals are artificial crystals composed of a periodic alternation of at last two different materials. These materials were first studied by Sigalas and Economou et al. [1] and Kushwaha et al. [2,3], in which one of the most important properties in the crystals, the band gap, was investigated. They found that for some periodic arrangement, the waves could be stopped by phononic crystals. This unique property was found to have potential to be applied in isolating incident waves, acoustics and vibrations in many studies [4,5,6]. Sigalas [7] studied the elastic wave propagation in phononic crystals by using the plane wave expansion method. They investigated the defect effects inside the crystal system. The analysis of the band gaps suggested that the wave propagation in phononic crystals can be controlled by the defect states. Other research regarding the defect effects on the wave propagation in phononic crystals was also presented in the literature [8,9]. These studies suggested that their design of defect modes can be used as high-Q narrow band pass acoustical filters.
The concept of locally resonant phononic crystal was first proposed by Liu et al. [10]. The materials, based on the simple realization and composited with locally resonant structural units, can exhibit effective negative elastic constants at certain frequency ranges in their study. The model and experiments showed the spectral gaps with lattice constants two orders of magnitude smaller than the relevant sonic wavelength. Their findings suggested that the locally resonant crystal is the key factor to control the wave propagation with large wave length by small-size phononic crystals. Liu [11] and Wang [12] performed systematic studies to analyze the factors that affect the locally resonant phononic crystals. They also found that the band gap can be manipulated by tuning the elastic modulus. The suitable range of the elastic modulus can be obtained through the design of the local structures in the periodic crystal [13,14,15]. For the locally resonant phononic crystals, the wave length is much larger than the lattice constant. Therefore, the crystals can be modelled as a homogeneous medium with effective elastic modulus and effective Poisson’s ratio [16,17].
Some of the above mentioned phononic crystals are the cellular structures, which have properties such as high stiffness and low density. Some cellular structures have unique properties, such as negative Poisson ratio. The material with the negative Poisson ratio is called the auxetics. Gibson et al. [18] investigated that some cellular structures, such as re-entrance honeycomb, have a negative effective Poisson ratio. Then, Almgren [19] proposed a 3D structure composed of honeycomb structure with a Poisson ratio of −1. The negative poisson’s ratios in composites with star-shaped inclusions by numerical homogenization approach was presented by Panagiotopoulus et al. [20]. The wave propagations in the auxetic materials were also investigated. Gonella and Ruzzene [21] considered the plane wave propagation in re-entrant and hexagonal honeycomb structures. Two-dimensional dispersion relations were analyzed to reveal peculiar properties of the re-entrant honeycomb structures. Liebold-Ribeiro and Körner [22] performed the eigenmode analysis to various periodic cellular materials. Then, the star-shaped honeycomb with varied Poisson ratio as presented by Meng et al. [23] and the wave propagation behaviors and the band gaps were also analyzed in their studies. Wojciechowski et al. [24] presented some investigations about the auxetics and other systems of anomalous characteristics. The abovementioned papers showed that the auxetic material with negative values of Poisson ratio have significant effects on the wave motion in elastic solid.
Self-collimation is another important phenomenon of wave propagation in phononic crystals. Qiu et al. [25] proposed a design of a highly-directional acoustic source, formed by placing a line acoustic source inside a planar cavity of two-dimensional phononic crystals. The propagation characteristics of elastic transverse waves emitted by line sources embedded inside two-dimensional phononic crystals were studied by Liu and Su [26]. Cicek and his co-workers [27,28] studied numerically the wave propagation characteristics and self-collimation phenomenon of phononic crystals. Besides this, some experimental results have been performed to validate this phenomenon [29,30,31].
Recently, we published some research regarding wave propagation in phononic crystals [32,33], in which we performed the wave propagation analysis and had some interesting findings, including self-collimation and negative reflection. In this paper, we propose a novel auxetic material composed of star-shaped honeycomb structure. The Poisson’s ratio with various negative value of the structure can be obtained by adjusting the star angle θ. The dispersion of the phononic crystal composed of auxetic star-shaped honeycomb will be analyzed by eigenmode analysis. The self-collimation characteristic of the star-shaped honeycomb is also investigated.

2. Materials and Methods

The dynamic equilibrium equations for an isotropic elastic solid are expressed as [34]:
( λ + μ ) u j , j i + μ u i , j j = ρ u ¨ i
where u i is the displacement vector component, λ and μ are the Lame constant, and ρ is the mass density, respectively. The Lame constants can be written in terms of Young’s modulus E and Poisson’s ratio ν , i.e., λ = E ν / ( 1 + ν ) ( 1 2 ν ) and μ = E / 2 ( 1 ν ) . The constitutive equation for a linear isotropic elastic material is denoted as:
τ i j = λ ε k k δ i j + 2 μ ε i j
in which, τ i j , ε i j , and δ i j are stress, strain tensor, and the unitary tensor, respectively. Then, the strain-displacement relation is as follows:
ε i j = 1 2 ( u i , j + u j , i )
The phononic crystal composed of square lattice of steel cylinder in an auxetic matrix is a periodic medium. The displacement field of the periodic structure can be expressed as follows by applying the Bloch theorem.
u ( r ) = u k ( r ) e i ( ω t k · r )
where u ( r ) is the displacement vector, ω is the natural frequency, r is the position vector, k is the wave vector in the first Brillouin zone, and u k ( r ) is the periodic vector function of r , respectively. The phononic crystal lattice and u k have the same periodicity. The wave propagation in the periodic structure can be solved by using Equations (1) and (2). In this paper, the solution procedures are done by the finite element method. The boundary conditions of the unit cell is:
u ( r + a ) = e i ( k · a ) u ( r )
Thus, we can obtain the discrete from of finite element eigenvalue system in the unit cell as following form by applying the boundary conditions:
( K ω 2 M ) u = 0
where K and M are stiffness and mass matrix of the unit cell and a is the lattice constant, respectively.

3. Results

In this section, numerical results will be presented to demonstrate the band structure of the phononic crystal composed of auxetic materials. The two types of phononic crystal systems in this paper are shown in Figure 1a,b, respectively. The lattice constant of phononic crystal is a = 0.8   m . The material used for the inclusion is steel with material properties E s = 200   G P a , ρ s = 7850   k g / m 3 and ν s = 0.33 . Two shapes of the cross sections of the inclusion, the circular shape with radius r = 0.3 a and rectangular shape with width w = 14 a / 15 and height h = a / 5 , are considered, as shown in Figure 1a. For rectangular shaped cross section as shown in Figure 1b, the incline angle ϕ of the inclusion with rectangular cross section is also considered as a parameter. Therefore, we can realize an auxetic material by using the star-shaped honeycomb structure. The effective Young’s modulus and effective Poisson’s ratio of the star-shaped honeycomb structure are considered as the material properties of the auxetic matrix.
Figure 2 presents band structures of the elastic wave in the phononic crystals with cylinder inclusions and auxetic matrix, of which three different Poisson’s ratios, ν = 0.5 , 0.7 , and 0.9 , are considered. The numerical results show that as the absolute value of the Poisson’s ratio, | ν | , increases, the center frequencies of the band gaps increase. The materials with negative Poisson’s ratio are rare in the natural materials, especially for high strength materials. Thus, the re-entrant and star-shaped honeycomb structures as shown in Figure 3 are proved to have negative values of Poisson’s ratio [23]. In this study, the shar-shaped honeycomb structure is adopted as the matrix in the phononic crystals. The scheme of the phononic crystal, composed of periodic steel inclusion with circular cross section and the star-shaped honeycomb matrix, is shown in Figure 3a. Figure 3b shows the unit cell of the star-shaped honeycomb, the material properties of which is steel, same as the material used for inclusions. The length of the ligament of the star is 0.05   m , the width of the beam is one-tenth of the length, and lattice constant of the star shaped unit cell is L = 0.16   m . Then, the effective material properties of the matrix are listed in Table 1 [35] and it can be observed that the Poisson’s ratio varies with star angle θ . For star angles smaller than 65 ° , the Poisson’s ratio is negative, and the matrix can be considered as an auxetic material.
Figure 4 shows the effects of the star angles on the band structures of the phononic crystal composed of steel cylinders in the star-shaped honeycomb matrix. Three different star angles, θ = 55 ° , 60 ° , and 70 ° are calculated and the corresponding effective Poisson’s ratios are 0.53 , 0.35 and 0.15 , respectively. The band structures of the phononic crystal system are computed and obtained by the finite element software C O M S O L ® . Three band gaps are found in the frequency range from 0   H z to 1.3   k H z . The frequency of the band gap increases if the absolute value of Poisson’s ratio decreases. Additionally, it is also observed that the widths of the band gap increase with the decrease of Poisson’s ratio. According to the figure, it is evident that the band structures strongly depend on the auxeticity of the matrix.
Then, the equi-frequency contour (EFC) of the third band of the star-shaped unit cell for the case of θ = 55 ° and ν = 0.53 is plotted in Figure 5. The dashed black line denotes the construction line and the red arrow indicates the gradient of the contour, which is the direction of the elastic wave propagation in the phononic crystal. It indicates that the negative refraction occurs at this band. The full wave finite element simulation result is shown in Figure 6. A plane elastic wave of 500   H z is incident from the left side of the phononic crystal composed of star-shaped honeycomb matrix with star angle θ = 55 ° . The width of the plane wave and the incident angle are 4 a and 30 ° , respectively. The simulation of the phononic crystal system is also computed and obtained by the finite element software C O M S O L ® .The negative refraction is shown clearly in Figure 6. The refraction angle is about 45 ° and the refraction ratio is about 0.707 .
Figure 7 presents the numerical results of wave propagation of the phononic crystal, composed of inclusion with rectangular shape and the star-shaped honeycomb as the matrix. The star angle is θ = 55 ° , and Poisson’s ratio of the matrix is ν = 0.53 , respectively. The geometric parameters of the inclusions are w = 14 a / 15 , h = a / 5 , and ϕ = 0 ° . Figure 7a is the band structure, in which the X M region exhibit approximately horizontal lines. The equi-frequency contours are shown in Figure 7b by taking the 8th band as an example. The frequency range of the 8th band is from 820   H z to 900   H z . In the case that the eigen-frequency is greater than 870   H z , the equi-frequency contours are more like the vertical straight lines. It suggests that the self-collimation phenomenon occurs with the wide range of the wave incident into the phononic crystal. After that, the simulation of the wave propagation for three different orientation angles of the rectangular inclusions is presented in Figure 8 and the finite element model is a 50 × 40 array of lattices. The incident wave with 885   H z is from the left boundary and the incident angle is fixed horizontally for all cases. It is obvious that the wave propagation direction aligns the orientation angle ϕ . In addition, Figure 9 shows that the simulation of the wave propagation for the model consists of gradually varied orientation angles. The orientation angle is ϕ = 0 ° at the left boundary, x = 0 , and linearly increases to ϕ = 50 ° at the right boundary, x   =   100 a . The simulation results reveal that when the wave is horizontally incident from the left boundary to the phononic crystal, the wave propagation direction is automatically collimated along the orientation angles of the rectangular inclusions, resulting in a curved path of the wave propagation. The numerical results shown in Figure 7 and Figure 8 validate the self-collimation phenomenon of the phononic crystal with star-shaped honeycomb matrix and rectangular inclusions.

4. Discussion and Conclusions

In this paper, the elastic wave propagations in phononic crystals composed of steel inclusions and auxetic material matrix are investigated. The star-shaped honeycomb structures are utilized and have been proven to have negative Poisson’s ratio. Two types of inclusions with circular and rectangular cross sections are considered and discussed. The band structures of the phononic crystals are calculated and obtained by the finite element method.
The auxeticity of the star-shaped honeycomb is found to have significant effect on the band structure of the phononic crystal. The bandwidths and the mid-frequencies of the band gap increase as the Poisson’s ratio decreases. In addition to this, the negative refraction is also found in the phononic crystals for the star-shaped honeycomb as the matrix. The self-collimation phenomenon is also studied. The wide angle and wide band collimation are observed in the phononic crystal with rectangular inclusions. The present study demonstrates the potential applications of the star-shaped honeycomb in phononic crystals, such as vibration isolation and the elastic waveguide.

Author Contributions

Formal analysis, S.-Y.C.; Methodology, S.-Y.C.; Resources, J.-Y.Y.; Supervision, L.-W.C.; Validation, J.-Y.Y.; Writing – review & editing, C.-D.C.

Acknowledgments

The financial support of Ministry of Science and Technology under the contract MOST 104-2221-E-006-113-MY3 is appreciated.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sigalas, M.M.; Economou, E.N. Elastic and Acoustic Wave Band Structure. J. Sound Vib. 1992, 158, 377–382. [Google Scholar] [CrossRef]
  2. Kushwaha, M.S.; Halevi, P.; Dobrzynski, L.; Djafari-Rouhani, B. Acoustic Band Structure of Periodic Elastic Composites. Phys. Rev. Lett. 1993, 71, 2022. [Google Scholar] [CrossRef] [PubMed]
  3. Kushwaha, M.S.; Halevi, P.; Martinez, G.; Dobrzynski, L.; Djafari-Rouhani, B. Theory of Acoustic Band Structure of Periodic Elastic Composites. Phys. Rev. B 1994, 49, 2313. [Google Scholar] [CrossRef]
  4. Martinezsala, R.; Sancho, J.; Sanchez, J.V.; Gomez, V.; Llinares, J.; Meseguer, F. Sound Attenuation by Sculpture. Nature 1995, 378, 241. [Google Scholar] [CrossRef]
  5. Thomas, E.L.; Gorishnyy, T.; Maldovan, M. Colloidal Crystals Go Hypersonic. Nat. Mater. 2006, 5, 773–774. [Google Scholar] [CrossRef] [PubMed]
  6. Wen, J.; Wang, G.; Yu, D.; Zhao, H.; Liu, Y. Theoretical and Experimental Investigation of Flexural Wave Propagation in Straight Beams with Periodic Structures: Application to a vibration isolation structure. J. Appl. Phys. 2005, 97, 114907. [Google Scholar] [CrossRef]
  7. Sigalas, M. Elastic Wave Band Gaps and Defect States in Two-dimensional Composites. J. Acous. Soc. Am. 1997, 101, 1256. [Google Scholar] [CrossRef]
  8. Wu, F.; Hou, Z.; Liu, Z.; Liu, Y. Point Defect States in Two-dimensional Phononic Crystals. Phys. Lett. A 2001, 292, 198–202. [Google Scholar] [CrossRef]
  9. Zhang, X.; Liu, Z.; Liu, Y.; Wu, F. Defect States in 2D Acoustic Band-gap Materials with Bend-shaped Linear Defects. Solid State Commun. 2004, 130, 67–71. [Google Scholar] [CrossRef]
  10. Liu, Z.; Zhang, X.; Mao, Y.; Zhu, Y.; Yang, Z.; Chan, C.T.; Sheng, P. Locally Resonant Sonic Materials. Science 2000, 289, 1734–1736. [Google Scholar] [CrossRef] [PubMed]
  11. Liu, Z.; Chan, C.T.; Sheng, P. Three-component Elastic Wave Band-gap Material. Phys. Rev. B 2002, 65, 165116. [Google Scholar] [CrossRef]
  12. Gang, W.; Li-Hui, S.; Yao-Zong, L.; Ji-Hong, W. Accurate Evaluation of Lowest Band Gaps in Ternary Locally Resonant Phononic Crystals. Chin. Phys. 2006, 15, 1843. [Google Scholar] [CrossRef]
  13. Larabi, H.; Pennec, Y.; Djafari-Rouhani, B.; Vasseur, J. Multicoaxial Cylindrical Inclusions in Locally Resonant Phononic Crystals. Phys. Rev. E 2007, 75, 066601. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, S.; Cheng, J. Existence of Broad Acoustic Bandgaps in Three-Component Composite. Phys. Rev. B 2003, 68, 245101. [Google Scholar] [CrossRef]
  15. Zhang, X.; Liu, Y.; Wu, F.; Liu, Z. Large Two-dimensional Band Gaps in Three-component Phononic Crystals. Phys. Lett. A 2003, 317, 144–149. [Google Scholar] [CrossRef]
  16. Fok, L.; Zhang, X. Negative Acoustic Index Metamaterial. Phys. Rev. B 2011, 83, 214304. [Google Scholar] [CrossRef]
  17. Wu, Y.; Lai, Y.; Zhang, Z.-Q. Elastic Metamaterials with Simultaneously Negative Effective Shear Modulus and Mass Density. Phys. Rev. Lett. 2011, 107, 105506. [Google Scholar] [CrossRef] [PubMed]
  18. Gibson, L.J.; Ashby, M.F.; Schajer, G.; Robertson, C.I. The Mechanics of Two-Dimensional Cellular Materials. Proc. R. Soc. Lond. A 1982, 382, 25–42. [Google Scholar] [CrossRef]
  19. Robert, F.A. An Isotropic Three-dimensional Structure with Poisson’s Ratio = −1. J. Elast. 1985, 15, 427. [Google Scholar]
  20. Theocaris, P.S.; Stavroulakis, G.E.; Panagiotopoulus, P.D. Negative Poisson’s ratios in composites with star-shaped inclusions: A numerical homogenization approach. Arch. Appl. Mech. 1997, 67, 274–286. [Google Scholar] [CrossRef]
  21. Gonella, S.; Ruzzene, M. Analysis of In-plane Wave Propagation in Hexagonal and Re-entrant Lattices. J. Sound Vib. 2008, 312, 125–139. [Google Scholar] [CrossRef]
  22. Liebold-Ribeiro, Y.; Körner, C. Phononic Band Gaps in Periodic Cellular Materials. Adv. Eng. Mater. 2014, 16, 328–334. [Google Scholar] [CrossRef]
  23. Meng, J.; Deng, Z.; Zhang, K.; Xu, X.; Wen, F. Band Gap Analysis of Star-shaped Honeycombs with Varied Poisson’s Ratio. Smart Mater. Struct. 2015, 24, 095011. [Google Scholar] [CrossRef]
  24. Wojciechowski, K.W.; Scarpa, F.; Grima, J.N. Auxetics and Other Systems of Anomalous Characteristics. Phys. Status Solidi B-Basic Solid Status Phys. 2017, 254, 1770266. [Google Scholar] [CrossRef] [Green Version]
  25. Qiu, C.; Liu, Z.; Shi, J.; Chan, C. Directional Acoustic Source Based on the Resonant Cavity of Two-dimensional Phononic Crystals. Appl. Phys. Lett. 2005, 86, 224105. [Google Scholar] [CrossRef]
  26. Liu, W.; Su, X. Collimation and Enhancement of Elastic Transverse Waves in Two-dimensional Solid Phononic Crystals. Phys. Lett. A 2010, 374, 29682971. [Google Scholar] [CrossRef]
  27. Cicek, A.; Kaya, O.A.; Ulug, B. Wide-band All-angle Acoustic Self-collimation by Rectangular Sonic Crystals with Elliptical Bases. J. Phys. D Appl. Phys. 2011, 44, 205104. [Google Scholar] [CrossRef]
  28. Cicek, A.; Kaya, O.A.; Ulug, B. Impacts of Uniaxial Elongation on the Bandstructures of Two-dimensional Sonic Crystals and Associated Applications. Appl. Acoust. 2012, 73, 28–36. [Google Scholar] [CrossRef]
  29. Cebrecos, A.; Romero-Garcia, V.; Pico, R.; Perez-Arjona, I.; Espinosa, V.; Sanchez-Morcillo, V.; Staliunas, K. Formation of Collimated Sound Beams by Three-dimensional Sonic Crystals. J. Appl. Phys. 2012, 111, 104910. [Google Scholar] [CrossRef]
  30. Morvan, B.; Tinel, A.; Vasseur, J.; Sainidou, R.; Rembert, P.; Hladky-Hennion, A.-C.; Swinteck, N.; Deymier, P. Ultra-directional Source of Longitudinal Acoustic Waves Based on a Two-dimensional Solid/solid Phononic Crystal. J. Appl. Phys. 2014, 116, 214901. [Google Scholar] [CrossRef]
  31. Soliveres, E.; Espinosa, V.; Pérez-Arjona, I.; Sánchez-Morcillo, V.J.; Staliunas, K. Self Collimation of Ultrasound in a Three-dimensional Sonic Crystal. Appl. Phys. Lett. 2009, 94, 164101. [Google Scholar] [CrossRef]
  32. Liu, G.-T.; Tsai, C.-N.; Chang, I.-L.; Chen, L.-W. Elastic Dispersion Analysis of Two Dimensional Phononic Crystals Composed of Auxetic Materials. In Proceedings of the ASME 2015 International Mechanical Engineering Congress and Exposition, Houston, TX, USA, 13–19 November 2015; American Society of Mechanical Engineers: New York, NY, USA; p. V013T16A021. [Google Scholar]
  33. Tsai, C.-N.; Chen, L.-W. The manipulation of self-collimated beam in phononic crystals composed of orientated rectangular inclusions. Appl. Phys. A 2016, 122, 659. [Google Scholar] [CrossRef]
  34. Graff, K.F. Wave Motion in Elastic Solids; Courier Corporation: Chelmsford, MA, USA, 2012. [Google Scholar]
  35. Tang, H.-W.; Chou, W.-D.; Chen, L.-W. Wave propagation in the polymer-filled star-shaped honeycomb periodic structure. Appl. Phys. A 2017, 123, 523. [Google Scholar] [CrossRef]
Figure 1. The illustration of phononic crystal. Two cross sections of the inclusions are considered. They are (a) circular cross section and (b) rectangular cross section.
Figure 1. The illustration of phononic crystal. Two cross sections of the inclusions are considered. They are (a) circular cross section and (b) rectangular cross section.
Crystals 09 00121 g001aCrystals 09 00121 g001b
Figure 2. The band structure of the elastic wave in the phononic crystal with auxetic matrix (a) ν = 0.5 (b) ν = 0.7 (c) ν = 0.9 .
Figure 2. The band structure of the elastic wave in the phononic crystal with auxetic matrix (a) ν = 0.5 (b) ν = 0.7 (c) ν = 0.9 .
Crystals 09 00121 g002aCrystals 09 00121 g002b
Figure 3. The scheme of the auxetic phononic crystals. (a) The unit cell of the phononic crystals composed of star-shaped honeycomb structure and circular shaped inclusion; (b) the unit cell structure of the star-shaped honeycomb. The parameter of star-shaped honeycomb is l = 0.05   m , e = 0.05   m , t l = 0.005   m and t a = 0.005   m . The filling ratio of cylinder are r = 0.3 a .
Figure 3. The scheme of the auxetic phononic crystals. (a) The unit cell of the phononic crystals composed of star-shaped honeycomb structure and circular shaped inclusion; (b) the unit cell structure of the star-shaped honeycomb. The parameter of star-shaped honeycomb is l = 0.05   m , e = 0.05   m , t l = 0.005   m and t a = 0.005   m . The filling ratio of cylinder are r = 0.3 a .
Crystals 09 00121 g003
Figure 4. The band structure of phononic crystals composed of steel cylinders and star-shaped honeycombs. (a) θ = 55 ° , ν = 0.53 (b) θ = 60 ° , ν = 0.35 (c) θ = 65 ° , ν = 0.15 .
Figure 4. The band structure of phononic crystals composed of steel cylinders and star-shaped honeycombs. (a) θ = 55 ° , ν = 0.53 (b) θ = 60 ° , ν = 0.35 (c) θ = 65 ° , ν = 0.15 .
Crystals 09 00121 g004
Figure 5. The equi-frequency contour of the third band in Figure 4a. The star angle and the effective Poisson’s ratio are θ = 55 ° and ν = 0.53 , respectively.
Figure 5. The equi-frequency contour of the third band in Figure 4a. The star angle and the effective Poisson’s ratio are θ = 55 ° and ν = 0.53 , respectively.
Crystals 09 00121 g005
Figure 6. The simulation of wave propagation for negative refraction. A plane elastic wave of 500 Hz is incident from the left side of phononic crystal. The width of the wave and the incident angle are 4 a and 30 ° , respectively. The refraction angle is about 45 ° and the refraction ratio is about 0.707 . The material on both sides of the structure is steel.
Figure 6. The simulation of wave propagation for negative refraction. A plane elastic wave of 500 Hz is incident from the left side of phononic crystal. The width of the wave and the incident angle are 4 a and 30 ° , respectively. The refraction angle is about 45 ° and the refraction ratio is about 0.707 . The material on both sides of the structure is steel.
Crystals 09 00121 g006
Figure 7. The numerical results of the wave propagation of the phononic crystals composed of inclusions with rectangular cross section and star-shaped honeycomb matrix. (a) the band structures; (b) the equi-frequency contour of the 8th band. The star angle is θ = 55 ° , Poisson’s ratio of the matrix is ν = 0.53 . The geometric parameters of the inclusions are w = 14 a / 15 , h = a / 5 and ϕ = 0 ° .
Figure 7. The numerical results of the wave propagation of the phononic crystals composed of inclusions with rectangular cross section and star-shaped honeycomb matrix. (a) the band structures; (b) the equi-frequency contour of the 8th band. The star angle is θ = 55 ° , Poisson’s ratio of the matrix is ν = 0.53 . The geometric parameters of the inclusions are w = 14 a / 15 , h = a / 5 and ϕ = 0 ° .
Crystals 09 00121 g007
Figure 8. The displacement field simulations of the wave propagation in phononic crystals composed of star-shaped honeycomb and rectangular inclusions for (a) ϕ = 0 ° , (b) ϕ = 25 ° and (c) ϕ = 60 ° . The star angle of the matrix is θ = 55 ° and the dimensions of the rectangular inclusions are w = 14 / 15 a , h = a / 5 .
Figure 8. The displacement field simulations of the wave propagation in phononic crystals composed of star-shaped honeycomb and rectangular inclusions for (a) ϕ = 0 ° , (b) ϕ = 25 ° and (c) ϕ = 60 ° . The star angle of the matrix is θ = 55 ° and the dimensions of the rectangular inclusions are w = 14 / 15 a , h = a / 5 .
Crystals 09 00121 g008aCrystals 09 00121 g008b
Figure 9. The displacement field simulation of the wave propagation in the phononic crystal composed of star-shaped honeycomb matrix and rectangular inclusions with gradually varied orientations in the x -direction. The orientation angle is ϕ = 0 ° at the left boundary, x = 0 , and linearly increases to ϕ = 50 ° at the right boundary, x = 100 a . The orientation angles ϕ are linearly varied by ϕ ( x ) = ( x / 100 a ) 50 ° .
Figure 9. The displacement field simulation of the wave propagation in the phononic crystal composed of star-shaped honeycomb matrix and rectangular inclusions with gradually varied orientations in the x -direction. The orientation angle is ϕ = 0 ° at the left boundary, x = 0 , and linearly increases to ϕ = 50 ° at the right boundary, x = 100 a . The orientation angles ϕ are linearly varied by ϕ ( x ) = ( x / 100 a ) 50 ° .
Crystals 09 00121 g009
Table 1. The effective material properties of steel star-shaped honeycomb, calculated by COMSOL FEM software [35].
Table 1. The effective material properties of steel star-shaped honeycomb, calculated by COMSOL FEM software [35].
STAR   ANGLE   θ 55 ° 60 ° 65 ° 70 °
Young’s module E ( G p a ) 0.51240.31230.23050.1806
Poisson’s ratio ν −0.53−0.35−0.150.05
Mass density ρ ( k g / m 3 ) 921921921921

Share and Cite

MDPI and ACS Style

Chang, S.-Y.; Chen, C.-D.; Yeh, J.-Y.; Chen, L.-W. Elastic Wave Propagation of Two-Dimensional Metamaterials Composed of Auxetic Star-Shaped Honeycomb Structures. Crystals 2019, 9, 121. https://doi.org/10.3390/cryst9030121

AMA Style

Chang S-Y, Chen C-D, Yeh J-Y, Chen L-W. Elastic Wave Propagation of Two-Dimensional Metamaterials Composed of Auxetic Star-Shaped Honeycomb Structures. Crystals. 2019; 9(3):121. https://doi.org/10.3390/cryst9030121

Chicago/Turabian Style

Chang, Shu-Yeh, Chung-De Chen, Jia-Yi Yeh, and Lien-Wen Chen. 2019. "Elastic Wave Propagation of Two-Dimensional Metamaterials Composed of Auxetic Star-Shaped Honeycomb Structures" Crystals 9, no. 3: 121. https://doi.org/10.3390/cryst9030121

APA Style

Chang, S. -Y., Chen, C. -D., Yeh, J. -Y., & Chen, L. -W. (2019). Elastic Wave Propagation of Two-Dimensional Metamaterials Composed of Auxetic Star-Shaped Honeycomb Structures. Crystals, 9(3), 121. https://doi.org/10.3390/cryst9030121

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