This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
|
|
- Augustus Cooper
- 5 years ago
- Views:
Transcription
1 This article appeared in a journal pulished y Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution and sharing with colleagues. Other uses, including reproduction and distriution, or selling or licensing copies, or posting to personal, institutional or third party wesites are prohiited. In most cases authors are permitted to post their version of the article (e.g. in Word or Tex form) to their personal wesite or institutional repository. Authors requiring further information regarding Elsevier s archiving and manuscript policies are encouraged to visit:
2 Wave Motion 49 (01) Contents lists availale at SciVerse ScienceDirect Wave Motion journal homepage: Multi-displacement microstructure continuum modeling of anisotropic elastic metamaterials A.P. Liu a,, R. Zhu, X.N. Liu a,, G.K. Hu a,, G.L. Huang, a School of Aerospace Engineering, Beijing Institute of Technology, Beijing, , China Department of Systems Engineering, University of Arkansas at Little Rock, AR, 704, USA a r t i c l e i n f o a s t r a c t Article history: Received 6 July 011 Received in revised form 0 Octoer 011 Accepted 6 Decemer 011 Availale online 8 January 01 Keywords: Metamaterial Multi-displacement Continuum Anisotropic mass density Wave propagation An elastic metamaterial made of lead cylinders coated with elliptical ruers in an epoxy matrix is considered, and its anisotropic effective dynamic mass density tensor is numerically determined and demonstrated. To capture oth dipolar resonant motion and microstructure deformation in the composite, a new multi-displacement microstructure continuum model is proposed. In the formulation, additional displacement and kinematic variales are introduced to descrie gloal and local deformations, respectively. The macroscopic governing equations of the two-dimensional anisotropic elastic metamaterial are explicitly derived through a simplified procedure. To verify the current model, wave dispersion curves from the current model are compared with those from the finite element simulation for oth longitudinal and transverse waves. Very good agreement is oserved in oth the acoustic and optic wave modes. This work could provide a enchmark of continuum modeling of elastic metamaterials with nonelementary microstructures. 011 Elsevier B.V. All rights reserved. 1. Introduction Metamaterials are a new class of ordered composites that exhiit exceptional electromagnetic, optic and mechanical properties not readily oserved in nature [1,]. Although electromagnetic (EM) and acoustic metamaterials have een studied intensively, elastic metamaterials have received much less attention, despite the fact that elastic metamaterials offer richer ehaviors as they enale oth longitudinal and transverse waves to propagate. Various novel concepts and engineering applications of elastic metamaterials have een successfully demonstrated such as mechanical filters, sound and viration isolators, elastic waveguides and energy harvesting [ 8]. To understand gloal wave mechanism in elastic metamaterials, the most important and efficient approach is homogenization. For isotropic elastic metamaterials, conventional homogenization will result in three independent effective parameters: effective mass density, effective ulk modulus and effective shear modulus. Until now, several homogenization methods have een developed to model this kind of materials as effective homogeneous elastic media. Those methods include the plane wave expansion (PWE) technique, the coherent potential approximation (CPA), the micromechanics analytical model or multiple scattering theory (MST) ased on the long-wavelength limit [9 16]. By requiring only the wavelength in the matrix to e much larger than the size of the microstructures, the extension of the CPA or effective medium theory (EMT) was further developed to predict the effective properties of the elastic metamaterials [17]. However, these analytical methods are mainly ased on the analytical scattering solution of microstructures with simple geometries (sphere or cylinder). To achieve the desirale effective properties, elastic metamaterials with complex microstructures or Corresponding authors. addresses: hugeng@it.edu.cn (G.K. Hu), glhuang@ualr.edu (G.L. Huang) /$ see front matter 011 Elsevier B.V. All rights reserved. doi: /j.wavemoti
3 41 A.P. Liu et al. / Wave Motion 49 (01) microstructure arrangements must e designed. For nonelementary microstructure geometries, an exact analytical model is not practicale and computational techniques should e suggested such as the multiple multipole (MMP) method [18], lumped mass method [19], and finite difference method [0]. For the acoustic metamaterial, the anisotropic effective mass density can e achieved y changing position or distriution array of the microstructure in the unit cell [1,]. However, for the elastic metamaterial, the microstructure geometry should e properly designed to otain the anisotropic effective mass density [ 6]. A D elastic resonator of aritrary geometry was systematically studied through a finite element modal analysis [7]. Gu et al. [8] investigated local resonance modes of elliptic cylinders coated with silicon ruer in a rigid matrix to otain the anisotropic effective mass density. However, most of these methods mainly focus on an understanding of the local resonant mechanism. Practically, one may want to predict dynamic ehavior of the elastic metamaterial y using an effective continuum model, wherein effects of the microstructure are included in effective properties and/or in macrodynamics equations. The otained macro-dynamics equations can e applied to prolems of time-dependent viration and transit wave propagation in the finite and infinite structures, which is important for potential engineering applications of metamaterial devices. However, continuum homogenization of elastic metamaterials with complex microstructures is a challenging task. Until now, few efforts have een made to study the effective properties eyond the quasi-static limit or in higher frequency regions, although some interesting wave phenomena occur only when the wavelength is comparale to the lattice scale. To descrie wave propagation prolems with wavelengths in an order of dimensions of microstructures, Wozniak [9] proposed a non-asymptotic approach of the macro-dynamics modeling of composite. Milton and Willis [0] and Milton [1] presented a rigorous theoretical foundation for the dynamic effective parameters of the elastic metamaterial eyond the quasi-static limit. The microscale effects due to non-local ehavior were explained y the modification of linear continuum elastodynamic equations. Similar efforts have een made through work on metamaterials with a macroscopic higher order gradient or non-local elastic response []. An alternative approach is to employ additional kinematic variales to descrie the nonhomogeneous local deformation in the microstructure of the solid. This approach leads to Cosserat continuum models [,4] or micropolar models [5,6] or strain gradient theory [7] or the microstructure continuum theory [8,9]. Microstructure continuum theory has recently een adopted to descrie gloal dynamic ehavior of isotropic elastic metamaterials with discrete microstructures (a mass-spring system) [40]. Accuracy of the model was verified y comparison with the results otained from the finite element method. However, the microstructure continuum model or high-order continuum model has difficulty in capturing the dipolar motion, in which the inner inclusion moves out of phase with respect to motion of the matrix in elastic metamaterials. Additional displacement variales are needed to capture the dipolar motion. In this paper, the anisotropic dynamic mass density tensor of an elastic metamaterial made of lead cylinders coated with elliptical ruers in an epoxy matrix is first numerically determined. To analytically otain gloal governing equations of motion in the elastic metamaterial, the complex continuous microstructure is simplified to a discrete mass-spring system ased on strain energy equivalence. Based on the simplified system, a new multi-displacement microstructure continuum theory is developed to capture oth dipolar local resonant motion and microstructure deformation y introducing oth multi-displacement variales and micro-deformation variales, which are not considered in the conventional continuum model. The developed model is verified and applied to investigate gloal wave propagation in D anisotropic elastic metamaterials under the plane stress assumption.. Anisotropic elastic metamaterials In this study, an anisotropic elastic metamaterial made of heavy cylinder cores coated with elliptical soft layer and emedded in a matrix is considered, as shown in Fig. 1(a), and its representative volume element (RVE) is identified in Fig. 1(). The microstructure is distriuted in a rectangular lattice array. In the figure, the hard inclusion core is laeled as medium 1, the soft coating medium is laeled as medium, and the matrix is laeled as medium. The isotropic material constants of the inclusion core, the coating layer and the matrix are ρ 1, E 1, ν 1 ; ρ, E, ν and ρ, E, ν, respectively. The radius of the core is a. The lattice constants along the X and X directions are denoted as d and d. The semimajor and semiminor axes of the ellipse are denoted as 1 and, respectively. In the unit cell (k, l), the position of the center point is given in the gloal coordinate X = X l and X = X k. For convenience, we define a local polar coordinate system (r, θ) as well as a local Cartesian coordinates (x, x ) in the unit cell with x = r cos θ and x = r sin θ, as shown in Fig. 1(). Strong anisotropic properties of the metamaterial along the X and X directions can e achieved y adjusting dimensions of the major and minor axes of the coating ellipse. Due to geometric complexity of the microstructure, analytical-ased methods cannot e applied directly for the determination of the effective dynamic mass density [8]. In this study, a numerical-ased method will e adopted [8]. Consider the unit cell in Fig. 1(), time-harmonic displacements u α = û α e iωt with α =, are prescried on the exterior oundary, and a harmonic analysis of the unit cell is performed y using the finite element method (FEM), from which resultant forces F α = ˆFα e iωt on the exterior oundary can e otained. The effective mass density tensor can then e numerically defined and determined y ˆF ρ eff = ω ρ eff A û c (1) ˆF ρ eff ρ eff û
4 A.P. Liu et al. / Wave Motion 49 (01) a Fig. 1. (a) An anisotropic elastic metamaterial made of cylinders coated with elliptical soft layer in a matrix; () The (k, l) element (RVE). Fig.. Coordinate systems of the unit cell with an aritrary direction. Tale 1 The microstructure geometry parameters. Lattice parameters Elliptical coat (mm) Circle coat (mm) d = d 0 0 a where A c = d d is the area of the unit cell. To demonstrate that the mass density follows the coordinate transformation law, an effective mass density tensor in an aritrary coordinate system (x, y), as shown in Fig., is similarly determined. In the figure, (x, x ) is a rectangular Cartesian coordinate system in principal directions and (x, y) is an aritrary Cartesian coordinate system. It should e mentioned that the unit cell in Fig. is specifically suggested for simplification of numerical calculation, in which the outer matrix seems to follow the unit cell transformation. This is ecause the outer host matrix is mainly functioned as an effective mass in the calculation of the effective mass density of the unit cell. Therefore, determination of the effective mass density is independent on the outer matrix geometry. Based on the numerical results of the unit cell, it can e proved that the effective mass density follows the coordinate transformation law as ρ eff xx ρ eff yx ρ eff xy ρ eff yy = C ρ eff S ρ eff CS(ρ eff ρeff ) CS(ρ eff ρeff ) C ρ eff S ρ eff where C = cos δ, S = sin δ, and δ is the angle of the transformation. Therefore, the effective mass density admits the second-order tensorial property. Fig. shows the effective mass density of the elastic metamaterial along a direction δ = 0 y using the FEM model, where ρ ave = (ρ 1 A 1 ρ A ρ A ) /A c is the average static mass density for the composite. The microstructure geometry parameters and the constituent material parameters are given in Tales 1 and, respectively. From the figure, we can clearly see that the effective mass density is frequency-dependent, anisotropic (ρ eff xx ρ eff yy, ρ eff xy 0) and ecomes negative around the resonant frequency. For the metamaterial with the current microstructure, two different resonant frequencies along x and x directions can e found as f 1 = 86 Hz and f = 950 Hz, respectively, which indicate anisotropy of the ()
5 414 A.P. Liu et al. / Wave Motion 49 (01) Tale The constituent material parameters. The parameters Core: lead Coating: ruer Matrix: epoxy Mass density ρ 1 : kg/m ρ : kg/m ρ : kg/m Young s modulus E 1 : N/m E : N/m E : N/m Poisson s ratio ν 1 : 0.45 ν : ν : 0.8 Area A 1 : πa A : π( 1 a ) A : d d π 1 Fig.. The anisotropic effective mass density of the elastic metamaterial predicted y the FEM model along a direction δ = 0. Fig. 4. The effective mass density with respect to the angle δ. metamaterial. The significant difference etween ρ eff xx and ρ eff yy can e oserved when the frequency is close to the resonant frequency. It is also interesting to note that the mass density ecome isotropic (ρ eff xx = ρ eff yy, ρ eff xy = 0) when the frequency is much less or much larger than the resonant frequency. Effective mass density for a D ternary phononic crystal with two resonators was also discussed y Gu et al. [8]. Fig. 4 shows the effective mass densities in functions of the angle δ at frequency f = 950 Hz. The strong anisotropy of the effective mass density can e clearly found (ρ eff xx ρ eff yy, ρ eff xy 0) with the angle change of the incident wave. It should e emphasized that the anisotropy of the effective mass density can e tuned y adjusting the geometry of the microstructure.
6 A.P. Liu et al. / Wave Motion 49 (01) Fig. 5. The equivalent RVE for the (k, l) element of the elastic metamaterial.. Multi-displacement microstructure continuum modeling of the elastic metamaterial In the previous section, the anisotropic effective mass density tensor of the elastic metamaterial and its mechanism are calculated and determined through the numerical analysis, from which the anisotropic negative mass density is illustrated. However, in order to understand gloal dynamic ehavior of the elastic metamaterial, an analytical-ased continuum model is needed to properly descrie the dispersive ehavior as well as and-gap ehavior. The otained macrodynamics equations can e applied directly to dynamic prolems of the elastic metamaterial with the finite and infinite structures, which is important for its potential engineering applications. In this section, a multi-displacement microstructure continuum model will e developed analytically for the first time. First, the analytical effective mass density tensor of the elastic metamaterial will e determined through a simplified model. Then additional displacement and kinematic variales will e introduced to capture the local dipolar motion and microstructure deformation..1. A simplified model of the elastic metamaterial For the elastic metamaterial with the complicated microstructure geometry, analytical-ased homogenization approaches have difficulties in finding the exact local scattering wave field. In this study, the RVE with the continuous microstructure is first simplified to the RVE with a discrete mass-spring microstructure with spring constants K and K, as shown in Fig. 5 [0], ased on the strain energy equivalence etween the two systems. The inclusion core is reduced as a rigid mass m 1 = ρ 1 πa and the coating material is replaced y the springs with the spring constants K and K ecause of the large stiffness mismatch etween the core and the coating material. The spring constants K and K can e numerically determined as K α = F α /γ α with α =,, where F α is the restoring force on the outer fixed oundary of the coating layer and u α = γ α is the applied displacement along x and x directions of the coating layer. For the general geometry of the coating layer, the prolem can e numerically solved y using the finite element method (FEM). Based on the analytical Lorentz model [0], four components of the effective density tensor of the composite system in Fig. 5 along principal direction can e otained analytically as ρ eff ρ eff ρ eff ρ eff = 1 A c m m 1K 0 K m 1 ω 0 m m 1K K m 1 ω () where A c = d d is the area of the RVE, and m = ρ A is the mass of the matrix in a RVE. The density tensor in an aritrary coordinate system (x, y) can e easily determined according to the second-order tensorial property. Fig. 6 shows a comparison of the effective mass density of the anisotropic elastic metamaterial with the microstructure along the direction δ = 0 from the numerical model and the current analytical Lorentz model. For the elliptical coat geometric parameters in Tale 1, the spring constants can e otained as K = kg/m and K = kg/m. Very good agreement etween these two models can e clearly oserved, which sustantiates the fact that the simplified discrete mass-spring model can e used as the dynamic equivalent system of the anisotropic elastic metamaterial in Fig. 1. A small difference around the resonant frequencies etween the two models is ecause the discrete Lorentz model ignores the mass of the coating layer and assumes the matrix to e rigid... A multi-displacement microstructure continuum model In this susection, a new homogeneous high-order continuum theory will e proposed to homogenize the simplified elastic anisotropic metamaterial, as shown in Fig. 7. For the elastic metamaterial, the inner mass will move out-of-phase
7 416 A.P. Liu et al. / Wave Motion 49 (01) a c Fig. 6. Comparison of the effective mass density predicted y the FEM model and the analytical Lorentz model. Fig. 7. Heterogeneous medium replaced y a homogeneous medium. with respect to the motion of the composite when the frequency is close to the resonant frequency of the inner mass, which is called the dipolar motion. The schematic illustration of the dipolar motion in the elastic metamaterial near the resonant frequency is shown in Fig. 8. The conventional continuum model cannot capture the dipolar motion, therefore, a continuum model with additional multi-displacement variales should e introduced. Specifically, attention is also paid to the prolem with wavelengths
8 A.P. Liu et al. / Wave Motion 49 (01) Fig. 8. A schematic illustration of dipolar motion in the metamaterial near the resonant frequency. in an order of the dimension of the unit-cell. To descrie the relative high-frequency dynamic ehaviors, microstructure kinetic variales are also needed. In the formulation, we assume the lead core is rigid ecause of its high modulus and the local displacements in the matrix can e approximated y linear series expansions in terms of quantities which are defined at the center of the element. The micro-deformations are assumed as follows: (1) In the inclusion of the cell (k, l), (r < a) u 1 = u (4a) () In the matrix of the cell (k, l), (r ) u = v Φ Θ r Φ Θ (4) where u 1 = u 1(k,l) u 1(k,l) T, u = u (k,l) u (k,l) T, (k,l) Φ (k,l) = (k,l) (k,l) u = u (k,l) u (k,l) T, v = v (k,l) T, Θ = cos θ sin θ T,, Φ = v (k,l) (k,l) (k,l) (k,l) (k,l), = 1 1 sin θ cos θ. where u, v, Φ and Φ are gloal variales in functions of (X, X, t), the local displacements u 1 and u are functions of (X, X, x, x, t), or in other words, functions of (X, X, r, θ, t). Physical interpretation of the terms in Eq. (4) is that gloal displacements u and v are the displacements of the center of the inclusion and the center of the matrix, respectively; while Φ represents micro-deformation in the area within the inner oundary of the matrix, and Φ represents micro-deformation in the matrix. It should e mentioned that the additional displacement components are necessary to capture the dipolar motion in the metamaterial, which is different from the microstructure continuum theory [9]. In the microstructure continuum theory, additional microstructure variales were only introduced to capture the micro-deformation in the microstructure. In principle, the oundary condition etween the unit-cell and the neighoring cells should e satisfied on every point of the oundaries. However, the approximated local field, assumed in Eq. (4), cannot exactly satisfy the point-to-point continuity at the oundaries. In this study, relaxation oundary continuity conditions, which are defined as averaged displacement continuity conditions at the interfaces of the cells, are suggested as: d / d / d / d / u (k1,l) α u (k,l1) α x = d x = d u (k,l) α u (k,l) α x = d x = d dx = 0 dx = 0 where uα (k1,l) and u (k,l1) α represent displacement components of the matrix in the cells (k 1, l) and (k, l 1). Eq. (5a) represents the averaged displacement continuity on the oundary etween the cell (k, l) and the cell (k 1, l); and Eq. (5) represents the averaged displacement continuity on the oundary etween the cell (k, l) and the cell (k, l 1). This approximate is much more reasonale for relatively low frequency cases. Sustituting the local displacement in the matrix Eq. (4) into Eq. (5a) results in (5a) (5)
9 418 A.P. Liu et al. / Wave Motion 49 (01) v (k1,l) α v (k,l) α 1d 1 1 ξ (k1,l) α d ξ (k1,l) α d (k1,l) α (k,l) α (k,l) α (k,l) α = 0 (6a) where ξ = 1d. d Similarly, another continuous condition can e otained from Eq. (5) as v (k,l1) α d v (k,l) α d (k,l1) α d (k,l) α ξ 1 ξ (k,l1) α (k,l1) α (k,l) α = 0. The corresponding local strain in the matrix can e otained as (k,l) α ε = u (7) (k,l) ε where is the differentiation operator in the local coordinate system (x, x ), ε ε (k,l) = is the strain tensor of the ε (k,l) ε (k,l) matrix and its components are given in Appendix A. Based on the displacement expressions in Eq. (4), the total kinetic energy density in the cell (k, l) can e calculated as (6) T (k,l) ave = 1 A c (T 1(k,l) T (k,l) ) where T 1(k,l) and T (k,l) are the kinetic energies in the inclusion and matrix, respectively. On the other hand, the strain deformation energy in the spring within the cell (k, l) can e otained as W s(k,l) = K v (k,l) u (k,l) 1 (k,l) 1 (k,l) K v (k,l) u (k,l) (k,l) (k,l) (8) (9) and the strain deformation energy in the matrix for the plane stress prolem is W m(k,l) E = ε (k,l) ε (k,l) ν ε (k,l) ε (k,l) ( ) A The total strain deformation energy averaged over the volume of cell (k, l) yields ( ) ε (k,l) da. (10) W (k,l) ave = 1 A c W s(k,l) W m(k,l). (11) To otain a continuum model, we now introduce fields that are continuous in the gloal coordinate system X and X, and the values at X = X l and X = X k coincide with those in the actual micro (local) field variales at the center of the cell. Therefore, we can consider the strain energy density W(X, X, t) and the kinetic energy density T(X, X, t) as continuous functions. Based on Eq. (8), the kinetic energy density T(X, X, t) in the continuum field is T (X, X, t) = T (k,l) ave = 1 A c (T 1(k,l) T (k,l) ) (1) and the strain energy density W(X, X, t) is W (X, X, t) = W (k,l) ave = 1 A c W s(k,l) W m(k,l) (1) where energy densities T (X, X, t) and W (X, X, t) are given in Appendix B. At the same time, the continuity conditions in Eq. (6) can e written in the terms of continuous variales y considering the field variales as continuous functions of X and X, as S α (X, X, t) = v α 1 X d α d 1 1 ξ ξ α α d α d X α α X = 0 (14a) X
10 A.P. Liu et al. / Wave Motion 49 (01) and S α (X, X, t) = v α ξ 1 ξ α X d α d α d X α X α d α = 0 (14) X where α =,. Considering a fixed regular region V of the medium, the displacement equations of motion can e otained y employing Hamilton s principle for independent variales in V and a specified time interval t 0 t t 1 as t1 t1 δ FdtdV δw 1 dt = 0 (15) t 0 t 0 V where F = T W, δw 1 is the variation of the work done y external forces, and dv is the scalar volume element. For the current model, there are no existing external forces and the continuity conditions can e considered as susidiary conditions through the use of Lagrangian multipliers, so the prolem can e redefined as t1 δ FdtdV = 0 (16) with t 0 V F = T W (Γ α S α Γ α S α ) α= where the Lagrangian multipliers Γ α and Γ α are functions of X, X and t. Since the function F as given in Eq. (17) depends only on the gloal field variales (X, X, t) and their first order derivatives, the system of the Euler equations can e written as r=1 p r F F = 0 (18) fs f s p r where f s (s = 1,,..., 16) represent the sixteen dependent variales u α, v α, α, α, α, α, Γ α and Γ α (α =, ); and p r (r = 1,, ) are the spatial variales X, X and time variale t. A system of sixteen governing equations of motion can e otained from the Eq. (18), which is given in Appendix C. (17) 4. Numerical simulation and discussions 4.1. Model verification To verify the proposed continuum model, let us first consider a two-dimensional elastic metamaterial with a cylinder heavy core coated with a circular soft layer and emedded in matrix in a square lattice array. The microstructure geometry parameters are listed in Tale 1 and the material parameters are provided in Tale. Based on the simplified model for the continuous systems, spring constants for the D plane stress prolem with the circle coated metamaterial can e numerically determined as K = K = N/m. It should e mentioned that for the current microstructure geometry, when the coating material is incompressile (the Poisson s ratio of the coat material ν is aout 0.5), the analytical solution for the spring constants [41] is also availale as K = K = 40πµ (a ) 5(a ) (/a) 9( a ), where µ is the shear modulus of the coating material µ = E. The good agreement aout the spring constant s prediction etween the two methods shows the (1ν ) accuracy of the current numerical method. For a longitudinal wave propagation along X direction in the metamaterial with infinite dimension in X direction, the availale kinematic variales include multi-displacement variales u, v, the microstructure field variales, and the Lagrangian multiplier Γ. Therefore, five governing equations can e otained as ρ 1 A 1 ü K (u v ) = 0 (19) ρ A v K (v u ) A c Γ X = 0 (0) ρ J 1 ρ (J 60 J 1 ) K E E J6 (1 ν ) Ac d J4 J7 ξ d 1 ξ Γ A c d ξ 1 ξ Γ X = 0 (1)
11 40 A.P. Liu et al. / Wave Motion 49 (01) Fig. 9. Comparison of the normalized dispersion curves y the current model, the reduced single-displacement model and the FE simulation for longitudinal wave propagation. ρ J1 I m 60 J ρ (J 60 J 1 ) E A J 4 E J 7 J 7 d A c d v X d (1 ν ) ξ 1 ξ d J6 Γ ξ 1 ξ d Ac d X ξ 1 ξ 1 = 0 () X d X d = 0 X where ξ = d for the circle-coated elastic metamaterial. For the longitudinal wave propagating in the X d direction, the continuum wave fields can e defined as u, v,,, Γ = {B1, B, B, B 4, B 5 } exp [i (q X ωt)] (4) Γ () where i = 1, B 1, B, B, B 4 and B 5 are constant amplitudes, q is the wave-numer in the X direction, and ω is the angular frequency. Sustitution of (4) into the Eqs. (19) () yields five homogeneous equations for B 1, B, B, B 4 and B 5. For a nontrivial set of solutions the determinant of the coefficients must vanish to yield the dispersion relation. The exact wave dispersion curves can also e calculated y using commercial finite element (FE) software, ANSYS In the FE model, Plane-8 element is used to model the elastic medium, and Mass-1 and Comine-14 elements are used to model the spring-mass system. After properly applying oundary conditions for the model, modal analysis is conducted to otain the natural frequencies, from which the dispersion curve can e otained. Fig. 9 shows comparison of the normalized dispersion curves of the longitudinal wave predicted y the current multi-displacement microstructure continuum model, the reduced single-displacement continuum model, and the FE simulation, where ω 0L = K /m 1 = rad/s and the resonant frequency is f 0L = ω 0L /(π) = 797 Hz. The single-displacement model is reduced from the current model y removing the additional gloal displacement variale. It is found that the multi-displacement microstructure continuum model can give an excellent prediction for the acoustic wave mode and is also reasonaly good for the optic wave mode for q d < 0.9. The and gap ehavior of the elastic metamaterial can e accurately predicted y using the current model, which is very useful for the design of the desirale elastic metamaterial. It is also very interesting to note that the and gap ehavior cannot e captured when only a single-displacement (conventional) continuum model is adopted. This is understandale ecause no additional degrees of freedom are used to capture the local dipolar resonance phenomenon. In addition, it is noticed that the current model has difficulty to capture the optic wave mode when the frequency is around ω/ω 0L = 8.48 (f = 6757 Hz). Based on the eigenmode analysis at this frequency, it can e found that the resonance deformation is focus on the coating layer and no motions occur in the core and host matrix, which is similar to the quadrupolar resonance. To capture the quadrupolar resonance in the metamaterial, additional microstructure kinematic variales in the coating layer may e needed.
12 A.P. Liu et al. / Wave Motion 49 (01) Fig. 10. The displacement amplitude ratio of the matrix to the core in function of the frequency. Fig. 11. Comparison of the normalized dispersion curves otained y the current model and the FE simulation for transverse wave propagation. Fig. 10 shows the wave amplitude ratio of the matrix to the core in a function of the frequency predicted y the current model and the finite element method (FEM). The very good agreement etween the current model and FEM is oserved for the local displacement prediction in the elastic metamaterial, which is impossile for the conventional continuum models. Around f = 797 Hz, there exists a resonance in the metamaterial in which the displacement of the core is very large and experience a very sharp change from the in-phase state to the out-of-phase state. A large displacement ratio enhancement can e oserved in the in-phase resonance region. From eigenmode analysis, we can find that the core inclusion moves opposite to the motion in the matrix in the and gap frequency range, which cannot e captured through the conventional continuum model. Fig. 11 shows a comparison of the normalized dispersion curves predicted y the current model and the FE simulation for transverse shear wave propagation. The material constants are the same as those used in Fig. 9, where ω 0T = K /m 1 = rad/s = ω 0L for the isotropic elastic metamaterial. The very reasonale agreement etween the current model and the FE simulation is also oserved for oth the acoustic wave mode and the optic wave mode. Also, it is noticed that the current model is invalid to capture the quadrupolar resonant mechanism when the frequency is around ω/ω 0L = 8.48 (f = 6757 Hz).
13 4 A.P. Liu et al. / Wave Motion 49 (01) Fig. 1. Comparison of the normalized dispersion curves for the anisotropic metamaterial for transverse wave propagation along X direction. Fig. 1. Comparison of the normalized dispersion curves for the anisotropic metamaterial for transverse wave propagation along X direction. 4.. Wave propagation in anisotropic elastic metamaterials It is well known that for isotropic metamaterials, wave propagation ehavior is identical for any wave propagation direction. For anisotropic elastic metamaterials, wave propagation ehavior is different along different directions and is direction-dependent. Figs. 1 and 1 show the comparison of the normalized dispersion curves for the transverse shear wave propagation along X and X directions for the anisotropic elastic metamaterial, respectively. In the figures, the material properties are the same as in Tale. The microstructure geometry parameters are shown in Tale 1 with an elliptical coating medium, from which the spring constants are numerically otained as K = kg/m, K = kg/m. The resonant angular frequencies are ω 0T = K /m 1 and ω 0T = K /m 1, respectively. It is found that the current model can provide good prediction of dispersion curves for oth the acoustic wave mode and the optic wave mode along different wave propagation directions. The anisotropic wave propagation ehavior can e oserved from the different wave dispersion curves along different directions. From the numerical simulation, it is also concluded that the and gap frequency regime is within the same frequency when the effective mass density ecomes negative, which confirms wave mechanism in the current anisotropic metamaterial is caused y the dipolar wave motion.
14 A.P. Liu et al. / Wave Motion 49 (01) Conclusions In this paper, an anisotropic elastic metamaterial made of lead cylinders coated with elliptical ruers in an epoxy matrix is considered. An anisotropic effective mass density tensor is first numerically determined. To descrie gloal dynamic ehavior in the elastic metamaterial, a new multi-displacement microstructure continuum model is developed to otain the macroscopic governing equations of the anisotropic elastic metamaterial. The current model is verified through comparison of wave dispersion curves predicted y the current model and the finite element simulation for oth longitudinal and transverse shear waves. Very good agreement is oserved in oth the acoustic and optic wave modes. The proposed model may provide an efficient tool for modeling of the elastic metamaterial with complex microstructures. Acknowledgments This work was supported y the National Natural Science Foundation of China (110701, 10800, ) and Natural Science Foundation EAGER program (107569). Appendix A. Components of the local stain tensor The detailed components of the local strain tensor are: ε (k,l) = (k,l) ε (k,l) = (k,l) ( cos θ) x ( cos θ) x (k,l) (k,l) (k,l) (k,l) ε (k,l) = ε (k,l) = 1 (k,l) (k,l) 1 x 1 ( sin θ) (k,l) (k,l) 1 x ( sin θ) x ( sin θ) ( cos θ) x ( cos θ) x (k,l) (k,l) (k,l) (k,l) (k,l) (k,l) (k,l) (k,l) 1 ( sin θ) x (k,l) (k,l). (A.1) (A.) (A.) Appendix B. Kinetic energy density and strain energy density The kinetic energy density T(X, X, t) in the continuum fields is T (X, X, t) = T (k,l) ave 1 ρ = 1 A c (T 1(k,l) T (k,l) ) = 1 A c 1 ρ 1A 1 u u α= A ( v α ) J 1 I m J 60 ( α ) J I m J 70 ( α ) J 1 ( α ) J ( α ) (J 60 J 1 ) α α (J 70 J ) α α (B.1) where I m = x A da, I m = x da, J 1 = cos θ da, A A J = sin θ da, J 60 = r cos θ da, J 70 = r sin θ da. A A A And the strain energy density W(X, X, t) in the continuum fields is (k, l) W (X, X, t) = W ave = 1 W s(k,l) W m(k,l) A c = 1 K (v u ) 1 1 K (v u ) A c E A A J 4 ( ) J 5 J 7 J 6 J 6 J8
15 44 A.P. Liu et al. / Wave Motion 49 (01) E ν A J 7 J 8 J 6 E 1 ν A 4 J 4 J 6 J6 1 4 J 5 J 6 J 6 J 6 J7 J 8 J 7 J 8 (B.) where J 4 = J 7 = A A r sin4 θ da, J 5 = A r cos4 θ da, J 6 = r sin θ da, J 8 = r cos θ da. A A r sin θ cos θ da, Appendix C. Governing equations The sixteen governing equations of motion are: ρ 1 A 1 ü K (u v ) = 0 ρ 1 A 1 ü K (u v ) = 0 (C.1) (C.) ρ A v K (v u ) A c Γ X A c Γ X = 0 (C.) ρ A v K (v u ) A c Γ X A c Γ X = 0 (C.4) ρ J 1 ρ (J 60 J 1 ) 1 K E E ν J7 J 6 E A c d (1 ν ) ξ 1 ξ d Γ A c d J4 J7 J6 J6 ξ 1 ξ Γ J6 ρ J 1 ρ (J 60 J 1 ) K E E ν J6 E J4 (1 ν ) J 7 J 7 Ac ξ 1 ξ d Γ A c ρ J ρ (J 70 J ) 1 K E E ν J6 E (1 ν ) d J 8 J 8 1 Ac d 1 1 ξ ξ X = 0 (C.5) J6 d J6 J5 ξ Γ 1 ξ J6 1 d Γ A c d J5 J8 1 1 ξ ρ J ρ (J 70 J ) K E E ν J8 J 6 E J6 (1 ν ) J6 ξ X = 0 (C.6) Γ X = 0 (C.7)
16 A c 1 d 1 1 ξ ξ A.P. Liu et al. / Wave Motion 49 (01) Γ A c 1 d d ρ J1 I m 60 J ρ (J 60 J 1 ) E J 7 J 7 E E ν J6 J6 (1 ν ) A c ξ 1 ξ d 1 1 ξ ξ A J 4 Γ X = 0 (C.8) A J 7 J 8 J6 1 Γ A c d d ρ J1 I m J 60 ρ (J 60 J 1 ) E ξ 1 ξ d Γ X = 0 (C.9) J6 E ν J6 E A (1 ν ) J4 J 6 J7 J 7 J 7 J 8 A c ξ d 1 ξ 1 Γ A c d d ρ J I m J 70 ρ (J 70 J ) E ξ 1 ξ d ] [J 6 E ν [J 6 E ] A (1 ν ) J5 J 6 J8 J 8 J 8 J ξ 1 d 1 1 ξ A c 1 Γ Ac d ξ d ξ ρ J I m 70 J ρ (J 70 J ) E J 8 J 8 E ν E J6 (1 ν ) J ξ A c 1 Γ Ac d ξ A J 5 Γ X = 0 (C.10) d A J 8 J 7 J6 v ξ 1 ξ X d d d X v X d v 1 X d v 1 X d ξ 1 ξ d 1 1 ξ ξ 1 1 ξ ξ 1 d 1 1 ξ d ξ d d X d d X d d X X X X X d Γ X = 0 (C.11) Γ X = 0 (C.1) X = 0 (C.1) d = 0 (C.14) X d = 0 (C.15) X d = 0. (C.16) X References [1] X. Hu, C.T. Chan, J. Zi, Two-dimensional sonic crystals with Helmholtz resonators, Phys. Rev. E 71 (005) (R). [] X. Hu, C.T. Chan, Refraction of water waves y periodic cylinder arrays, Phys. Rev. Lett. 95 (005) [] Z. Liu, X. Zhang, Y. Mao, Y.Y. Zhu, Z. Yang, C.T. Chan, P. Sheng, Locally resonant sonic materials, Science 89 (000) [4] S. Yang, J.H. Page, Z. Liu, M.L. Cowan, C.T. Chan, P. Sheng, Ultrasound tunneling through D phononic crystals, Phys. Rev. Lett. 88 (00)
17 46 A.P. Liu et al. / Wave Motion 49 (01) [5] Z. Hou, X. Fu, Y. Liu, Calculational method to study the transmission properties of phononic crystals, Phys. Rev. B 70 (004) [6] N. Fang, D. Xi, J. Xu, M. Amati, W. Srituravanich, C. Sun, X. Zhang, Ultrasonic metamaterials with negative modulus, Nature Mater. 5 (006) [7] S. Gonella, A.C. To, W.K. Liu, Interplay etween phononic andgaps and piezoelectric microstructures for energy harvesting, J. Mech. Phys. Solids 57 (009) [8] X.N. Liu, G.K. Hu, G.L. Huang, C.T. Sun, Wave propagation characterization and design of two-dimensional elastic chiral metacomposite, J. Vi. Acoust. 0 (011). [9] M. Kafesaki, E.N. Economou, Multiple scattering theory for D periodic acoustic composites, Phys. Rev. B 60 (1999) [10] X. Hu, K.M. Ho, C.T. Chan, J. Zi, Homogenization of acoustic metamaterials of Helmholtz resonators in fluid, Phys. Rev. B 77 (008) [11] P. Sheng, Introduction to Wave Scattering, Localization, and Mesoscopic Phenomena, Academic, San Diego, [1] Z. Liu, C.T. Chan, P. Sheng, Analytic model of phononic crystals with local resonances, Phys. Rev. B 71 (005) [1] J. Mei, Z. Liu, W. Wen, P. Sheng, Effective mass density of fluid solid composites, Phys. Rev. Lett. 96 (006) [14] Y. Wu, Y. Lai, Z.Q. Zhang, Effective medium theory for elastic metamaterials in two dimensions, Phys. Rev. B 76 (007) 051. [15] Z. Yang, J. Mei, M. Yang, N.H. Chan, P. Sheng, Memrane-type acoustic metamaterial with negative dynamic mass, Phys. Rev. Lett. 101 (008) [16] X.M. Zhou, G.K. Hu, Analytic model of elastic metamaterials with local resonances, Phys. Rev. B 79 (009) [17] Y. Ding, Z. Liu, C. Qiu, J. Shi, Metamaterial with simultaneously negative ulk modulus and mass density, Phys. Rev. Lett. 99 (007) [18] M.G. Imhof, Multiple multipole expansions for elastic scattering, J. Acoust. Soc. Am. 100 (1996) [19] G. Wang, X. Wen, J. Wen, L. Shao, Y. Liu, Two-dimensional locally resonant phononic crystals with inary structures, Phys. Rev. Lett. 9 (004) [0] A. Khelif, P.A. Deymier, B. Djafari-Rouhani, J.O. Vasseur, L. Dorzynski, Two-dimensional phononic crystal with tunale narrow pass and: application to a waveguide with selective frequency, J. Appl. Phys. 94 (00) 108. [1] D. Torrent, J. Sánchez-Dehesa, Anisotropic mass density y two-dimensional acoustic metamaterials, New J. Phys. 10 (008) [] X. Ao, C.T. Chan, Far-field image magnification for acoustic waves using anisotropic acoustic metamaterials, Phys. Rev. E 77 (008) [] J.R. Willis, Exact effective relations for dynamics of a laminated ody, Mech. Mater. 41 (009) [4] J.R. Willis, Effective constitutive relations for waves in composites and metamaterials, Proc. R. Soc. A (011), doi: /rspa [5] A.L. Shuvalov, A.A. Kutsenko, A.N. Norris, O. Poncelet, Effective Willis constitutive equations for periodically stratified anisotropic elastic media, Proc. R. Soc. A 467 (011) [6] S. Nemat-Nasser, J.R. Willis, A. Srivastava, A.V. Amirkhizi, Homogenization of periodic elastic composites and locally resonant sonic materials, Phys. Rev. B 8 (011) [7] M. Hirsekorn, P.P. Delsanto, Elastic wave propagation in locally resonant sonic material: comparison etween local interaction simulation approach and modal analysis, J. Appl. Phys. 99 (006) [8] Y. Gu, X. Luo, H. Ma, Low frequency elastic wave propagation in two dimensional locally resonant phononic crystal with asymmetric resonator, J. Appl. Phys. 105 (009) [9] Cz. Wozniak, Refined macro-dynamics of periodic structures, Arch. Mech. 45 (199) [0] G.W. Milton, J.R. Willis, On modifications of Newton s second law and linear continuum elastodynamics, Proc. R. Soc. A 46 (007) [1] G.W. Milton, New metamaterials with macroscopic ehavior outside that of continuum elastodynamics, New J. Phys. 9 (007) 59. [] J.J. Aliert, P. Seppecher, F. Dell Isola, Truss modular eams with deformation energy depending on higher displacement gradients, Math. Mech. Solids 8 (00) [] E. Cosserat, F. Cosserat, Théorie des corps déformales, A. Hermann et Fils, Paris, [4] A.C. Eringen, E.S. Suhui, Nonlinear theory of simple micro-elastic solids-i, Int. J. Eng. Sci. (1964) [5] R.A. Toupin, Elastic materials with couple-stresses, Arch. Ration. Mech. Anal. 11 (196) [6] A.C. Eringen, Microcontinuum Field Theories I: Foundations and Solids, Springer-Verlag, New York, [7] R.D. Mindlin, Micro-structure in linear elasticity, Arch. Ration. Mech. Anal. 16 (1964) [8] C.T. Sun, J.D. Achenach, G. Herrmann, Continuum theory for a laminated medium, J. Appl. Mech. 5 (1968) [9] G.L. Huang, C.T. Sun, Continuum modeling of solids with micro/nanostructures, Philos. Mag. 87 (007) [40] R. Zhu, H.H. Huang, G.L. Huang, C.T. Sun, Microstructure continuum modeling of an elastic metamaterial, Int. J. Eng. Sci. 49 (011) [41] J.R. Barer, Elasticity, third ed., Springer, Dordrecht, 010.
Band Gaps in a Multiresonator Acoustic Metamaterial
G. L. Huang Department of Systems Engineering, University of Arkansas at Little Rock, Little Rock, AR 7 e-mail: glhuang@ualr.edu C. T. Sun School of Aeronautics and Astronautics, Purdue University, W.
More informationWave Motion. On the limit and applicability of dynamic homogenization. Ankit Srivastava a,, Sia Nemat-Nasser b. h i g h l i g h t s.
Wave Motion ( ) Contents lists available at ScienceDirect Wave Motion journal homepage: www.elsevier.com/locate/wavemoti On the limit and applicability of dynamic homogenization Ankit Srivastava a,, Sia
More informationTheoretical study of subwavelength imaging by. acoustic metamaterial slabs
Theoretical study of subwavelength imaging by acoustic metamaterial slabs Ke Deng,2, Yiqun Ding, Zhaojian He, Heping Zhao 2, Jing Shi, and Zhengyou Liu,a) Key Lab of Acoustic and Photonic materials and
More informationRESONANCE-COUPLING EFFECT ON BROAD BAND GAP FORMATION AND SOUND ABSORPTION IN LOCALLY RESONANT SONIC METAMATERIALS WITH WOODPILE STRUCTURE
RESONANCE-COUPLING EFFECT ON BROAD BAND GAP FORMATION AND SOUND ABSORPTION IN LOCALLY RESONANT SONIC METAMATERIALS WITH WOODPILE STRUCTURE Wang Yuren Key Laboratory of Microgravity, Institute of Mechanics,
More informationBand gaps in a phononic crystal constituted by cylindrical dots on a homogeneous plate
Band gaps in a phononic crystal constituted by cylindrical dots on a homogeneous plate B. Djafari-Rouhani, Y. Pennec, H. Larabi, J. Vasseur and A.-C. Hladky IEN, UR CNRS 852, avenue Poincaré, BP 669, 59652
More informationTwo-dimensional ternary locally resonant phononic crystals with a comblike coating
Two-dimensional ternary locally resonant phononic crystals with a comblike coating Yan-Feng Wang, Yue-Sheng Wang,*, and Litian Wang Institute of Engineering Mechanics, Beijing Jiaotong University, Beijing,
More informationNoise Shielding Using Acoustic Metamaterials
Commun. Theor. Phys. (Beijing, China) 53 (2010) pp. 560 564 c Chinese Physical Society and IOP Publishing Ltd Vol. 53, No. 3, March 15, 2010 Noise Shielding Using Acoustic Metamaterials LIU Bin ( Ê) and
More informationAn Elastic Metamaterial Beam for Broadband Vibration Suppression
An Elastic Metamaterial Beam for Broadband Vibration Suppression R. Zhu a, G. K. Hu b, M. Reynolds c and G. L. Huang *,a a Department of Systems Engineering, University of Arkansas at Little Rock, Little
More informationControlling elastic wave with isotropic transformation materials
Controlling elastic wave with isotropic transformation materials Zheng Chang, Jin Hu, a, Gengkai Hu, b, Ran Tao and Yue Wang School of Aerospace Engineering, Beijing Institute of Technology, 0008,Beijing,
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationMethod for retrieving effective properties of locally resonant acoustic metamaterials
Method for retrieving effective properties of locally resonant acoustic metamaterials Vladimir Fokin, Muralidhar Ambati, Cheng Sun, and Xiang Zhang* Nano-scale Science and Engineering Center, University
More informationSuper-resolution imaging by resonant tunneling in anisotropic acoustic metamaterials
Super-resolution imaging by resonant tunneling in anisotropic acoustic metamaterials Aiping Liu and Xiaoming Zhou a) Key Laboratory of Dynamics and Control of Flight Vehicle, Ministry of Education and
More informationAcoustic pressure characteristic analysis in cavity of 2-D phononic crystal
Journal of Engineering Technology and Education, Vol. 9, No. June 1, pp. 115-11 Acoustic pressure characteristic analysis in cavity of -D phononic crystal Jia-Yi Yeh 1, Jiun-Yeu Chen 1 Department of Information
More informationA Lumped Model for Rotational Modes in Phononic Crystals
A Lumped Model for Rotational Modes in Phononic Crystals Pai Peng, Jun Mei and Ying Wu Division of Mathematical and Computer Sciences and Engineering, King Abdullah University of Science and Technology
More informationPHYSICAL REVIEW B 71,
Coupling of electromagnetic waves and superlattice vibrations in a piezomagnetic superlattice: Creation of a polariton through the piezomagnetic effect H. Liu, S. N. Zhu, Z. G. Dong, Y. Y. Zhu, Y. F. Chen,
More informationA direct evaluation of the Fabric Tensor in anisotropic porous media
A direct evaluation of the Fabric Tensor in anisotropic porous media Maria Cristina Pernice 1, Luciano Nunziante 1, Massimiliano Fraldi 1,2 1 Department of Structural Engineering, University of Naples
More informationGraduate School of Engineering, Kyoto University, Kyoto daigaku-katsura, Nishikyo-ku, Kyoto, Japan.
On relationship between contact surface rigidity and harmonic generation behavior in composite materials with mechanical nonlinearity at fiber-matrix interface (Singapore November 2017) N. Matsuda, K.
More informationarxiv: v3 [cond-mat.mtrl-sci] 19 Jan 2009
Negative Effective Density in An Acoustic Metamaterial Sam Hyeon Lee, 1 Choon Mahn Park, 2 Yong Mun Seo, 3 Zhi Guo Wang, 4 and Chul Koo Kim 1 1 Institute of Physics and Applied Physics, Yonsei University,
More informationEffects of mass layer dimension on a finite quartz crystal microbalance
University of Neraska - Lincoln DigitalCommons@University of Neraska - Lincoln Mechanical & Materials Engineering Faculty Pulications Mechanical & Materials Engineering, Department of Summer 8-2011 Effects
More informationMetamaterials with tunable dynamic properties
Metamaterials with tunable dynamic properties Varvara Kouznetsova Marc Geers 6 October 2015 Mechanics of Materials Project aim development of new generation mechanical metamaterials with adaptive, tunable
More informationAnalytical coupled vibroacoustic modeling of membranetype acoustic metamaterials: membrane model
Analytical coupled vibroacoustic modeling of membranetype acoustic metamaterials: membrane model Yangyang Chen and Guoliang Huang a) Department of Systems Engineering, University of Arkansas at Little
More informationSIMILARITY METHODS IN ELASTO-PLASTIC BEAM BENDING
Similarity methods in elasto-plastic eam ending XIII International Conference on Computational Plasticity Fundamentals and Applications COMPLAS XIII E Oñate, DRJ Owen, D Peric and M Chiumenti (Eds) SIMILARIT
More informationA simple plane-strain solution for functionally graded multilayered isotropic cylinders
Structural Engineering and Mechanics, Vol. 24, o. 6 (2006) 000-000 1 A simple plane-strain solution for functionally graded multilayered isotropic cylinders E. Pan Department of Civil Engineering, The
More informationFlute-Model Acoustic Metamaterials with Simultaneously. Negative Bulk Modulus and Mass Density
Flute-Model Acoustic Metamaterials with Simultaneously Negative Bulk Modulus and Mass Density H. C. Zeng, C. R. Luo, H. J. Chen, S. L. Zhai and X. P. Zhao * Smart Materials Laboratory, Department of Applied
More informationON THE MODELING AND SIMULATION OF AN ACOUSTIC CLOAK
5 th International Conference Advanced Composite Materials Engineering COMAT 014 16-17 October 014, Braşov, Romania ON THE MODELING AND SIMULATION OF AN ACOUSTIC CLOAK Veturia Chiroiu 1, Rodica Ioan 1,,
More informationPhononic Crystals: Towards the Full Control of Elastic Waves propagation OUTLINE
Phononic Crystals: Towards the Full Control of Elastic Waves propagation José Sánchez-Dehesa Wave Phenomena Group, Department of Electronic Engineering, Polytechnic University of Valencia, SPAIN. OUTLINE
More informationResonant modal group theory of membrane-type acoustical metamaterials for low-frequency sound attenuation
Eur. Phys. J. Appl. Phys. (2015) 71: 30504 DOI: 10.1051/epjap/2015150310 Resonant modal group theory of membrane-type acoustical metamaterials for low-frequency sound attenuation Fuyin Ma, Jiu Hui Wu,
More informationNegative effective dynamic mass-density and stiffness: Micro-architecture and phononic transport in periodic composites
AIP ADVANCES 1, 041502 (2011) Negative effective dynamic mass-density and stiffness: Micro-architecture and phononic transport in periodic composites Sia Nemat-Nasser and Ankit Srivastava a Department
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationBuckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads
NASA Technical Paper 3659 Buckling Behavior of Long Symmetrically Laminated Plates Sujected to Shear and Linearly Varying Axial Edge Loads Michael P. Nemeth Langley Research Center Hampton, Virginia National
More informationA Stable and Convergent Finite Difference Scheme for 2D Incompressible Nonlinear Viscoelastic Fluid Dynamics Problem
Applied and Computational Mathematics 2018; (1): 11-18 http://www.sciencepulishinggroup.com/j/acm doi: 10.11648/j.acm.2018001.12 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Stale and Convergent
More informationDevelopment and Application of Acoustic Metamaterials with Locally Resonant Microstructures
Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures AFOSR grant #FA9550-10-1-0061 Program manager: Dr. Les Lee PI: C.T. Sun Purdue University West Lafayette, Indiana
More informationCRACK-TIP DIFFRACTION IN A TRANSVERSELY ISOTROPIC SOLID. A.N. Norris and J.D. Achenbach
CRACK-TIP DIFFRACTION IN A TRANSVERSELY ISOTROPIC SOLID A.N. Norris and J.D. Achenbach The Technological Institute Northwestern University Evanston, IL 60201 ABSTRACT Crack diffraction in a transversely
More informationACOUSTIC PERFORMANCE OF PERIODIC COMPOSITE MATERIALS
ACOUSTIC PERFORMANCE OF PERIODIC COMPOSITE MATERIALS Abstract Gyani Shankar Sharma 1, Daniel Eggler 1, Herwig Peters 1, Nicole Kessissoglou 1, Alex Skvortsov 2, Ian MacGillivray 2 1 School of Mechanical
More informationApplication and analysis of phononic crystal energy harvesting devices
J. Eng. Technol. Educ. (013) 10(1): 18-6 March 013 Application and analysis of phononic crystal energy harvesting devices Department of Information Management, Chung Hwa University of Medical Technology.
More informationInterfacial effects in electromagnetic coupling within piezoelectric phononic crystals
Acta Mech Sin (29) 25:95 99 DOI 1.17/s149-8-21-y RESEARCH PAPER Interfacial effects in electromagnetic coupling within pieoelectric phononic crystals F. J. Sabina A. B. Movchan Received: 14 July 28 / Accepted:
More informationDesign and Control of Advanced Functional Systems with Non- Conventional Material Properties
Design and Control of Advanced Functional Systems with Non- Conventional Material Properties Dr Simon Pope Department of Automatic Control and Systems Engineering The University of Sheffield s.a.pope@sheffield.ac.uk
More informationPhononic Crystals. J.H. Page
Phononic Crystals J.H. Page University of Manitoba with Suxia Yang and M.L. Cowan at U of M, Ping Sheng and C.T. Chan at HKUST, & Zhengyou Liu at Wuhan University. We study ultrasonic waves in complex
More informationAdd-on unidirectional elastic metamaterial plate cloak
Add-on unidirectional elastic metamaterial plate cloak Min Kyung Lee *a and Yoon Young Kim **a,b a Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro, Gwanak-gu, Seoul,
More informationDesignable hybrid sonic crystals for transportation and division of acoustic images
Designable hybrid sonic crystals for transportation and division of acoustic images Zhaojian He 1a), Ke Deng 1b), Heping Zhao 1, and Xiaochun Li 2, 3 1. Department of Physics, Jishou University, Jishou
More informationRedirection of flexural waves in platonic crystal slabs
Redirection of flexural waves in platonic crystal slabs Penglin Gao Center for Composite Materials, Harbin Institute of Technology, Harbin, China. Wave Phenomena Group, Department of Electronic Engineering,
More informationExploiting pattern transformation to tune phononic band gaps in a two-dimensional granular crystal
Exploiting pattern transformation to tune phononic band gaps in a two-dimensional granular crystal The Harvard community has made this article openly available. Please share how this access benefits you.
More informationNegative-refraction imaging with two-dimensional phononic crystals
Negative-refraction imaging with two-dimensional phononic crystals Manzhu Ke, Zhengyou Liu,* Chunyin Qiu, Wengang Wang, and Jing Shi Department of Physics, Wuhan University, Wuhan 430072, People s Republic
More informationA method of the forward problem for magneto-acousto-electrical tomography
Technology and Health Care 24 (216) S733 S738 DOI 1.3233/THC-16122 IOS Press S733 A method of the forward problem for magneto-acousto-electrical tomography Jingxiang Lv a,b,c, Guoqiang Liu a,, Xinli Wang
More informationOverview. 1. What range of ε eff, µ eff parameter space is accessible to simple metamaterial geometries? ``
MURI-Transformational Electromagnetics Innovative use of Metamaterials in Confining, Controlling, and Radiating Intense Microwave Pulses University of New Mexico August 21, 2012 Engineering Dispersive
More informationEffect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue Ver. III (Jan. - Fe. 06), 5-65 www.iosrjournals.org Effect of Uniform Horizontal Magnetic Field on Thermal Instaility
More informationMoving screw dislocations in piezoelectric bimaterials
phys stat sol (b) 38 No 1 10 16 (003) / DOI 10100/pssb00301805 Moving screw dislocations in piezoelectric bimaterials Xiang-Fa Wu *1 Yuris A Dzenis 1 and Wen-Sheng Zou 1 Department of Engineering Mechanics
More informationPrecise method to control elastic waves by conformal mapping
THEORETICAL & APPLIED MECHANICS LETTERS 3, (3) Precise method to control elastic waves by conformal mapping Yongquan Liu, Wei Liu, and Xianyue Su a) LTCS and Department of Mechanics and Aerospace Engineering,
More informationBand gaps and the electromechanical coupling coefficient of a surface acoustic wave in a two-dimensional piezoelectric phononic crystal
Band gaps and the electromechanical coupling coefficient of a surface acoustic wave in a two-dimensional piezoelectric phononic crystal Tsung-Tsong Wu* Zin-Chen Hsu and Zi-ui Huang Institute of Applied
More informationSUPPLEMENTARY FIGURES
SUPPLEMENTARY FIGURES Supplementary Figure 1 Tuning the eigenfrequencies of the translational modes. A cutaway drawing of the metamaterial unit cell is shown in a. Here, h s denotes the silicone rubber
More informationMetadamping: An emergent phenomenon in dissipative metamaterials
Metadamping: An emergent phenomenon in dissipative metamaterials Mahmoud I. Hussein and Michael J. Frazier Department of Aerospace Engineering Sciences, University of Colorado Boulder, Boulder, CO 80309
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 9, 23 http://acousticalsociety.org/ ICA 23 Montreal Montreal, Canada 2-7 June 23 Structural Acoustics and Vibration Session psa: Measurement and Modeling of
More informationSurface resonant states and superlensing in acoustic metamaterials
Surface resonant states and superlensing in acoustic metamaterials Muralidhar Ambati, Nicholas Fang, Cheng Sun, and Xiang Zhang* Center for Scalable and ntegrated Nanomanufacturing (SNAM), University of
More informationSmall-Scale Effect on the Static Deflection of a Clamped Graphene Sheet
Copyright 05 Tech Science Press CMC, vol.8, no., pp.03-7, 05 Small-Scale Effect on the Static Deflection of a Clamped Graphene Sheet G. Q. Xie, J. P. Wang, Q. L. Zhang Abstract: Small-scale effect on the
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationUNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume, Numer 4/0, pp. 9 95 UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION
More informationPhysical and Biological Properties of Agricultural Products Acoustic, Electrical and Optical Properties and Biochemical Property
Physical and Biological Properties of Agricultural Products Acoustic, Electrical and Optical Properties and Biochemical Property 1. Acoustic and Vibrational Properties 1.1 Acoustics and Vibration Engineering
More informationBulletin of the Transilvania University of Braşov Vol 10(59), No Series III: Mathematics, Informatics, Physics, 83-90
Bulletin of the Transilvania University of Braşov Vol 10(59), No. 1-2017 Series III: Mathematics, Informatics, Physics, 83-90 GENERALIZED MICROPOLAR THERMOELASTICITY WITH FRACTIONAL ORDER STRAIN Adina
More informationComplete band gaps in two-dimensional phononic crystal slabs
Complete band gaps in two-dimensional phononic crystal slabs A. Khelif, 1 B. Aoubiza, 2 S. Mohammadi, 3 A. Adibi, 3 and V. Laude 1 1 Institut FEMTO-ST, CNRS UMR 6174, Université de Franche-Comté, Besançon,
More informationWAVE PROPAGATION IN PLATES WITH PERIODIC ARRAY OF IMPERFECT ACOUSTIC BLACK HOLES
WAVE PROPAGATION IN PLATES WITH PERIODIC ARRAY OF IMPERFECT ACOUSTIC BLACK HOLES Bing Han 1, Hongli Ji 2 and Jinhao Qiu 3 1 Yudao Street 29, Nanjing 210016, China, State Key Laboratory of Mechanics and
More informationMicrostructural Randomness and Scaling in Mechanics of Materials. Martin Ostoja-Starzewski. University of Illinois at Urbana-Champaign
Microstructural Randomness and Scaling in Mechanics of Materials Martin Ostoja-Starzewski University of Illinois at Urbana-Champaign Contents Preface ix 1. Randomness versus determinism ix 2. Randomness
More informationSurface-Plasmon-Polariton (SPP)-Like Acoustic Surface Waves on Elastic Metamaterials
Surface-Plasmon-Polariton (SPP)-Lie Acoustic Surface Waves on Elastic Metamaterials Ke Deng,2, Zhaojian He,2, Yiqun Ding, Heping Zhao 2, and Zhengyou Liu,* Key Lab of Artificial Micro- and Nano-structures
More information2018. Nonlinear free vibration analysis of nanobeams under magnetic field based on nonlocal elasticity theory
2018. Nonlinear free vibration analysis of nanobeams under magnetic field based on nonlocal elasticity theory Tai-Ping Chang National Kaohsiung First University of Science and Technology, Kaohsiung, Taiwan
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationModal analysis of waveguide using method of moment
HAIT Journal of Science and Engineering B, Volume x, Issue x, pp. xxx-xxx Copyright C 27 Holon Institute of Technology Modal analysis of waveguide using method of moment Arti Vaish and Harish Parthasarathy
More informationMotivation. Confined acoustics phonons. Modification of phonon lifetimes Antisymmetric Bulk. Symmetric. 10 nm
Motivation Confined acoustics phonons Modification of phonon lifetimes 0 0 Symmetric Antisymmetric Bulk 0 nm A. Balandin et al, PRB 58(998) 544 Effect of native oxide on dispersion relation Heat transport
More informationINTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 2, 2011
Volume, No, 11 Copyright 1 All rights reserved Integrated Pulishing Association REVIEW ARTICLE ISSN 976 459 Analysis of free virations of VISCO elastic square plate of variale thickness with temperature
More informationUC San Diego UC San Diego Electronic Theses and Dissertations
UC San Diego UC San Diego Electronic Theses and Dissertations Title Microstructurally Controlled Composites with Optimal Elastodynamic Properties Permalink https://escholarship.org/uc/item/9gw5w8tq Author
More informationPrediction of Elastic Constants on 3D Four-directional Braided
Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationSpatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Spatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media Karim G. Sabra, School of Mechanical Engineering,
More informationDual phononic and photonic band gaps in a periodic array of pillars deposited on a membrane
Excerpt from the Proceedings of the COMSOL Conference 010 Paris Dual phononic and photonic band gaps in a periodic array of pillars deposited on a membrane Y. Pennec, Y. El Hassouani, C. Li, B. Djafari
More informationComposite Plates Under Concentrated Load on One Edge and Uniform Load on the Opposite Edge
Mechanics of Advanced Materials and Structures, 7:96, Copyright Taylor & Francis Group, LLC ISSN: 57-69 print / 57-65 online DOI:.8/5769955658 Composite Plates Under Concentrated Load on One Edge and Uniform
More informationChapter 5. Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice
Chapter 5 Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice In chapter 3 and 4, we have demonstrated that the deformed rods, rotational rods and perturbation
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationPhysics Letters A 374 (2010) Contents lists available at ScienceDirect. Physics Letters A.
Physics Letters A 374 (2010) 1063 1067 Contents lists available at ScienceDirect Physics Letters A www.elsevier.com/locate/pla Macroscopic far-field observation of the sub-wavelength near-field dipole
More informationASEISMIC DESIGN OF TALL STRUCTURES USING VARIABLE FREQUENCY PENDULUM OSCILLATOR
ASEISMIC DESIGN OF TALL STRUCTURES USING VARIABLE FREQUENCY PENDULUM OSCILLATOR M PRANESH And Ravi SINHA SUMMARY Tuned Mass Dampers (TMD) provide an effective technique for viration control of flexile
More informationChapter 2 Fundamental Properties of Phononic Crystal
Chapter 2 Fundamental Properties of Phononic Crystal Yan Pennec and Bahram Djafari-Rouhani 2.1 Introduction to the Concept of Phononic Crystals and Their Band Structures The control and manipulation of
More informationLASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE
LASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE H. M. Al-Qahtani and S. K. Datta University of Colorado Boulder CO 839-7 ABSTRACT. An analysis of the propagation of thermoelastic waves
More information908. Dynamic responses of an elastic beam moving over a simple beam using modal superposition method
98. Dynamic responses of an elastic eam moving over a simple eam using modal superposition method Y. J. Wang, J. Shi, Y. Xia 3, School of Civil Engineering, Beijing Jiaotong University, Beijing 44, China,
More informationAcoustooptic Bragg Diffraction in 2-Dimensional Photonic Crystals
Acoustooptic Bragg Diffraction in 2-Dimensional Photonic Crystals Z.A. Pyatakova M.V. Lomonosov Moscow State University, Physics Department zoya.pyatakova@gmail.com Abstract. The paper shows that silicon-based
More informationVIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS
Journal of Engineering Science and Technology Vol. 12, No. 12 (217) 3398-3411 School of Engineering, Taylor s University VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS DILEEP
More informationNonlinear Free Vibration of Nanobeams Subjected to Magnetic Field Based on Nonlocal Elasticity Theory
Nonlinear Free Vibration of Nanobeams Subjected to Magnetic Field Based on Nonlocal Elasticity Theory Tai-Ping Chang 1 and Quey-Jen Yeh 1 Department of Construction Engineering, National Kaohsiung First
More informationSUPPLEMENTARY INFORMATION
Ultra-sparse metasurface for high reflection of low-frequency sound based on artificial Mie resonances Y. Cheng, 1,2 C. Zhou, 1 B.G. Yuan, 1 D.J. Wu, 3 Q. Wei, 1 X.J. Liu 1,2* 1 Key Laboratory of Modern
More informationThe effect of rigidity on torsional vibrations in a two layered poroelastic cylinder
Int. J. Adv. Appl. Math. and Mech. 3(1) (2015) 116 121 (ISSN: 2347-2529) Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics The effect of rigidity on
More informationSuper-reflection and Cloaking Based on Zero Index Metamaterial
Super-reflection and Cloaking Based on Zero Index Metamaterial Jiaming Hao, Wei Yan, and Min Qiu Photonics and Microwave ngineering, Royal Institute of Technology (KTH), lectrum 9, 164 4, Kista, Sweden
More informationNiiranen, Jarkko; Khakalo, Sergei; Balobanov, Viacheslav Strain gradient elasticity theories in lattice structure modelling
Powered by TCPDF (www.tcpdf.org) This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Niiranen, Jarkko; Khakalo, Sergei;
More informationEE 5344 Introduction to MEMS CHAPTER 6 Mechanical Sensors. 1. Position Displacement x, θ 2. Velocity, speed Kinematic
I. Mechanical Measurands: 1. Classification of main types: EE 5344 Introduction MEMS CHAPTER 6 Mechanical Sensors 1. Position Displacement x, θ. Velocity, speed Kinematic dx dθ v =, = ω 3. Acceleration
More informationIranian Journal of Mathematical Sciences and Informatics Vol.2, No.2 (2007), pp 1-16
Iranian Journal of Mathematical Sciences and Informatics Vol.2, No.2 (2007), pp 1-16 THE EFFECT OF PURE SHEAR ON THE REFLECTION OF PLANE WAVES AT THE BOUNDARY OF AN ELASTIC HALF-SPACE W. HUSSAIN DEPARTMENT
More informationSecond-gradient theory : application to Cahn-Hilliard fluids
Second-gradient theory : application to Cahn-Hilliard fluids P. Seppecher Laboratoire d Analyse Non Linéaire Appliquée Université de Toulon et du Var BP 132-83957 La Garde Cedex seppecher@univ-tln.fr Abstract.
More informationInternational Journal of Pure and Applied Sciences and Technology
Int. J. Pure Appl. Sci. Technol., 3() (0), pp. 7-39 International Journal of Pure and Applied Sciences and Technology ISSN 9-607 Available online at www.ijopaasat.in Research Paper Reflection of Quasi
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationSection 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)).
Difference Equations to Differential Equations Section 8.5 Applications: Pendulums Mass-Spring Systems In this section we will investigate two applications of our work in Section 8.4. First, we will consider
More informationU.S. South America Workshop. Mechanics and Advanced Materials Research and Education. Rio de Janeiro, Brazil. August 2 6, Steven L.
Computational Modeling of Composite and Functionally Graded Materials U.S. South America Workshop Mechanics and Advanced Materials Research and Education Rio de Janeiro, Brazil August 2 6, 2002 Steven
More informationElastic Waves Scattering without Conversion in Metamaterials with. Simultaneous Zero Indices for Longitudinal and Transverse waves
lastic Waves Scattering without Conversion in Metamaterials with Simultaneous Zero Indices for Longitudinal and Transverse waves Fengming Liu,,3 and Zhengyou Liu 2 School of Science, Hubei University of
More informationHot Topics in Physical Acoustics
Hot Topics in Physical Acoustics J.R. (Josh) Gladden Dept. of Physics and Astronomy University of ASA Fall Meeting November 12, 2008 Outline Sound waves in the early universe Nature of sound in a hot plasma
More informationA new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem
A new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem E Hanukah Faculty of Mechanical Engineering, Technion Israel Institute of Technology, Haifa
More informationBand structure and transmission of 3D chiral sonic crystals
Band structure and transmission of 3D chiral sonic crystals R. Pico, V. Romero-Garcia, V. Sanchez-Morcillo, A. Cebrecos and L.M. Garcia-Raffi Universidad Politecnica de Valencia, Paranimf 1, 46730 Gandia,
More information6th NDT in Progress Lamb waves in an anisotropic plate of a single crystal silicon wafer
6th NDT in Progress 2011 International Workshop of NDT Experts, Prague, 10-12 Oct 2011 Lamb waves in an anisotropic plate of a single crystal silicon wafer Young-Kyu PARK 1, Young H. KIM 1 1 Applied Acoustics
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More information