0
Research Papers

# Modulating Elastic Band Gap Structure in Layered Soft Composites Using Sacrificial InterfacesOPEN ACCESS

[+] Author and Article Information
Qianli Chen, Ahmed Elbanna

Department of Civil and Environmental
Engineering,
University of Illinois at Urbana Champaign,
Champaign, IL 61801

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received June 13, 2016; final manuscript received August 23, 2016; published online September 9, 2016. Assoc. Editor: Kyung-Suk Kim.

J. Appl. Mech 83(11), 111009 (Sep 09, 2016) (8 pages) Paper No: JAM-16-1293; doi: 10.1115/1.4034537 History: Received June 13, 2016; Revised August 23, 2016

## Abstract

A wide range of engineered and natural composites exhibit a layered architecture whereby individual building blocks are assembled layer by layer using cohesive interfaces. We present a novel mechanism for evolving acoustic band gap structure in a model system of these composites through patterning the microstructure in a way that triggers nonplanar interfacial deformations between the layers as they are stretched. Through the controlled deformation and growth of interlayer channels under macroscopic tension, we observe the emergence of multiple wide band gaps due to Bragg diffraction and local resonance. We describe these phenomena in details for three example microstructures and discuss the implications of our approach for harnessing controlled deformation in modulating band gap properties of composite materials.

<>

## Introduction

Layered composites are frequently found as building blocks of tough biological structures [14]. Among those, bone is a typical example exhibiting remarkable combination of stiffness and toughness [5]. The basic building block of bone at the microscale is the mineralized collagen fibril (MCF) in which the soft collagen matrix is reinforced by stiff hydrated calcium phosphates [6]. The MCFs are glued together using interfaces of noncollagenous proteins with sacrificial bonds and hidden length to form a layered structure [710]. Another example of a layered structure is seashell nacre which endows dramatic improvement of stiffness and toughness compared to its constituent phases independently [11]. The enhancement of toughness in those biological materials is partially attributed to wavy surfaces and cohesive interaction along the interfaces between bulk materials [1214]. These interfaces are sacrificial in the sense that they represent weak spots where energy gets dissipated by progressive accumulation of damage while the rest of the composite remains intact and elastic. The naturally optimized properties of stiffness, strength, and toughness, suggest that there is a link between geometric patterning, sacrificial interfaces, and improved mechanical performance. More recently, it has also been hypothesized that the origin of these improved properties, in a composite like bone, for example, may be linked to the ability of generating multiple band gaps [15].

Composite materials with periodic microstructure are capable of generating a band gap structure through which elastic and acoustic wave propagation may be controlled [1620]. Waves with frequencies falling into the band gap are barred by the periodic structure. The band gap phenomenon in composite material is caused by to two mechanisms: Bragg diffraction and local resonance [21,22]. The Bragg diffraction occurs when the elastic wavelength is comparable to the periodicity of the microstructure. The position and the width of the band gap may be controlled by careful design of the microstructure materials [2326]. Moreover, mechanically triggered large deformation [27] and instability-induced interfacial wrinkling [28,29] are found of capacity to tune the evolution of band gap profile. Local resonance happens when the phases of composite material have strong elastic properties contrast. Banded frequencies caused by local resonance are generally lower than those caused by Bragg diffraction [30] and may exist even in the absence of periodicity and symmetry [21,31]. Designing such band gap structures enables the application of wave filtering, wave guiding, acoustic mirrors and vibration isolators [3238].

We showed in an earlier paper that a bone-inspired composite, with patterned heterogeneity, exhibits the capacity of generating tunable corrugation under externally applied concentric tension [39]. By stacking 1D fibers, with the designed microstructure, to form layered composites, these corrugations collectively lead to the formation of interlayer channels with shapes and sizes tunable by the level of stretch (Fig. 1). In this paper, we demonstrate an application of harnessing this design in evolving the elastic band gap structure in these composites. A unique feature of our design is the dynamic evolution of the structural composition as a function of stretch and inclusion distribution. In particular, an additional phase, namely voided channels, gradually emerge (disappear) with increasing (decreasing) stretch despite being absent at zero stretch.

## Model and Results

We consider layered composites made up of stacks of fibrils that are glued along their longer dimension as shown in Fig. 1. We choose fibril thickness $h=30 μm$, length of inclusion $Li=75 μm$, inclusion thickness $t=7.5 μm$, inclusion spacing $S=40 μm$, and inclusion eccentricity from center line of fibril $c=4.5μm$ [39]. Three examples of inclusion patterns are discussed: (a) nonstaggered aligned inclusions (Fig. 1(a)), (b) staggered identical inclusions (Fig. 1(b)), and (c) staggered symmetric inclusions (Fig. 1(c)). An extensive analysis of the effect of geometry and material properties on the response of these composites under longitudinal (along x-axis in Fig. 1) stretch are reported in our previous paper [39]. In general, the corrugation amplitude increases with the stiffness of inclusion and its eccentricity from the local tension axis.

The unit cell is discretized using eight-node bi-quadrilateral elements (Q8). We run the basic stretch analysis and band gap calculations in an in-house matlab code for finite-deformation elasticity. The dynamic wave propagation verification simulation is run by the commercial finite-element software abaqus. The two phases are modeled using Neo-Hookean hyperelastic material model. The silicon inclusion has Young's modulus of $Ei=170 GPa$, Poisson's ratio $νi=0.064$, and density $ρi=2300 kg/m3$. The matrix material is modeled as polydimethylsiloxane (PDMS) with initial Young's modulus of $Em=750 kPa$, Poisson's ratio of 0.495, and density $ρm=1000 kg/m3$. This combination of materials were used in our previous experiment [40]. Reduced integration is applied for quadrature and the nonlinear problem is solved with full Newton Raphson method. Elastic band structure calculations are executed for a unit cell in the deformed configuration, and thus, the dimensions of the Brillouin zone are updated accordingly with Craig–Bampton mode decomposition, and Bloch boundary conditions are then applied [41]. The formulas for finite element method (FEM) and block wave calculations are summarized in the Appendix.

The interaction between two composite fibrils is simulated by a traction–separation law with equivalent Young's modulus of $Ecoh=1 kPa$, Poisson's ratio $νcoh=0.495$, and initial thickness $tcoh=1μm$. Then, we have equivalent normal stiffness $Kn=Ecohtop/tcoh$ and shear stiffness $Ks=Ecohtop/[2tcoh(1+νcoh)]$, where $top$ is the out-of-plane thickness of fibril. In this paper, we use the bilinear law (Eq. (A7)) for simplicity but the qualitative nature of the results depends weakly on the specific form and properties of the cohesive law. To investigate the influence of the interfaces, the responses of both monolithic (i.e., no interfaces) and layered systems (with interfaces) under uniaxial stretch along x-axis in Fig. 1 are simulated and compared.

Three monolithic composites with inclusions are modeled and served as reference groups for comparisons. Unlike the layered system there are no interfaces in the monolithic case. The monolithic nature of matrix material stiffens the overall composite response. The Bloch wave analysis reveals that no elastic band-structure behavior is found except in case (c) where a narrow band gap of width 0.008 MHz appears near 0.26 MHz (Fig. 2).

The introduction of interfaces enables more complex deformation patterns and richer elastic band gap structure. In the case of the nonstaggered inclusion pattern (case a), the fibrils shrink in the transverse direction, within the region of pure PDMS, more than in the regions where the silicon inclusions exist. In cases (b) and (c), corrugation develops due to the eccentricity of inclusions. Channels, of different shapes and orientations, form as stretch level increases because of different stacking patterns. All three cases show deformed wavy surface under uniaxial stretch along x-axis in Fig. 1, and this alters the dispersion property of unit cell.

Unlike in the monolithic cases, the layered composites have band gaps at the undeformed state and are able to develop richer behavior after stretch.

###### Case (a) Composite With Nonstaggered Inclusions.

Figure 3(a) reveals that there exist multiple wide elastic band gaps at zero strain. As stretch increases, some of the smaller bandgaps disappear and other bandgaps become even wider (Fig. 4(a)). To understand the origin of the band gaps in this case, it is illuminating to analyze the eigenmode deformation. For that purpose, selected eigenmodes are plotted to the right of the dispersion curves. Most of these eigenmodes suggest that the deformation is mainly concentrated within soft matrix material representing a local resonance behavior. Few exceptions exist. For example modes D, F, and D’ show distributed deformations where the inclusions also vibrate due to Brag scattering. While in mode F’ the deformation is localized in the matrix but the influence of the channel is apparent: the deformation is concentrated close to the channel tips. Furthermore, the dispersion lines along the vertical propagation direction $Y−Γ$ are almost flat signifying very low group velocities. This is due to the existence of the horizontal cohesive interfaces with low stiffness values that slow down wave propagation across them. Finally, eigenmodes of the composite in this inclusion pattern keep their general shape as stretch grows but the corresponding eigenfrequencies change.

###### Case (b) Composite With Identical Staggered Inclusions.

Figure 3(b) shows examples of the band gap structure in this case. Fewer band gaps exist compared to case (a) and they mostly arise due to local resonance in the PDMS matrix. Once again, the effect of channel development is apparent where deformations are banded parallel and normal to their locations. Furthermore, the dispersion relations are altered significantly by deformation, unlike in case (a). This is primarily due to the complex growth of the interfacial channels.

###### Case (c) Composite With Symmetrically Staggered Inclusions.

Figure 3(c) shows examples of the band gap structure in this case. Unlike the previous cases, no band gaps develop at zero strain. Also the dispersion lines along the vertical propagation direction $Y−Γ$ are not as flat as in cases (a) and (b) above. This is attributed to the development of contact between the layers upon bending. The composite is thus continuous through contact region and waves may propagate across the horizontal interfaces at the speed allowed by the matrix material. We also note that the eigenmode analysis reveals that the deformation at the boundaries of the band gap is localized near the nonplanar features of the channels, signifying the influence of the complex evolution of the interfacial separation.

The band gap evolution as a function of stretch is shown in Fig. 4. As previously discussed, multiple band gaps develop as stretch level increases. The composite with nonstaggered inclusions (case (a)), shows multiple wide band gaps above 0.13 MHz at all levels of stretch. For case (b), four clusters of band gaps initially exists near frequencies 0.16, 0.23, 0.33, and 0.35 MHz, respectively. As stretch increases, the band gap near 0.16 MHz becomes wider and moves to higher frequency while the other three close. Near stretch levels of 1.15 and 1.17, two band gaps open at frequency 0.26 and 0.28 MHz, respectively, but they are generally smaller than the lower frequency band gap. For case (c), a band gap appears after stretch of 1.05 near 0.2 MHz and opens up to width of 0.04 MHz as the stretch level increases.

The opening and closure of band gaps as a function of deformation may be explained as follows. Due to the hyperelastic nature of the matrix material, the matrix material undergoes stiffness softening as stretch grows. Stiffness changes nonuniformly among the bulk constituents with stretch since the matrix experiences the most changes while the silicon inclusion remains almost linear elastic with constant stiffness. The variability of the band gap width is due to the competition between stiffness changes as the hyperelastic material is changed and the evolving geometry due to the nonuniform deformation of the interfaces which lead to complex channel shapes and scattering response.

As a verification of band gap behavior, we examine the transmission plots for the composite with staggered inclusion pattern (case (b)) as an example. Nine unit cells are connected together in the longitudinal direction, and periodic boundary condition is applied to the upper and lower boundaries of those cells to create a semifinite composite system. We connect one end of the composite system to infinite elements to absorb waves and apply stretch and vibration at the other end. The finite-element analysis proceeds as follows. We first stretch the composite to the required level. Then, we apply a small amplitude vibration (less than 1% of the stretch amplitude) at the stretched end. We compute the transmission ratio as the ratio of the vibration amplitude at the interface of the composite with the infinite elements to the applied vibration amplitude at the stretched end. Two scenarios are presented in Fig. 5, one at stretch of 1.10 where no band gap has shown up between 0.2 MHz and 0.3 MHz, and one at stretch of 1.30, where there are two band gaps near 0.26 MHz and 0.28 MHz. From Fig. 5(b), we may observe that for both stretch levels, shear waves are attenuated above 0.25 MHz. However, for the longitudinal waves (Fig. 5(a)), the attenuation is most noticeable below 0.24 MHz and above 0.29 MHz for the stretch level of 1.10 but it is coincident with the band gap region at stretch level of 1.30. That is, while longitudinal and shear waves may each attenuate at frequencies outside the band gap, it is only within the band gap that the amplitudes of both types of waves are significantly reduced simultaneously.

## Discussion

We have numerically demonstrated the capability of modulating band gap structure in soft composites using sacrificial interfaces. While several approaches have been used previously to achieve this purpose, including designing of hierarchical composite structure [23,37], triggered instabilities in monolithic periodic material [27,42], and instability-induced interfacial wrinkling [28,29], this is the first time, to the best of our knowledge, that sacrificial interfaces with controlled damage behavior is used to tune the band gap behavior.

The introduction of interfaces to a monolithic composite enriches the deformability of the system and provides a mean to create and alter the band gap structure. The interfacial deformation along composite fibrils is controlled by the inclusion pattern. Nonstaggered (case a) and staggered symmetrical (case c) inclusion patterns result in the formation of parallel horizontal array of holes in aligned and staggered manners, respectively, while the staggered inclusion pattern (case b) leads to oblique orientation of holes. In all cases, the composite geometry is evolving gradually. The layered nature of the composite makes the group velocity of vertical wave approach zero in case (a) and (b). In case (c), the contact between the fibrils promote relatively faster group velocities.

The eigenmodes of unit cell indicate that the band gap evolution in our system is facilitated by both Bragg diffraction and local resonance. The interfacial channels represent a new periodic inclusion of degrading stiffness and irregular boundary geometry. The periodicity of the channels provides additional diffraction paths for wave propagations. Moreover, wave reflection and refraction at the channel/matrix interface is opposite in polarity to wave dynamics at the matrix/stiffer inclusion interface (channel stiffness < matrix stiffness < inclusion stiffness) and this contribute to local resonance.

In the current study, the composite is elastic throughout the deformation, and damage is only localized to the interfaces between fibrils. With the use of self-healing interfaces [43,44], the phenomena we are describing here are completely reversible. Furthermore, even if the interface remains damaged, the channels nucleate and grow to complex irregular shapes that interfere with the wave propagation and cause more complex behavior and possibly more frequency band gaps than composites with regular holes or cuts [42,45,46].

In general, band gaps of the composite fibril system shows rich and controllable band gap behavior than the monolithic composite where no band gap shows up even with the same inclusion pattern as the layered composite (except for a tiny band gap in case (c)). It is the interaction between the inclusion pattern and the sacrificial interfaces which manipulate wave diffraction and interference leading to the creation of multiple frequency band gaps. For the cases considered in this paper, we show that the band gaps exist at initial stretch-free state and can be modulated through stretch level. Multiple wide band gaps are observed in nonstaggered inclusion pattern (case a) making it very efficient for wave attenuation. For the other two cases, band gap are also wide but not as rich as the case (a). For all three cases, the introduction of interfaces and deformation-caused stiffness redistribution plays essential roles in the evolution of these band gaps.

Future extension of this study will involve investigation of further applications and optimization of the tunable band gap structure. For example, the emergence of multiple band gaps under large stretch may have applications in wave filtering, sound isolation, and vibration damping. Furthermore, the width of the band gaps may be maximized through topology optimization so that the tunable band gaps, observed in the current study, may extend over wider regions of frequencies. Finally, the connection between band gap structure and the effective toughness of the composite material is a topic that requires further exploration. In particular, it may be possible, through careful designing of the inclusion pattern and sacrificial interfaces, to be able to control crack nucleation and propagation by varying the band gap structure in the composite as a function of deformation.

## Acknowledgements

This work was supported by the National Science Foundation through Grant No. CMMI-1435920 and UIUC Start-up funds.

## Appendices

###### Appendix: Constitutive Laws

The material deformation is described as a mapping from reference configuration $X$ to the deformed configuration $x$Display Formula

(A1)$x=χ(X,t)$

The displacement of a material point form is Display Formula

(A2)$u=x−X=χ(X,t)−X$

The deformation gradient tensor is defined as Display Formula

(A3)$F=∂x∂X=I+∂u∂X$

where $I$ is the identity matrix, and its Jacobian is $J=det(F)$.

Two different soft materials are building up the composite (Fig. 1). Each material is assumed to be homogeneous and isotropic with hyperelastic response. We use the generalized Neo-Hookean solid [47,48] model with strain energy density proposed by Display Formula

(A4)$W=μ2(I¯1−3)+K2(J−1)2$

where $I¯1=I1/J2/3$ making it more convenient for nearly incompressible materials, $μ$ and $K$ are the shear modulus and bulk modulus of the material, respectively.

The stress–strain relation and tangent stiffness follow as Display Formula

(A5)$σ=μJ5/3(B−13I1I)+K(J−1)I$

(A6)$Cijkl=μJ2/3[δikBjl+δjkBil−23(δklBij+δijBkl)+23I13δijδkl]+K(2J−1)Jδijδkl$

###### Cohesive Interaction for the Layered Composite.

The interaction between two composite fibrils is simulated by a traction–separation relation. In this paper, we use the bilinear law for simplicity but the qualitative nature of the results depends weakly on the specific form of the cohesive law. The two-dimensional bilinear cohesive law is defined by [49,50] Display Formula

(A7)$Tn,s={Kn,sδn,s, δ≤δcKn,sδcδ−δfδc−δf, δc≤δ≤δf 0, δ≥δf$

where $Tn,s$, $Kn,s$, $δn,s$ is the normal or tangent traction, stiffness, and separation, respectively. $δ$ is the magnitude of total separation given by $δ=δn2+δs2$. $δc$ and $δf$ determine the damage initialization and complete failure, respectively (Fig. 6).

Figure 7 shows the effect of varying the cohesive interface equivalent Young's modulus on the band gap profile of the composite with the nonstaggered inclusion pattern (Case a) at a stretch-free state. In this case of zero elongation, since the cohesive material remains intact along the whole interface, the band gap profile is only affected by the elastic normal and shear stiffness of the layer and not its softening response. The results in Fig. 7 suggest that the band gap width decreases as the equivalent Young's modulus of the cohesive interface increases The dispersion lines plot converges to the monolithic composite case (Fig. 2(a)) as the modulus increases. Indeed, when $Ecoh$ is infinitely large (which is equivalent to imposing a rigid constraint between the two points on the opposite side of the interface), we will get exactly the same result in Fig. 2(a) where interface no longer exists.

###### Wave Propagation.

The equations of motion are given by Display Formula

(A8)$∂σ∂x=ρ∂2u∂t2$

in absence of body force, where $σ$ is the Cauchy stress, and $ρ$ is the density which varies for different material in composite.

The propagation of elastic wave in periodic composites may be analyzed by Bloch's theorem [51] and further developed in the following way [52]. We start by decomposing the solution as Display Formula

(A9)$u(X,t)=U(X)e−iwt$

where $U(X)$ is a complex valued function of $X$.

Accordingly, we also have the following decomposition for stress Display Formula

(A10)$σ(X,t)=Σ(X)e−iwt$

Equation (A8) now becomes Display Formula

(A11)$∂Σ∂x=ρω2U$

Due to periodicity, the displacement is constrained by Display Formula

(A12)$U(X)=U(X+L)e−ikL$

where $L$ is the vector connecting equivalent points in periodic structure, $k$ is the wave vector contain the information of wavelength and direction.

Similarly, the stress also satisfies Display Formula

(A13)$Σ(X)=Σ(X+L)e−ikL$

which is satisfied by applying the constrain of displacement under FEM discretization.

The wave vector is chosen from the first Brillouin zone in reciprocal lattice for the representative unit cell as shown in Fig. 7. The reciprocal lattice is defined by base vectors $R1$ and $R2$ and Display Formula

(A14)$R1=2πA2×ea, R2=2πe×A1a$

where $A1$ and $A2$ are base vector in direct lattice (Fig. 1), $a=‖A1×A2‖$, and $e=(A1×A2)/a$ so that Display Formula

(A15)$RiAj=2πδij$

It is suggested in the literature that the band gap information may be obtained by surfing only the points along boundaries of this zone [53]. While this is not rigorously proven, we adopt the same procedure in the current paper (Fig 8).

The wave propagation response is analyzed by solving the eigenvalue problem of Display Formula

(A16)$K−ω2M=0$

The matrices $K$ and $M$ are the tangent stiffness and mass matrices generated from the standard finite-element nonlinear analysis at different stretch level. The displacement constrain of Eq. (A12) should be applied to both $K$ and $M$ before solving the eigenvalue problem. Mode decomposition can then be used to accelerate the computation process [41].

## References

Espinosa, H. D. , Rim, J. E. , Barthelat, F. , and Buehler, M. J. , 2009, “ Merger of Structure and Material in Nacre and Bone—Perspectives on de Novo Biomimetic Materials,” Prog. Mater. Sci., 54(8), pp. 1059–1100.
Luz, G. M. , and Mano, J. F. , 2010, “ Mineralized Structures in Nature: Examples and Inspirations for the Design of New Composite Materials and Biomaterials,” Compos. Sci. Technol., 70(13), pp. 1777–1788.
Dunlop, J. W. C. , and Fratzl, P. , 2010, “ Biological Composites,” Annu. Rev. Mater. Res., 40(1), pp. 1–24.
Ritchie, R. O. , 2011, “ The Conflicts Between Strength and Toughness,” Nat. Mater., 10(11), pp. 817–822. [PubMed]
Weiner, S. , and Wagner, H. D. , 1998, “ The Material Bone: Structure-Mechanical Function Relations,” Annu. Rev. Mater. Sci., 28(1), pp. 271–298.
Jäger, I. , and Fratzl, P. , 2000, “ Mineralized Collagen Fibrils: A Mechanical Model With a Staggered Arrangement of Mineral Particles,” Biophys. J., 79(4), pp. 1737–1746. [PubMed]
Elbanna, A. E. , and Carlson, J. M. , 2013, “ Dynamics of Polymer Molecules With Sacrificial Bond and Hidden Length Systems: Towards a Physically-Based Mesoscopic Constitutive Law,” PLoS One, 8(4), p. e56118. [PubMed]
Lieou, C. K. C. , Elbanna, A. E. , and Carlson, J. M. , 2013, “ Sacrificial Bonds and Hidden Length in Biomaterials: A Kinetic Constitutive Description of Strength and Toughness in Bone,” Phys. Rev. E: Stat. Nonlinear, Soft Matter Phys., 88, p. 012703.
Wang, W. , and Elbanna, A. , 2014, “ Crack Propagation in Bone on the Scale of Mineralized Collagen Fibrils: Role of Polymers With Sacrificial Bonds and Hidden Length,” Bone, 68, pp. 20–31. [PubMed]
Fantner, G. E. , Hassenkam, T. , Kindt, J. H. , Weaver, J. C. , Birkedal, H. , Pechenik, L. , Cutronni, J. A. , Cidade, G. A. G. , Stucky, G. D. , Morse, D. E. , and Hansma, P. K. , 2005, “ Sacrificial Bonds and Hidden Length Dissipate Energy as Mineralized Fibrils Separate During Bone Fracture,” Nat. Mater., 4(8), pp. 612–616. [PubMed]
Zhang, N. , Yang, S. , Xiong, L. , Hong, Y. , and Chen, Y. , 2016, “ Nanoscale Toughening Mechanism of Nacre Tablet,” J. Mech. Behav. Biomed. Mater., 53, pp. 200–209. [PubMed]
Li, B.-W. , Zhao, H.-P. , Qin, Q.-H. , Feng, X.-Q. , and Yu, S.-W. , 2012, “ Numerical Study on the Effects of Hierarchical Wavy Interface Morphology on Fracture Toughness,” Comput. Mater. Sci., 57, pp. 14–22.
Chan, E. P. , Smith, E. J. , Hayward, R. C. , and Crosby, A. J. , 2008, “ Surface Wrinkles for Smart Adhesion,” Adv. Mater., 20(4), pp. 711–716.
Janarthanan, V. , Garrett, P. D. , Stein, R. S. , and Srinivasarao, M. , 1997, “ Adhesion Enhancement in Immiscible Polymer Bilayer Using Oriented Macroscopic Roughness,” Polymer (Guildf), 38(1), pp. 105–111.
Davies, B. , King, A. , Newman, P. , Minett, A. , Dunstan, C. R. , and Zreiqat, H. , 2014, “ Hypothesis: Bones Toughness Arises From the Supression of Elastic Waves,” Sci. Rep., 4, p. 7538. [PubMed]
Maldovan, M. , and Thomas, E. L. , 2008, Periodic Materials and Interference Lithography, Wiley-VCH Verlag GmbH, New York.
Lee, J.-H. , Singer, J. P. , and Thomas, E. L. , 2012, “ Micro-/Nanostructured Mechanical Metamaterials,” Adv. Mater., 24(36), pp. 4782–810. [PubMed]
Aly, A. H. , and Mehaney, A. , 2012, “ Enhancement of Phononic Band Gaps in Ternary/Binary Structure,” Phys. B Condens. Matter, 407(21), pp. 4262–4268.
Maldovan, M. , 2013, “ Sound and Heat Revolutions in Phononics,” Nature, 503(7475), pp. 209–217. [PubMed]
Sigalas, M. M. , and Economou, E. N. , 1995, “ Comment on Acoustic Band Structure of Periodic Elastic Composites,” Phys. Rev. Lett., 75(19), p. 3580. [PubMed]
Goffaux, C. , and Sánchez-Dehesa, J. , 2003, “ Two-Dimensional Phononic Crystals Studied Using a Variational Method: Application to Lattices of Locally Resonant Materials,” Phys. Rev. B, 67(14), p. 144301.
Achaoui, Y. , Khelif, A. , Benchabane, S. , Robert, L. , and Laude, V. , 2011, “ Experimental Observation of Locally-Resonant and Bragg Band Gaps for Surface Guided Waves in a Phononic Crystal of Pillars,” Phys. Rev. B, 83(10), p. 104201.
Chen, Y. , and Wang, L. , 2014, “ Tunable Band Gaps in Bio-Inspired Periodic Composites With Nacre-Like Microstructure,” J. Appl. Phys., 116(6), p. 63506.
Goffaux, C. , and Vigneron, J. P. , 2001, “ Theoretical Study of a Tunable Phononic Band Gap System,” Phys. Rev. B, 64(7), p. 75118.
Vatanabe, S. L. , Paulino, G. H. , and Silva, E. C. N. , 2014, “ Maximizing Phononic Band Gaps in Piezocomposite Materials by Means of Topology Optimization,” J. Acoust. Soc. Am., 136(2), pp. 494–501. [PubMed]
Sigmund, O. , and Jensen, J. S. , 2003, “ Systematic Design of Phononic Band-Gap Materials and Structures by Topology Optimization,” Philos. Trans. A. Math. Phys. Eng. Sci., 361(1806), pp. 1001–1019. [PubMed]
Bertoldi, K. , and Boyce, M. C. , 2008, “ Wave Propagation and Instabilities in Monolithic and Periodically Structured Elastomeric Materials Undergoing Large Deformations,” Phys. Rev. B: Condens. Matter Mater. Phys., 78(18), p. 184107.
Rudykh, S. , and Boyce, M. C. , 2014, “ Transforming Wave Propagation in Layered Media Via Instability-Induced Interfacial Wrinkling,” Phys. Rev. Lett., 112, pp. 1–5.
Bayat, A. , and Gordaninejad, F. , 2015, “ Switching Band-Gaps of a Phononic Crystal Slab by Surface Instability,” Smart Mater. Struct., 24(7), p. 75009.
Liu, Z. , 2000, “ Locally Resonant Sonic Materials,” Science, 289(5485), pp. 1734–1736. [PubMed]
Wang, G. , Wen, X. , Wen, J. , Shao, L. , and Liu, Y. , 2004, “ Two-Dimensional Locally Resonant Phononic Crystals With Binary Structures,” Phys. Rev. Lett., 93(15), p. 154302. [PubMed]
Wang, X.-P. , Jiang, P. , Chen, T.-N. , and Zhu, J. , 2015, “ Tuning Characteristic of Band Gap and Waveguide in a Multi-Stub Locally Resonant Phononic Crystal Plate,” AIP Adv., 5(10), p. 107141.
Kushwaha, M. S. , 1997, “ Stop-Bands for Periodic Metallic Rods: Sculptures That can Filter the Noise,” Appl. Phys. Lett., 70(24), pp. 3218–3220.
Khelif, A. , Deymier, P. A. , Djafari-Rouhani, B. , Vasseur, J. O. , and Dobrzynski, L. , 2003, “ Two-Dimensional Phononic Crystal With Tunable Narrow Pass Band: Application to a Waveguide With Selective Frequency,” J. Appl. Phys., 94(3), p. 1308.
Vasseur, J. O. , Deymier, P. A. , Khelif, A. , Lambin, P. , Djafari-Rouhani, B. , Akjouj, A. , Dobrzynski, L. , Fettouhi, N. , and Zemmouri, J. , 2002, “ Phononic Crystal With Low Filling Fraction and Absolute Acoustic Band Gap in the Audible Frequency Range: A Theoretical and Experimental Study,” Phys. Rev. E. Stat. Nonlinear Soft Matter Phys., 65(5), p. 56608.
Kushwaha, M. S. , Halevi, P. , Dobrzynski, L. , and Djafari-Rouhani, B. , 1993, “ Acoustic Band Structure of Periodic Elastic Composites,” Phys. Rev. Lett., 71(13), pp. 2022–2025. [PubMed]
Chen, Y. , and Wang, L. , 2015, “ Multiband Wave Filtering and Waveguiding in Bio-Inspired Hierarchical Composites,” Extrem. Mech. Lett., 5, pp. 18–24.
Khelif, A. , Choujaa, A. , Benchabane, S. , Djafari-Rouhani, B. , and Laude, V. , 2004, “ Guiding and Bending of Acoustic Waves in Highly Confined Phononic Crystal Waveguides,” Appl. Phys. Lett., 84(22), p. 4400.
Chen, Q. , and Elbanna, A. , 2015, “ Tension-Induced Tunable Corrugation in Two-Phase Soft Composites: Mechanisms and Implications,” Extrem. Mech. Lett., 4, pp. 26–37.
Yang, Z. , Chen, Q. , Elbanna, A. E. , and Kim, S. , 2016, “ Transfer Printing Enabled Soft Composite Films for Tunable Surface Topography,” Extrem. Mech. Lett., 7, pp. 145–153.
Krattiger, D. , and Hussein, M. I. , 2014, “ Bloch Mode Synthesis: Ultrafast Methodology for Elastic Band-Structure Calculations,” Phys. Rev. E—Stat. Nonlinear, Soft Matter Phys., 90, p. 063306.
Bertoldi, K. , and Boyce, M. C. , 2008, “ Mechanically Triggered Transformations of Phononic Band Gaps in Periodic Elastomeric Structures,” Phys. Rev. B, 77(5), p. 52105.
Shen, Z. L. , Dodge, M. R. , Kahn, H. , Ballarini, R. , and Eppell, S. J. , 2008, “ Stress-Strain Experiments on Individual Collagen Fibrils,” Biophys. J., 95(8), pp. 3956–3963. [PubMed]
Ahn, B. K. , Lee, D. W. , Israelachvili, J. N. , and Waite, J. H. , 2014, “ Surface-Initiated Self-Healing of Polymers in Aqueous Media,” Nat. Mater., 13(9), pp. 867–872. [PubMed]
Vasseur, J. O. , Deymier, P. A. , Djafari-Rouhani, B. , Pennec, Y. , and Hladky-Hennion, A.-C. , 2008, “ Absolute Forbidden Bands and Waveguiding in Two-Dimensional Phononic Crystal Plates,” Phys. Rev. B, 77(8), p. 85415.
Tang, Y. , Lin, G. , Han, L. , Qiu, S. , Yang, S. , and Yin, J. , 2015, “ Design of Hierarchically Cut Hinges for Highly Stretchable and Reconfigurable Metamaterials With Enhanced Strength,” Adv. Mater., 27(44), pp. 7181–7190. [PubMed]
Treloar, L. R. G. , 1948, “ Stresses and Birefringence in Rubber subjected to General Homogeneous Strain,” Proc. Phys. Soc., 60(2), pp. 135–144.
Bower, A. F. , “ Constitutive Laws—3.5 Hyperelasticity,” Applied Mechanics of Solids, CRC Press, Boca Raton, FL, Chap. 3.
Swaminathan, S. , Pagano, N. J. , and Ghosh, S. , 2006, “ Analysis of Interfacial Debonding in Three-Dimensional Composite Microstructures,” ASME J. Eng. Mater. Technol., 128(1), pp. 96–106.
Lee, S. C. , and Kim, K. T. , 2008, “ A Densification Model for Powder Materials Under Cold Isostatic Pressing—Effect of Adhesion and Friction of Rubber Molds,” Mater. Sci. Eng. A, 498(1–2), pp. 359–368.
Brillouin, L. , 1953, Wave Propagation in Periodic Structures; Electric Filters and Crystal Lattices, Dover Publications, Mineola, NY.
Aberg, M. , 1997, “ The Usage of Standard Finite Element Codes for Computation of Dispersion Relations in Materials With Periodic Microstructure,” J. Acoust. Soc. Am., 102(4), p. 2007.
Hussein, M. I. , 2009, “ Reduced Bloch Mode Expansion for Periodic Media Band Structure Calculations,” Proc. Math. Phys. Eng. Sci., 465(2109), pp. 2825–2848.
View article in PDF format.

## References

Espinosa, H. D. , Rim, J. E. , Barthelat, F. , and Buehler, M. J. , 2009, “ Merger of Structure and Material in Nacre and Bone—Perspectives on de Novo Biomimetic Materials,” Prog. Mater. Sci., 54(8), pp. 1059–1100.
Luz, G. M. , and Mano, J. F. , 2010, “ Mineralized Structures in Nature: Examples and Inspirations for the Design of New Composite Materials and Biomaterials,” Compos. Sci. Technol., 70(13), pp. 1777–1788.
Dunlop, J. W. C. , and Fratzl, P. , 2010, “ Biological Composites,” Annu. Rev. Mater. Res., 40(1), pp. 1–24.
Ritchie, R. O. , 2011, “ The Conflicts Between Strength and Toughness,” Nat. Mater., 10(11), pp. 817–822. [PubMed]
Weiner, S. , and Wagner, H. D. , 1998, “ The Material Bone: Structure-Mechanical Function Relations,” Annu. Rev. Mater. Sci., 28(1), pp. 271–298.
Jäger, I. , and Fratzl, P. , 2000, “ Mineralized Collagen Fibrils: A Mechanical Model With a Staggered Arrangement of Mineral Particles,” Biophys. J., 79(4), pp. 1737–1746. [PubMed]
Elbanna, A. E. , and Carlson, J. M. , 2013, “ Dynamics of Polymer Molecules With Sacrificial Bond and Hidden Length Systems: Towards a Physically-Based Mesoscopic Constitutive Law,” PLoS One, 8(4), p. e56118. [PubMed]
Lieou, C. K. C. , Elbanna, A. E. , and Carlson, J. M. , 2013, “ Sacrificial Bonds and Hidden Length in Biomaterials: A Kinetic Constitutive Description of Strength and Toughness in Bone,” Phys. Rev. E: Stat. Nonlinear, Soft Matter Phys., 88, p. 012703.
Wang, W. , and Elbanna, A. , 2014, “ Crack Propagation in Bone on the Scale of Mineralized Collagen Fibrils: Role of Polymers With Sacrificial Bonds and Hidden Length,” Bone, 68, pp. 20–31. [PubMed]
Fantner, G. E. , Hassenkam, T. , Kindt, J. H. , Weaver, J. C. , Birkedal, H. , Pechenik, L. , Cutronni, J. A. , Cidade, G. A. G. , Stucky, G. D. , Morse, D. E. , and Hansma, P. K. , 2005, “ Sacrificial Bonds and Hidden Length Dissipate Energy as Mineralized Fibrils Separate During Bone Fracture,” Nat. Mater., 4(8), pp. 612–616. [PubMed]
Zhang, N. , Yang, S. , Xiong, L. , Hong, Y. , and Chen, Y. , 2016, “ Nanoscale Toughening Mechanism of Nacre Tablet,” J. Mech. Behav. Biomed. Mater., 53, pp. 200–209. [PubMed]
Li, B.-W. , Zhao, H.-P. , Qin, Q.-H. , Feng, X.-Q. , and Yu, S.-W. , 2012, “ Numerical Study on the Effects of Hierarchical Wavy Interface Morphology on Fracture Toughness,” Comput. Mater. Sci., 57, pp. 14–22.
Chan, E. P. , Smith, E. J. , Hayward, R. C. , and Crosby, A. J. , 2008, “ Surface Wrinkles for Smart Adhesion,” Adv. Mater., 20(4), pp. 711–716.
Janarthanan, V. , Garrett, P. D. , Stein, R. S. , and Srinivasarao, M. , 1997, “ Adhesion Enhancement in Immiscible Polymer Bilayer Using Oriented Macroscopic Roughness,” Polymer (Guildf), 38(1), pp. 105–111.
Davies, B. , King, A. , Newman, P. , Minett, A. , Dunstan, C. R. , and Zreiqat, H. , 2014, “ Hypothesis: Bones Toughness Arises From the Supression of Elastic Waves,” Sci. Rep., 4, p. 7538. [PubMed]
Maldovan, M. , and Thomas, E. L. , 2008, Periodic Materials and Interference Lithography, Wiley-VCH Verlag GmbH, New York.
Lee, J.-H. , Singer, J. P. , and Thomas, E. L. , 2012, “ Micro-/Nanostructured Mechanical Metamaterials,” Adv. Mater., 24(36), pp. 4782–810. [PubMed]
Aly, A. H. , and Mehaney, A. , 2012, “ Enhancement of Phononic Band Gaps in Ternary/Binary Structure,” Phys. B Condens. Matter, 407(21), pp. 4262–4268.
Maldovan, M. , 2013, “ Sound and Heat Revolutions in Phononics,” Nature, 503(7475), pp. 209–217. [PubMed]
Sigalas, M. M. , and Economou, E. N. , 1995, “ Comment on Acoustic Band Structure of Periodic Elastic Composites,” Phys. Rev. Lett., 75(19), p. 3580. [PubMed]
Goffaux, C. , and Sánchez-Dehesa, J. , 2003, “ Two-Dimensional Phononic Crystals Studied Using a Variational Method: Application to Lattices of Locally Resonant Materials,” Phys. Rev. B, 67(14), p. 144301.
Achaoui, Y. , Khelif, A. , Benchabane, S. , Robert, L. , and Laude, V. , 2011, “ Experimental Observation of Locally-Resonant and Bragg Band Gaps for Surface Guided Waves in a Phononic Crystal of Pillars,” Phys. Rev. B, 83(10), p. 104201.
Chen, Y. , and Wang, L. , 2014, “ Tunable Band Gaps in Bio-Inspired Periodic Composites With Nacre-Like Microstructure,” J. Appl. Phys., 116(6), p. 63506.
Goffaux, C. , and Vigneron, J. P. , 2001, “ Theoretical Study of a Tunable Phononic Band Gap System,” Phys. Rev. B, 64(7), p. 75118.
Vatanabe, S. L. , Paulino, G. H. , and Silva, E. C. N. , 2014, “ Maximizing Phononic Band Gaps in Piezocomposite Materials by Means of Topology Optimization,” J. Acoust. Soc. Am., 136(2), pp. 494–501. [PubMed]
Sigmund, O. , and Jensen, J. S. , 2003, “ Systematic Design of Phononic Band-Gap Materials and Structures by Topology Optimization,” Philos. Trans. A. Math. Phys. Eng. Sci., 361(1806), pp. 1001–1019. [PubMed]
Bertoldi, K. , and Boyce, M. C. , 2008, “ Wave Propagation and Instabilities in Monolithic and Periodically Structured Elastomeric Materials Undergoing Large Deformations,” Phys. Rev. B: Condens. Matter Mater. Phys., 78(18), p. 184107.
Rudykh, S. , and Boyce, M. C. , 2014, “ Transforming Wave Propagation in Layered Media Via Instability-Induced Interfacial Wrinkling,” Phys. Rev. Lett., 112, pp. 1–5.
Bayat, A. , and Gordaninejad, F. , 2015, “ Switching Band-Gaps of a Phononic Crystal Slab by Surface Instability,” Smart Mater. Struct., 24(7), p. 75009.
Liu, Z. , 2000, “ Locally Resonant Sonic Materials,” Science, 289(5485), pp. 1734–1736. [PubMed]
Wang, G. , Wen, X. , Wen, J. , Shao, L. , and Liu, Y. , 2004, “ Two-Dimensional Locally Resonant Phononic Crystals With Binary Structures,” Phys. Rev. Lett., 93(15), p. 154302. [PubMed]
Wang, X.-P. , Jiang, P. , Chen, T.-N. , and Zhu, J. , 2015, “ Tuning Characteristic of Band Gap and Waveguide in a Multi-Stub Locally Resonant Phononic Crystal Plate,” AIP Adv., 5(10), p. 107141.
Kushwaha, M. S. , 1997, “ Stop-Bands for Periodic Metallic Rods: Sculptures That can Filter the Noise,” Appl. Phys. Lett., 70(24), pp. 3218–3220.
Khelif, A. , Deymier, P. A. , Djafari-Rouhani, B. , Vasseur, J. O. , and Dobrzynski, L. , 2003, “ Two-Dimensional Phononic Crystal With Tunable Narrow Pass Band: Application to a Waveguide With Selective Frequency,” J. Appl. Phys., 94(3), p. 1308.
Vasseur, J. O. , Deymier, P. A. , Khelif, A. , Lambin, P. , Djafari-Rouhani, B. , Akjouj, A. , Dobrzynski, L. , Fettouhi, N. , and Zemmouri, J. , 2002, “ Phononic Crystal With Low Filling Fraction and Absolute Acoustic Band Gap in the Audible Frequency Range: A Theoretical and Experimental Study,” Phys. Rev. E. Stat. Nonlinear Soft Matter Phys., 65(5), p. 56608.
Kushwaha, M. S. , Halevi, P. , Dobrzynski, L. , and Djafari-Rouhani, B. , 1993, “ Acoustic Band Structure of Periodic Elastic Composites,” Phys. Rev. Lett., 71(13), pp. 2022–2025. [PubMed]
Chen, Y. , and Wang, L. , 2015, “ Multiband Wave Filtering and Waveguiding in Bio-Inspired Hierarchical Composites,” Extrem. Mech. Lett., 5, pp. 18–24.
Khelif, A. , Choujaa, A. , Benchabane, S. , Djafari-Rouhani, B. , and Laude, V. , 2004, “ Guiding and Bending of Acoustic Waves in Highly Confined Phononic Crystal Waveguides,” Appl. Phys. Lett., 84(22), p. 4400.
Chen, Q. , and Elbanna, A. , 2015, “ Tension-Induced Tunable Corrugation in Two-Phase Soft Composites: Mechanisms and Implications,” Extrem. Mech. Lett., 4, pp. 26–37.
Yang, Z. , Chen, Q. , Elbanna, A. E. , and Kim, S. , 2016, “ Transfer Printing Enabled Soft Composite Films for Tunable Surface Topography,” Extrem. Mech. Lett., 7, pp. 145–153.
Krattiger, D. , and Hussein, M. I. , 2014, “ Bloch Mode Synthesis: Ultrafast Methodology for Elastic Band-Structure Calculations,” Phys. Rev. E—Stat. Nonlinear, Soft Matter Phys., 90, p. 063306.
Bertoldi, K. , and Boyce, M. C. , 2008, “ Mechanically Triggered Transformations of Phononic Band Gaps in Periodic Elastomeric Structures,” Phys. Rev. B, 77(5), p. 52105.
Shen, Z. L. , Dodge, M. R. , Kahn, H. , Ballarini, R. , and Eppell, S. J. , 2008, “ Stress-Strain Experiments on Individual Collagen Fibrils,” Biophys. J., 95(8), pp. 3956–3963. [PubMed]
Ahn, B. K. , Lee, D. W. , Israelachvili, J. N. , and Waite, J. H. , 2014, “ Surface-Initiated Self-Healing of Polymers in Aqueous Media,” Nat. Mater., 13(9), pp. 867–872. [PubMed]
Vasseur, J. O. , Deymier, P. A. , Djafari-Rouhani, B. , Pennec, Y. , and Hladky-Hennion, A.-C. , 2008, “ Absolute Forbidden Bands and Waveguiding in Two-Dimensional Phononic Crystal Plates,” Phys. Rev. B, 77(8), p. 85415.
Tang, Y. , Lin, G. , Han, L. , Qiu, S. , Yang, S. , and Yin, J. , 2015, “ Design of Hierarchically Cut Hinges for Highly Stretchable and Reconfigurable Metamaterials With Enhanced Strength,” Adv. Mater., 27(44), pp. 7181–7190. [PubMed]
Treloar, L. R. G. , 1948, “ Stresses and Birefringence in Rubber subjected to General Homogeneous Strain,” Proc. Phys. Soc., 60(2), pp. 135–144.
Bower, A. F. , “ Constitutive Laws—3.5 Hyperelasticity,” Applied Mechanics of Solids, CRC Press, Boca Raton, FL, Chap. 3.
Swaminathan, S. , Pagano, N. J. , and Ghosh, S. , 2006, “ Analysis of Interfacial Debonding in Three-Dimensional Composite Microstructures,” ASME J. Eng. Mater. Technol., 128(1), pp. 96–106.
Lee, S. C. , and Kim, K. T. , 2008, “ A Densification Model for Powder Materials Under Cold Isostatic Pressing—Effect of Adhesion and Friction of Rubber Molds,” Mater. Sci. Eng. A, 498(1–2), pp. 359–368.
Brillouin, L. , 1953, Wave Propagation in Periodic Structures; Electric Filters and Crystal Lattices, Dover Publications, Mineola, NY.
Aberg, M. , 1997, “ The Usage of Standard Finite Element Codes for Computation of Dispersion Relations in Materials With Periodic Microstructure,” J. Acoust. Soc. Am., 102(4), p. 2007.
Hussein, M. I. , 2009, “ Reduced Bloch Mode Expansion for Periodic Media Band Structure Calculations,” Proc. Math. Phys. Eng. Sci., 465(2109), pp. 2825–2848.

## Figures

Fig. 2

Dispersion plots for three monolithic composites corresponding to the three inclusion patterns illustrated in Fig. 1 ((a) nonstaggered, (b) staggered identical, and (c) staggered symmetric) and corresponding Brillouin zones

Fig. 3

Dispersion relations (blue lines), band gap (gray region), and critical eigen-modes of the three composites fibril systems corresponding to the three inclusion patterns ((a) nonstaggered, (b) staggered identical, and (c) staggered symmetric) at stretch level of 1.00/1.01 and 1.20

Fig. 4

Band gap evolutions corresponding to three inclusion patterns ((a) nonstaggered, (b) staggered identical, and (c) staggered symmetric)

Fig. 5

Transmission ratio of composite fibril system corresponding to the case (b) of staggered identical inclusion pattern under: (a) longitudinal wave and (b) shear wave

Fig. 6

Bilinear cohesive law of inter-fibril interaction

Fig. 7

Band gap evolution and dispersion plots corresponding to inclusion patterns (a) of nonstaggered pattern in terms of cohesive layer equivalent Young's modulus

Fig. 8

First Brillouin zone of periodic composite fibril unit cell

Fig. 1

Representative assembly of (a) composite with nonstaggered aligned inclusion pattern, (b) composite with staggered identical inclusion pattern, and (c) composite with staggered symmetric inclusion pattern (cohesive interfaces do not exist in monolithic models). The corresponding unit cells are shown to the right, while the deformed shapes of the composite under uniaxial stretch along x-axis in Fig. 1 are shown below the undeformed ones. Overall dimensions are marked in the unit cell in Fig. 1(a). Inclusion dimensions and spacing are the same in all cases.

## Discussions

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections