Microsoft Word - numero_29_art_1 A. Bacigalupo et alii, Frattura ed Integrità Strutturale, 29 (2014) 1-8; DOI: 10.3221/IGF-ESIS.29.01 1 Focussed on: Computational Mechanics and Mechanics of Materials in Italy A micropolar model for the analysis of dispersive waves in chiral mass-in-mass lattices A. Bacigalupo University of Trento andrea.bacigalupo@unitn.it L. Gambarotta University of Genoa luigi.gambarotta@unige.it ABSTRACT. The possibility of obtaining band gap structures in chiral auxetic lattices is here considered and applied to the case of inertial locally resonant structures. These periodic materials are modelled as beam-lattices made up of a periodic array of rigid rings, each one connected to the others through elastic slender ligaments. To obtain low-frequency stop bands, elastic circular resonating inclusions made up of masses located inside the rings and connected to them through an elastic surrounding interface are considered and modeled. The equations of motion are obtained for an equivalent homogenized micropolar continuum and the overall elastic moduli and the inertia terms are given for both the hexachiral and the tetrachiral lattice. The constitutive equation of the beam lattice given by the Authors [15] are then applied and a system of six equations of motion is obtained. The propagation of plane waves travelling along the direction of the lines connecting the ring centres of the lattice is analysed and the secular equation is derived, from which the dispersive functions may be obtained. KEYWORDS. Auxetic materials; Chirality; Cellular materials; Mass-in-mass dynamic systems; Dispersive waves. INTRODUCTION n recent years a considerable interest in acoustic metamaterials was witnessed by several research (see for reference [1]), most of them focused to the design of artificial materials having periodic microstructure conceived to get complete sound attenuation for a certain frequency range, namely acoustic wave spectral gap. Sonic crystals with spectral gaps [2] have been developed on the realization that composites with locally resonant structural units may exhibit effective negative elastic constants at certain frequency ranges, as shown by Liu et al. [3] and by Huang et al. [4,5]. Recently, Bigoni et al. [6] proposed a periodic metamaterial with internal locally resonant structures that supports tunable low- frequency stop bands. This effect is associated with localized rotational modes obtained from a chiral microstructure of the periodic cell. Tee et al. [7] have recently obtained a band gap structure in periodic tetrachiral materials. Spadoni et al. [8] have obtained analogous results with reference to hexachiral lattices. These kind of materials are auxetics [9] and their structure is I A. Bacigalupo et alii, Frattura ed Integrità Strutturale, 29 (2014) 1-8; DOI: 10.3221/IGF-ESIS.29.01 2 characterized by a periodic array of the rings connected by four or six elastic ligaments. Lakes [10], with reference to the planar isotropic case of hexachiral lattice, proposed this microstructure geometry firstly. Later Prall and Lakes [11] showed for this material a Poisson’s ratio of -1 under the hypothesis of ignoring the axial strain of the ligaments. Further studies on the homogenization of hexachiral auxetic materials have been carried out by Spadoni and Ruzzene [12] and by Liu et al. [13]. Liu et al. [14] proposed a hexachiral metacomposite by integrating a two-dimensional hexachiral lattice with elastic resonating inclusions to obtain low-frequency band gap. This metacomposite has been analysed through a numerical model where the ligaments are modelled as multi-beam elements and the inclusions as a two dimensional FEM model. In this paper, hexa- and tetra-chiral beam lattices are considered having local resonators at the nodes of the periodic array. The model is developed in closed-form and is based on a micropolar homogenization of the lattice. This approach partially relies on the results by Bacigalupo and Gambarotta [15], which developed and compared the results from the micropolar and a second displacement gradient homogenization for both the hexachiral and the tetrachiral periodic cells in order to evaluate the validity limits of the beam lattice model. The possibility of obtaining band gap structures in chiral auxetic lattices is here considered and applied to the case of inertial locally resonant structures. These structures are obtained following Huang et al. [4] and Liu et al. [14], by the insertion of a circular mass connected through an elastic surrounding interface to each ring of the microstructure. The equations of motion are given within a micropolar continuum model and the overall elastic moduli and the inertia terms are obtained for both the hexachiral and the tetrachiral lattice. The constitutive equation of the beam lattice given in [15] are then applied and a system of six equations of motion is obtained. The propagation of plane waves travelling along the direction of the lines connecting the ring centres of the lattice is analysed and the secular equation is derived, from which the dispersive functions may be obtained. (a) (b) Figure 1: (a) Hexachiral lattice; (b) tetrachiral lattice. CHIRAL MASS-IN-MASS PERIODIC MATERIAL: MICROPOLAR HOMOGENIZATION he periodic materials shown in Figure 1 are considered as beam-lattices made up of a periodic array of rigid rings, each one connected to the others through n elastic slender ligaments rigidly connected to the rings. In Figure 2.a the periodic hexachiral cell (n=6) is shown, while in Figure 1.b is shown the tetrachiral cell (n=4). Each ligament is tangent to the joined rigid rings and has length l measured between the connection points, with section width t, thickness d and Young’s modulus sE . Both the ligaments and the rings have mass density s , so that the mass and the rotation inertia of the rings are 1M 2 s rt  and 2 1 1J M r , respectively. The two dimensional composite materials are auxetic. The hexachiral lattice is isotropic and its Poisson’s ratio becomes negative when increasing the chirality angle  . On the other hand, the tetrachiral material is strongly anisotropic and presents a directional dependency of the chirality on the direction of the applied stress, as shown in [15]. To obtain low-frequency stop bands, a metamaterial inclusion consisting of a softly heavy disk is located inside the ring as shown in Figure 3.a. This inclusion plays the role of low-frequency resonator with mass and rotation inertia denoted with 2M and 2J , respectively. The motion of the rigid ring of the beam lattice is denoted by the displacement vector u and the rotation  , respectively (see Figure 3.b), while the motion of the mass of the resonator is denoted by the displacement vector v and the rotation  (see Figure 3.c). T A. Bacigalupo et alii, Frattura ed Integrità Strutturale, 29 (2014) 1-8; DOI: 10.3221/IGF-ESIS.29.01 3 (a) (b) Figure 2: Periodic cell of the (a) hexachiral lattice; (b) tetrachiral lattice. (a) (b) (c) Figure 3: (a) Internal mass with external elastic thick layer; (b) rigid ring and related dofs; (c) internal mass and related dofs. A soft elastic interface connects the internal mass to the rigid ring. To get a simplified formulation, the constitutive equation of the interface is assumed in the form    , c = dk k      f v u (1) f being the force exerted by the rigid ring on the internal mass and c the corresponding couple (see Figure 4). The parameters dk and k are the isotropic translational stiffness and the rotational stiffness, respectively. Figure 4: Contact force and couple between the rigid ring and the internal mass. The lattice model is here approximated as a micropolar continuum model resulting from a homogenization process based on the macro-homogeneity criterion, involving both the total potential energy and the kinetic energy through the Hamilton principle (see [13, 16] for reference). The equations of motion for the homogenized beam-lattice are A. Bacigalupo et alii, Frattura ed Integrità Strutturale, 29 (2014) 1-8; DOI: 10.3221/IGF-ESIS.29.01 4             11,1 12,2 1 1 1 1 21,1 22,2 2 2 1 2 1,1 2,2 21 12 1 1 1 2 1 2 2 2 2 2 ˆ ˆ ˆ ˆ ˆ ˆ d d d d k v u u k v u u m m k I k u v v k u v v k I                                                  (2) where the stiffnesses ˆd d cellk k A and ˆ cellk k A  are introduced, cellA being the area of the periodic cell, together with the mass densities M11 cellA   and M 22 cellA   , and the micro-inertia terms J11 cell I A and J 2 2 cell I A , respectively. In Eq. (2) ij and im are the overall stress components, namely the asymmetric stress components and the micro- couples, respectively. In case of hexachiral lattice, one obtains 2 22 3 coscell lA   and the following inertial parameters are obtained  1 3 sin 2 24 s       2 21 3 tan sin 24 sI l       22 2 3 sin 24 i cell M A       2 2 222 3 sin tan 192 i cell J I l A       For the tetrachiral lattice one obtains 2 2coscell lA   and the resulting inertia parameters are  1 sin 2 24 s       , 2 2 1 sin 48 sI l      , 222 sin 48 i cell M A       , 2 2 222 sin tan 192 i cell J I l A       . To obtain the displacement formulation of the equations of motion, the compatibility equations involving the macrostrain components 11 1,1u  , 22 2,2u  , 12 1,2u   ,   21 2,1u and the curvatures 1 ,1   and 2 ,2  have to be considered together with the constitutive equation. MICROPOLAR CONSTITUTIVE EQUATION FOR HEXA- AND TETRA-CHIRAL LATTICES he constitutive equation of hexachiral honeycomb corresponds to that obtained in [13, 15] and is written as follows 11 11 22 22 12 12 21 21 1 1 2 2 2 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 A A A A A A k k A A k k m S m S                                                                             (3) in which the five elastic moduli T A. Bacigalupo et alii, Frattura ed Integrità Strutturale, 29 (2014) 1-8; DOI: 10.3221/IGF-ESIS.29.01 5            2 2 2 2 2 2 2 2 2 2 2 2 3 1 4 3 1 cos sin 4 3 sin cos 2 3 3sin 4 cos 12 3 1 sin cos 2 s s s s s E E k E S E a A E                                 (4) depend on the Young’s modulus, the slenderness ratio t l  and the angle  of inclination of the ligaments. The constitutive Eq. (3) show the coupling between the extensional strains 11 and 22 and the asymmetric strains 12 and 21 through the elastic constant A. The elastic moduli, with the exception of parameter  , depend on the parameter of chirality  , but only the constant A is an odd function of this parameter, i.e. it reverses its sign when the handedness of the material pattern is flipped over. In case of symmetric macro-strain fields, the fourth order elastic tensor for the hexagonal system corresponds to that of the transversely isotropic system whose elastic moduli in the plane of the lattice are:          2 4 2 2 3 2 hom 2 2 2 2 hom 2 4 2 2 2 hom 3 cos sin 2 3 1 cos sin 1 3 cos sin 3 = 1 4 s s E E E                           G (5) For the tetrachiral lattice, the constitutive equation is written [15]: 11 11 22 22 12 12 21 21 1 1 2 2 2 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 B B B B m S m S                                                               (6) where four elastic moduli are related to the lattice parameters as follows         2 2 2 2 2 2 2 2 2 2 2 cos sin 2 sin cos 1 sin cos 1 3sin 4 cos 12 s s s s E k E B E S E a                           (7) Similarly to the hexachiral honeycomb, a coupling is obtained between the extensional strains and the asymmetric strains through the elastic modulus B which is an odd function of the parameter of chirality  , while the other elastic moduli are even functions. In case of symmetric macro-strain fields, the resulting classical fourth order elasticity tensor has the elastic A. Bacigalupo et alii, Frattura ed Integrità Strutturale, 29 (2014) 1-8; DOI: 10.3221/IGF-ESIS.29.01 6 moduli of the tetragonal system. The elasticity tensor depends on the chirality parameter, but unlike the hexachiral honeycomb, some elasticities are odd functions of  [15]. PLANE WAVES PROPAGATION he free wave motion in the hexachiral material is written in terms of the components of the generalized displacement field    1 2 1 2 T u u v v  U x in the following form             1,11 1,12 1,22 2 ,11 2,12 2,22 ,1 ,2 1 1 1 1 1,11 1,12 1,22 2,11 2,12 2,22 ,1 ,2 2 2 1 2 ,11 ,22 1,1 2,2 1,2 2 ( 2 ) 2 2 2 2 ( 2 ) 2 2 2 2 2 2 d d u Au k u Au k u Au A k k v u u Au k u Au k u Au u k A k v u u S S Au Au ku ku                                                                ,1 1 1 1 2 1 2 2 2 2 2 4 ˆ , ˆ , ˆ . d d k k I k u v v k u v v k I                                  (8) If a harmonic plane wave propagating along axis 1x in an infinite planar micropolar medium is admitted, the generalized displacement field at a point is assumed in the following form  1ˆ expT i qx t   U U , where q and  denote the wave number and the circular frequency, respectively, and  1 2 1 2ˆ ˆˆ ˆ ˆ ˆ ˆT u u v v  U is the vector of the amplitudes  1i   . Substituting the assumed generalized displacement field in the equation of motion (8) one obtains the secular equation system for the equivalent continuum model     2 2 2 1 2 2 2 1 2 2 1 2 2 2 2 2 2 2 ˆ2 0 0 ˆ ˆ2 0 0 ˆ 4 ˆ2 2 0 0 ˆ ˆ ˆ0 0 0 0 ˆ ˆ0 0 0 0 ˆ ˆ0 0 0 0 d d d d d d d d q Aq iAq k k k q Aq ikq k k Sq k iAq ikq k k I k k k k k k I                                                                       1 2 1 2 ˆ ˆ ˆ ˆ ˆ ˆ u u v v                            0 (9) The solution of the eigenproblem (9), characterized by a hermitian matrix, provides dispersion functions  q defined in the domain ,q a a        , a being the characteristic size of the cell as shown in Figure 2. Since the model has six degrees of freedom, it follows that in the general case six plane waves may propagate, each characterized by a dispersion function. It follows that this model exhibits plane waves more complex than those obtained and discussed by Parfitt and Eringen [17] for isotropic achiral two dimensional micropolar solids and more recently by Khurana and Tomar [18] for 3D chiral T A. Bacigalupo et alii, Frattura ed Integrità Strutturale, 29 (2014) 1-8; DOI: 10.3221/IGF-ESIS.29.01 7 micropolar media. Since the matrix of the coefficient in the eigenproblem (9) is hermitian, the eigenvalues are real, but not necessarily positive. Therefore, for some values of the model parameters is expected a reduced number of dispersive functions. In case of a tetrachiral lattice with internal resonant masses the free wave equation of motion takes the form:           1,11 1,22 2,11 2,22 ,1 ,2 1 1 1 1 1,11 1,22 2,11 2,22 ,1 ,2 2 2 1 2 ,11 ,22 1,1 2 ,2 1,2 2,1 1 1 1 2 1 2 2 2 2 2 2 2 ˆ ˆ ˆ d d d d u ku Bu Bu B k k v u u Bu Bu ku u k B k v u u S S Bu Bu ku ku k k I k u v v k u v v                                                         2 k I                 (10) In this case the secular equation system is written as follows: 2 2 2 1 2 2 2 1 2 2 1 2 2 2 2 2 2 ˆ2 ˆ 0 0 ˆ ˆ0 0 2 ˆ0 0 ˆ ˆ ˆ0 0 0 0 ˆ ˆ0 0 0 0 ˆ ˆ0 0 0 0 d d d d d d d d q k Bq iBq k kq k Bq ikq k Sq k iBq ikq k k I k k k k k k I                                                                         1 2 1 2 ˆ ˆ ˆ ˆ ˆ ˆ u u v v                         0 (11) Again the solution of the eigenproblem (11) with hermitian matrix provides dispersion functions  q , defined in the domain ,q a a       . The analytical formulation here derived, which is based on a micropolar continuum model enriched with additional degrees of freedom of the resonating masses located inside the rings, will be developed in future research to analyze the wave propagation and to appreciate the influence of the geometrical and mechanical parameters of the microstructure and of the resonant devices on the frequency spectra and on the conditions of wave attenuation. The reliability of the results has to be evaluated from the comparison with the rigorous solutions obtained from a Floquet-Bloch analysis of the generalized beam lattice proposed in this paper. Finally, it should be noted that the validity of this model relies on the rigidity assumption of the rings [15], a condition that may be easily implemented in the physical model. ACKNOWLEDGEMENTS he Authors acknowledge financial support of the (MURST) Italian Department for University and Scientific and Technological Research in the framework of the research MIUR Prin09 project XWLFKW, Multi-scale modelling of materials and structures, coordinated by prof. A. Corigliano. Andrea Bacigalupo gratefully thanks financial T A. Bacigalupo et alii, Frattura ed Integrità Strutturale, 29 (2014) 1-8; DOI: 10.3221/IGF-ESIS.29.01 8 support of the Italian Ministry of Education, University and Research in the framework of the FIRB project 2010, “Structural mechanics models for renewable energy applications". REFERENCES [1] Lu, M.H., Feng, L., Chen, Y.-F., Phononic crystals and acoustic metamaterials, MaterialsToday, 12 (2009) 34-42. [2] Pennec, Y., Vasseur, J.O., Djafari-Rouhani, B., Dobrzyński, L., Deymier, P.A., Two-dimensional phononic crystals: Examples and applications, Surface Science Reports, 65 (2010) 229–291. [3] Liu, Z., Zhang, X., Mao, Y., Zhu, Y., Yang, Z., Chang, C.T., Sheng, P., Locally Resonant Sonic Materials, Science, 289 (2000) 1734-1736. [4] Huang, H.H., Sun, C.T., Huang, G.L., On the negative effective mass density in acoustic metamaterials, Int. J. of Engineering Science, 47 (2009) 610-617. [5] Huang, H.H., Sun, C.T., Wave attenuation mechanism in an acoustic metamaterial with negative effective mass density, New Journal of Physiscs, 11 (2009) 013003. [6] Bigoni, D., Guenneau, S., Movchan, A. B., Brun M., Elastic metamaterials with inertial locally resonant structures: Application to lensing and localization, Physical Review, B 87 (2013) 174303. [7] Tee, K.F., Spadoni, A., Scarpa, F., Ruzzene, M., Wave propagation in auxetic tetrachiral honeycombs, J. of Vibration and Acoustics ASME, 132 (2010) 031007-1/8. [8] Spadoni, A., Ruzzene, M., Gonella, S., Scarpa, F., Phononic properties of hexagonal chiral lattices, Wave Motion, 46 (2009) 435-450. [9] Prawoto, Y., Seeeing auxetic materials from the mechanics point of view: A structural review on the negative Poisson's ratio, Computational Material Science, 58 (2012) 140-153. [10] Lakes, R. S.,Foam structures with a negative Poisson’s ratio, Science, 235 (1987) 1038-1040. [11] Prall, D., Lakes, R.S., Properties of chiral honeycomb with a Poisson ratio of -1, Int. J. Mechanical Sciences, 39 (1997) 305-314. [12] Spadoni, A., Ruzzene, M., Elasto-static micropolar behavior of a chiral auxetic lattice, J. Mechanics and Physics of Solids, 60 (2012) 156-171. [13] Liu, X.N., Huang, G.L., Hu, G.K., Chiral effect in plane isotropic micropolar elasticity and it’s application to chiral lattices, J. Mechanics and Physics of Solids, 60 (2012) 1907-1921. [14] Liu, X.N., Hu, G.K., Sun, C.T., Huang, G.L., Wave propagation characterization and design of two-dimensional elastic chiral metacomposite, J. of Sound and Vibration, 330 (2011) 2536-2553. [15] Bacigalupo, A., Gambarotta, L., Homogenization of periodic hexa- and tetrachiral cellular solids, Composite Structures, DOI 10.1016/j.compstruct.2014.05.033, 2014.to appear. [16] Stefanou, I., Sulem, J., Vardoulakis, I., Three-dimensional Cosserat homogenization of masonry structures: elasticity, Acta Geotechnica, 3 (2008) 71–83. [17] Parfitt, V.R., Eringen, A.C., Reflection of Plane Waves from the Flat Boundary of a Micropolar Elastic Half-Space, J. Acoustical Society of America, 45 (1969) 1258-1272. [18] Khurana, A., Tomar, S.K., Longitudinal wave response of a chiral slab interposed between micropolar solid half- spaces, Int. J. of Solids and Structures, 46 (2009) 135–150.